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PREFACE 


During my years working as a chemical engineer in development laboratories, pro¬ 
cess engineering groups, and plant startup and support operations, the most frequently 
referenced documents were process flow diagrams (PFDs), which contain the material 
and energy balances and the basic process design information. Equally important were 
process and instrument diagrams (P&IDs), which contain details of all equipment, all 
controls, all instruments, and all lines (i.e., process, instrument, and utilities). Process 
simulation software is an excellent tool for producing high-quality PFDs, and when 
integrated with computer-aided design software, facilitates the production of P&IDs. 
There are several process simulation software systems available to the chemical engi¬ 
neering community, and Aspen Plus is arguably the most popular. 

Teach Yourself the Basics of Aspen Plus™ evolved from two graduate courses that I 
taught at Stevens Institute of Technology over the past 20-odd years. The first course, 
ChE662, is an introduction to steady-state chemical process simulation, which is usually 
taken by graduate students and is organized around a series of workshops that intro¬ 
duce Aspen Plus functionality. Occasionally, undergraduates are enrolled and typically 
experience difficulties in the thermodynamics of phase equilibrium and parameter esti¬ 
mation, due to limitations in their undergraduate courses. The second course, ChE612, 
deals with the analysis and design of complex equilibrium stage processes and with 
difficult multicolumn problems such as, extractive distillation systems. Over time, the 
course evolved from the use of stand-alone two- and three-phase flashes, decantation, 
and two-phase distillation software, to their equivalent blocks in Flowtran and later, 
Aspen Plus. 

The idea for this book originated from my observations of students in these courses. 
I noted that after an initial period dedicated to learning the basics of how to navigate, 
locate material, and enter data into Aspen Plus, students could proceed through the 
exercises, within the workshops, mostly on their own. I would give an introductory 
lecture for each subject studied, show examples, and provide individual help on the 

xiii 
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workshops when needed. It was a rare student who did not finish all of the workshops 
during the course. The book is organized in the same manner, 

If you expect to teach yourself Aspen Plus by reading this book, you will be disap¬ 
pointed. Aspen Plus is a complex process simulator, and in my opinion, the best way to 
learn is with hands-on experience, by attempting each exercise within each workshop, 
and when difficulties are encountered, referring to the problem setup and solutions on 
the accompanying CD as well as the notes in the workshop section at the end of each 
chapter. 

The accompanying CD contains the input and solutions to all the examples and 
workshops. There is a root folder for each chapter, within which there are subfolders 
named Examples and Workshops. Each example and workshop exercise is provided 
in .bkp, Aspen Plus format and .txt format. The .bkp files are set up as input files to 
view details, and may be executed. The .txt files are solutions and may be viewed with 
Notepad by a reader who does not have access to Aspen Plus. References to material 
on the CD within each chapter of the book are by subfolder/filename: for example, 
Chapter Four Examples/Rubin. Some of the workshops were developed using earlier 
versions of Aspen Plus, and when attempting to use them from the CD, a message to 
that effect may appear; however, all have been executed successfully with version 7.0, 
which currently resides on the server at Stevens Institute of Technology. I recommend 
that while reading the text, Aspen Plus be used simultaneously to execute and review 
each example. If Aspen Plus is not available, the .txt solutions may be reviewed. 

The book is designed to be used by undergraduates, graduate students, and practic¬ 
ing chemical engineers. The first section of the book explains the basic structure of the 
software and leads the student through a hands-on introduction to the various features 
of the software designed to facilitate the setup of simple problems. Features such as 
the material-balance-only option, access to Aspen Plus documentation through Help, 
the Next button, menu navigation, and the report function are introduced. The remain¬ 
der of the book is organized in a series of sections that focus on particular types of 
operations: for example, a two-phase flash. Each chapter is accompanied by the equiv¬ 
alent of lecture material that describes the equations being solved, various limitations, 
potential sources of error, and a set of workshops containing exercises that the students 
should solve to gain experience with the particular subject. Some of the exercises are 
designed to produce errors that students need to analyze in order to complete their 
experience. Much of this part of the book is suitable for undergraduates, although 
some will be limited by courses in their curriculum that have not yet been taken (e.g., 
exposure to the thermodynamics of phase equilibrium). Undergraduates should limit 
their exposure to Chapters Six and Fourteen. Chapter Six deals with phase equilibrium 
and provides exposure to the most popular thermodynamic equations as well as mate¬ 
rial on parameter fitting. Chapter Fourteen addresses advanced problems in distillation. 
Graduate students and practicing engineers who undertake these sections should have 
had exposure to undergraduate equilibrium stage operations and, preferably, a graduate 
course in thermodynamics. 

This book is not intended to be a self-study guide for all the features of Aspen Plus. 
For example, material on some reactor blocks, batch blocks, and the thermodynamics 
of electrolytes is not covered. Many subjects not addressed can be found by selecting 
the Help button on the main Aspen Plus display. The philosophy of the book is based 
on the idea that once a chemical engineer becomes thoroughly facile in the use of the 
software and has a good understanding of the basic blocks, he or she should be able 
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to learn to use many of the unaddressed functions by applying the same philosophy 
as that of the text itself: namely, to study appropriate sections in chemical engineering 
textbooks that describe the subject matter and to familiarize oneself with the function’s 
implementation by reading Aspen Plus’s documentation and attempting a sample prob¬ 
lem. As an example, to understand the Aspen Plus electrolyte methodology, it would 
be useful to read the section on electrolyte equilibrium in Molecular Thermodynamics 
of Fluid-Phase Equilibria , 3rd ed., by J. M. Prausnitz, R. N. Lichtenthaler, and G. M. 
de Azevedo (Prentice Hall, 1999) and in Aspen Plus’s Help, and follow up with the 
section entitled “Generating Electrolyte Components.” 

I have made an effort to provide the describing equations of most of the models 
(blocks) referred to here, and if not possible because of the proprietary nature of 
the software, I have described the functionality. One should recognize that Aspen 
Plus is proprietary software and that the source code and implementation details are 
not available. Additionally, there are frequently several ways to solve the equations 
that describe the blocks, and there is no way to ascertain these details since Aspen 
Technology does not provide them. 

I wish to acknowledge the help provided by Aspen Technology’s academic support 
group, especially for the loan of an Aspen Engineering stand-alone license for use 
while I was out of the United States and unable to access the Stevens Institute of 
Technology server. 


Ralph Schefflan 


CHAPTER ONE 


INTRODUCTION TO ASPEN PLUS 


Aspen Plus is based on techniques for solving flowsheets that were employed by 
chemical engineers many years ago. Computer programs were just beginning to be 
used, were of the stand-alone variety, and were typically used for designing single 
units. The solution of even the simplest flowsheet without recycle required an engineer 
to design each unit one at a time and, manually, introduce the solution values of a 
previously designed unit into the input of the next unit in the flowsheet. When it 
became necessary to deal with a recycle, the calculations began with a guess of the 
recycle values, and calculations ended when the values produced by the last unit in 
the loop matched the guesses. This involved much repetitive work and convergence 
was not guaranteed. This procedure evolved through the construction of rating models 
of units, as opposed to design models, which could be tied together by software in a 
way that emulated the procedure above and employed robust mathematical methods 
to converge the recycle elements of the process. This type of system is termed a 
sequential modular simulator. An excellent example of such software was Monsanto 
Corporation’s Flowtran (1974), which eventually became the kernel upon which Aspen 
Plus was built. 

Subsequently, Aspen Plus, although still basically a sequential modular simulator, 
has grown considerably and has many advanced functionalities, such as links to a 
variety of specialized software, such as detailed heat exchanger design, dynamic sim¬ 
ulation, batch process modeling, and many additional functions. It also has a facility 
for using an equation-based approach in some of its models, which permits convenient 
use of design specifications in process modeling. 

The Aspen Engineering Suite, which incorporates Aspen Plus, can be installed 
in a variety of ways using network servers or on a stand-alone personal computer. 
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2 INTRODUCTION TO ASPEN PLUS 


Installation is the responsibility of either the user, with tools provided by Aspen Tech¬ 
nology, or the information technology department, which services the user. This is done 
only once and modified, typically annually, with future releases. Whether the user’s 
installation is by network downloads or by CDs, it is necessary that the user select 
desired modules, at a minimum Aspen Plus and its required add-ons and associated 
documentation. No other modules are necessary. 

1.1 STARTING ASPEN PLUS 

When the Aspen Plus User Interface icon on the desktop is clicked or, alternatively, 
the sequence All Programs/Aspen Tech/Process Modeling V7.0/Aspen Plus/Aspen Plus 
User Interface is selected from the Start menu, the Aspen Plus Startup display shown 
in Figure 1.1 appears. 

A selection from a list of existing applications may also be made, or other stored 
files may by selected by choosing the appropriate radio button. If the Template button 
is selected, a list of application-oriented possibilities is shown (Figure 1.2). Although 
the flowsheet selection appears in the Run Type box in the lower right-hand corner of 
the display, other selections are available from the associated drop-down list. 



Figure 1.1 Aspen Plus startup. 
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Figure 1.2 Preconfigured selections. 


If the option “blank simulation” is chosen, an application can be custom configured. 
After selecting the desired option and clicking OK, the Connect to Engine screen 
appears, and upon selecting OK, a blank workplace that facilitates the graphic users 
interface appears if “flowsheet” was chosen. If the option chosen is anything else, the 
first required input form appears. 


1.2 GRAPHIC USERS INTERFACE 

The graphic users interface (GUI) is the means by which a flowsheet is defined. The 
process consists of placing blocks and streams on the workplace and connecting them. 
Aspen Plus assigns generic names such as B1 to the blocks. The user may change these 
names by right-clicking on the element of interest and using the menu that appears. 
Blocks are selected by choosing a category tab from the model library—for example, 
Mixers/Splitters—and clicking on the icon that represents the block desired. After 
clicking, a movable + sign appears on the open area of the display. After positioning 
it on the screen, a left click will place the block. The + sign remains and can be 
moved to insert another instance of the same block. This function ceases when the 
arrow button at the lower left is selected. In identical fashion, streams can be placed 
on the flowsheet. Material, heat, and work streams may be selected. When a stream 
input is selected and the cursor is moved onto the workplace, the ports to which the 
streams may be connected are shown. The connection is made by moving the active 































4 INTRODUCTION TO ASPEN PLUS 



Figure 1.3 Connecting streams. 


cursor over an open port and clicking. An example of connecting streams to ports is 
shown in Figure 1.3. 

All icons, block names, and stream names can be selected and moved using standard 
Windows techniques. Similarly, streams can be moved, rerouted, disconnected, and 
reconnected. Selecting and right-clicking on any of the objects displays a menu that 
provides many useful functions for manipulating the graphics. These include changing 
icons, rotating objects, renaming, deleting, and aligning the graphics. 

1.3 NEXT BUTTON 

Aspen Plus provides the user with a mechanism for filling out forms in an orderly 
fashion. At any point after the flowsheet has been fully defined with the GUI, the 
user may select the Next button, which appears as the symbol N on all forms. 
The Next button moves to the next form required. On occasion, after using the Next 
button, Aspen Plus will prompt the user to select from a choice of actions to be taken. 
The Next button provides only the minimum required input. As an example, when an 
activity coefficient equation to be used in the simulation is chosen, Aspen Plus will use 
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a default data source, such as a vapor-liquid equilibrium (VLE) source; however, if 
the simulation involves liquid-liquid equilibrium (LLE), it is the users’ responsibility 
to select the appropriate data source. Aspen Plus will not open the appropriate displays 
by using the Next button. 

1.4 SETUP SPECIFICATIONS DISPLAY 

After the flowsheet has been defined with the GUI, pushing the Next button brings 
up the Setup Specifications—Data Browser display. The Data Browser panel of the 
display should show a list of all possible categories that can be chosen for selecting 
various options and for entering appropriate data. If the setup information shown in 
the left panel of the display does not appear, choosing the upper menu selection, Data, 
and entry Setup will present the Data Browser panel, shown in Figure 1.4. If the list 
of process model icons shown at the bottom of Figure 1.4 does not appear, clicking 
the View menu Model Library will display them. 

When starting a problem, the Setup Specification display provides a drop-down list 
associated with the entry box Run type, which shows the six primary functions that 
Aspen Plus is capable of: 
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Figure 1.4 Setup specifications for data browser. 
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1. Data regression: fitting data to models 

2. Flowsheet: process simulation 

3. Property display: showing properties of a component in Aspen Plus’s database 

4. Property analysis: estimating physical and thermodynamic properties 

5. Assay data analysis 

6. Property plus 


The last two functions are not considered here. The user selects the required function 
to initiate data input for a specific requirement. 

Note that the Data Browser panel shows a blue check next to various items. This 
indicates that the default values are sufficient to proceed with data input; however, this 
is the minimum data required and the values may be modified to meet the requirements 
for any problem. The red elements on the browser panel list indicate that user input 
is required. A red element may be deleted by right-clicking on it and selecting Delete 
from the menu that appears, or if Aspen Plus permits, and sometimes by using the 
Windows delete key. 

There are drop-down lists for all the categories in the data-entry boxes under the 
Global tab, and all may be changed. For example, the value Mole has been selected for 
the option Flow basis. Additionally, under units of measurement, the input and output 
units have engineering (ENG) as default values, while the associated drop-down list 
offers Metric and SI units as other possibilities. 

Aspen Plus provides a Help button, on the topmost menu, for accessing information 
by subject. Additionally, help with any entry on any display is available by moving 
the cursor to the entry and pushing FI or selecting the \? button on the tool bar. 

1.5 SIMULATION OPTIONS 

If the Simulation Options category is selected, the Simulation Options display shown 
in Figure 1.5 appears. All of the default values that appear under the various tabs need 
not be changed except under the tab Calculations, where the option not to use energy 
balances in the calculations is available. This is an important option for preliminary 
calculations. If simulations do not involve solids or electrolytes, the appropriate options 
may be unchecked. 

1.6 UNITS 

Aspen Plus provides a user with a choice of units: engineering, metric, and the inter¬ 
national system of units, SI. An important option is the ability to select mixed units; 
for example, the choice of engineering units with mmHg and degrees C as temperature 
and pressure is not uncommon in some pharmaceutical applications. To accomplish 
this, Units-Sets is selected from the Data Browser menu, which produces the Object 
Manager shown in Figure 1.6. Then, selecting New produces an Identifier for the new 
unit set, in this case US-1, and a choice as to whether or not to assign this as a global 
data set. Then Figure 1.7 appears. 


UNITS 



Figure 1.5 Simulation options. 
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Figure 1.6 Units-sets and the object manager. 


































































































8 INTRODUCTION TO ASPEN PLUS 


inputexample.bkp - Aspen Plus V7.0 - a 


File Edit View Data Tools Run Plot Library 

DMlal a|a| Nb| n?| H 


Costing Help 


Jd*J 


Ej M| oj [QOlOll @|| |l| Costing: Uninil 


»|alBl B-|lHg| 


Setup Units-Sets US-1 - Data Browse 


JnjxJ - 


|© US-1 


~ZI <H~H <<||lnput d>>| 


N»| * k|x| 


B-SI Units-Sets 3 

; © ENG 
© MET 
I © SI 

© US-1 

| © Custom Units 
1 © Report Options 
3~Gfl Components 

© Specifications 
Q Assay/Blend 
O Light-End Propel 
El O Petro Characters 
O Pseudocompone 
O Attr-Comps 
Cl Henry Comps 
O Moisture Comps 
O UNIFAC Groups 
Cl Comp-Groups B 


'/Standard| Heat J Transport | Concentration | Size J Currency | Miscellaneous | 


T3 


Flow related 


Mass flow: 

flb/hr 

Mole flow: 

| Ibmol/hr ▼] 

Volume flow: 

| cuft/hr 

Flow: 

Ib/hr 

Flux: 

cuft/sqft-se ▼ | 

Mass flux: 

| Ib/sqft-hr ▼ I 


T emperature related 
T emperature: 

Delta T: 

Inverse temperature: 


F~jO 

I 1 R jlI 



Pressure: 

| psia 

”3] 

Delta P: 

psi 

B 

Delta P / Height: 

in-water/ft ^ | 

Head: 

| ft-lbf/lb 

T3 

Inverse pressure: 

11 /psi 



* 

-CH 

Material 

STREAMS 


Mixers/Splitters | Separators | Heat Exchangers | Columns | Reactors | Pressure Changers | Manipulators | Solids | User Models | 


Mixer FSplit S Split 


B^ 

Cono r < I ► 


For Help, press FI 


—rr 


| |c:\...en Plus Textbook\proposal | |NUM | |Flowsheet Not Complete 


Figure 1.7 New units. 


The Copy from entry permits selection of which unit set to use as a base. Note the 
selection of degrees C from the temperature drop-down list. Each unit’s entry has an 
appropriate drop-down list, which will enable a user to customize units. 

Custom units for a specific variable—for example, special composition units—may 
be defined by selecting the Custom Units entry under the Setup list. 

Aspen Plus input forms, displays, and reports are generated using the units selected. 
Input displays typically have a drop-down list adjacent to input boxes which permits 
the user to select the units of the input data required. This does not, however, affect 
the units used for output. 

Selection of the Report Options entry under the Setup list permits customization of 
the Aspen Plus results displays and the. txt report, which can optionally be generated 
when results are available. Details are given in Section 1.11. 


1.7 COMPONENTS 

Selecting the Components item on the Data Browser panel produces Figure 1.8. Com¬ 
ponents may be entered by either component name or chemical formula. Aspen Plus 
also requires that a nickname, Component ID, which is used in all stream reports, be 
provided for all components. 
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Figure 1.8 Component selection. 


An entry in the box Component ID is a user-provided component short name 
employed by Aspen Plus for report purposes and in some cases, such as water, is 
recognized as the component water. An entry is always required. Alternatively, the 
user may enter a proper component name or component formula. If neither is rec¬ 
ognized as an entry in the database, the user may select the Find button and Aspen 
Plus will display a set of names or formulas that incorporate the entry. For example, 
entering the formula C 7 H 8 gives the results shown in Figure 1.9. Upon selection of 
the component of interest, pressing the Add selected compounds button enters the 
component into the data associated with the current Aspen Plus run. 

If a component does not exist in the Aspen Plus database, choosing User Defined 
in Figure 1.8 produces the User Defined Component Wizard shown in Figure 1.10. 

Note that several properties are required. After entering the Component ID and 
formula, the Basic Data and Molecular Structure displays shown in Figure 1.11 appear. 
Known experimental data and structural information are entered. As the data-entry 
process proceeds, an option for Aspen Plus to estimate any missing data appears. 
Chapter Two concerns methods for estimating data which may be employed selectively 
for supplying the missing data above rather than permiting Aspen Plus to provide 
missing data using default methods. 
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Figure 1.9 Component search. 


1.8 PROPERTIES 

The first display when selecting properties for a simulation is given in Figure 1.12, 
where the choice of a global property method that applies to an entire flowsheet may 
be selected. If, however, the simulation involves a situation where more than prop¬ 
erty method is required, for example, a process involving distillation (vapor-liquid 
equilibrium) and extraction (liquid-liquid equilibrium), each block in the flowsheet is 
identified as part of a unique flowsheet section when the flowsheet is created, and is 
identified with a particular property method by choosing the Flowsheet Sections tab. 
A display similar to Figure 1.12 appears; the properties method for each flowsheet 
section can be identified and a global property method is not used. The Tools menu 
contains the Property Method Selection Assistant, which can be used as a guide for 
selection of a property method for specific applications; alternatively, the small box 
to the left of the Uniquac selection in this example, although not explicitly labeled, 
performs the same function. 

For this example process, sections S-l and S-2 are assigned the property methods 
Uniquac and Uniq-2. Aspen Plus sets up the required Uniquac binary parameters for 
each selection. Care must be taken that the appropriate database is accessed for each set 
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Figure 1.10 User defined component. 



Figure 1.11 User defined component data and molecular structure. 
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Figure 1.12 Global property method. 


of parameters. As an example, if the Data Browser topic Properties/Parameters/Binary 
Parameters is selected and the listing Uniq-1 is clicked, Figure 1.13 appears. This shows 
the values of Aspen Plus-supplied parameters for each binary pair and the source of 
the data used in the regression. In Figure 1.13, note that the source of equilibrium 
data for each binary pair is identified. A selection of data sources is available by 
choosing the tab Databanks, which produces Figure 1.14. In the current example, 
section S-l contains the extractor, and it is therefore important that the data source be 
liquid-liquid equilibrium data. Section S-2 involves a distillation column; hence the 
binary parameters source is vapor-liquid equilibrium data. 

1.9 STREAMS 

Pushing the Next button moves the input sequence to stream data. All feed streams 
are defined using a display such as that given in Figure 1.15. The data entry is very 
straightforward and provides a user with the possibility of changing the units of both 
the material flow and the state variables. For a stream of a single component or several 
components, the number of degrees of freedom calculated by the Gibbs phase rule 
applies: 


F = C-P +2 
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Figure 1.13 UNIQ-1 parameters from Aspen Plus databank. 



Figure 1.14 Available databanks. 
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where F is the degrees of freedom, C the number of components, and P the num¬ 
ber of phases. For a single-phase system, F = 2; therefore, two specifications are 
required. For an n -component system, F = n + 1. Since n — 1 mole fractions need 
to be specified and the last calculated by the sum of the mole fractions equal to 1, 
only two additional specifications are required. In both cases these are usually, but not 
necessarily, temperature and pressure. 

When two phases in equilibrium are involved, only one degree of freedom is avail¬ 
able. For example, for a one-component stream such as a saturated liquid, specification 
of the temperature fixes the (vapor) pressure. But in such circumstances it is neces¬ 
sary to state the fraction of the mixed stream that is vapor (or liquid). For a saturated 
liquid the V/F specification would be a very small number, such as 0.00001. For a 
multicomponent stream the situation is identical and it is necessary to specify either 
temperature or pressure and the vapor or liquid fraction. 

If the process contains tear (recycle) streams, they will not be treated as required 
input, and pushing the Next button may not suffice. Typically, Aspen Plus will assign 
zeros as starting values to the variables that are to be converged, but if a user wishes 
to provide starting values, the stream name under the streams list in the Data Browser 
can be clicked and a display analogous to Figure 1.15 will appear. 

1.10 BLOCKS 

When all stream input has been completed, pushing the Next button will result in 
the appearance of the first input form for a block that requires data. Details of block 
input are addressed in other chapters. After the data input forms for the first block are 
completed, pushing the Next button will produce the forms for the next block in the 
process until all the block data have been entered. 

The browser list in Figure 1.15 shows some additional topics that may be appropriate 
for a user’s simulation. These will not appear when the Next button is pushed; however, 
they may be selected by clicking on the subject. For example, clicking the Convergence 
entry permits the selection of the convergence method and parametric default values. 
These subjects are addressed in other chapters. When data entry for all blocks and 
supplementary data entry is completed, selecting the Next button produces a dialogue 
to enable execution. All input can be reviewed, prior or after execution, by selecting 
the or keys above the Flash Options tab near the top of Figure 1.15 and selecting 
input from the drop-down list between the double arrows. 

1.11 VIEWING RESULTS 

In preparation for viewing the results, prior to execution of a simulation, a user may 
customize reports by selecting Report Options, as shown in Figure 1.16. Customization 
is available under each of the tabs. The stream tab illustrates the variety of options 
available for selecting the stream information to be displayed. 

At the conclusion of a simulation, the run control panel is usually displayed. If 
not, it can be viewed from the drop-down list under the main menu item View. This 
contains a brief summary of the execution of each block and a list of error messages. It 
is important that users check and correct any errors that might have occurred. When a 
process contains tear streams, errors that occur during the iterations will be presented. 
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Figure 1.15 Stream specifications. 
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Figure 1.16 Report options. 
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Figure 1.17 Run control panel. 


It is not uncommon that errors occur before convergence, but when the process has 
converged, there should be none. Figure 1.17 displays an example of the run control 
panel. The results of the simulation can be viewed by selecting the Results button, 
which contains a check symbol overlaid on a file symbol on the fourth row from the 
top of the display. 

Selection of the Results button produces a display similar to Figure 1.15, with the 
central drop-down list displaying results. Selection of the or keys permits paging 
through all the results. If All is selected from the drop-down list, one may view both 
input and results in order. 

When a simulation has been completed, selection of Report from the main menu 
item View will produce a. txt report of the complete simulation, in a Notepad window 
that can be copied and pasted. 

The file Chapter One Examples/inputexample.bkp was used to produce all of the 
figures in this chapter as well as the inputexample.txt file. 

1.12 OBJECT MANAGER 

In certain data input situations there are specialized setups that involve the specification 
of material that requires a special set of dialogues to define what is required. Examples 
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Figure 1.18 Object manager. 


are: definition of a set of properties to be displayed, regression of a data set to a 
thermodynamic model, and specification of the value of a block output variable. This 
type of input is managed through use of the Object Manager. These situations typically 
begin with the identification of a function, such as a regression, followed by a series 
of choices and specifications. An example of the use of the Object Manager for a 
Calculator block, which is used for in-line calculations, is given in Figure 1.18. The 
user creates an ID such as C-l, which is followed by all the input forms required, 
including association of calculation variables to flowsheet variables and definition, in 
Fortran, of the calculation details. 


1.13 PLOTTING RESULTS 

Aspen Plus has a plotting facility accessed through the Plot menu. The facility either 
provides the selection of independent and dependent variables manually or by means 
of a Plot Wizard, which permits the user to select a preconfigured plot. As an example, 
after obtaining a converged solution to a distillation column, selecting the Plot Wiz¬ 
ard, and pushing the Next button, a collection of distillation-oriented plots, shown in 
Figure 1.19 is available. Selecting the Comp plot produces Figure 1.20, which shows 
the change in vapor composition through the column. For situations where preconfig¬ 
ured plots are not available but a display of tabular output is available, such as the 
results of a sensitivity analysis (details in Chapter Five), one may select a column of 
data and assign it to either the v or y coordinate. It is also possible to assign two vari¬ 
ables to the y coordinate by holding down the control key when selecting the columns 
of data that are to be employed. If the scale or title of an axis or title of a plot is 
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Figure 1.19 Plot wizard. 
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Figure 1.20 Sample plot. 
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not suitable, one may select it by clicking, and a display that offers editing options 
appears. 


REFERENCES 

Aspen Plus version 7.0 documentation. 

Monsanto Corporation, Flowtran Simulation: An Introduction, 1974. 


CHAPTER TWO 


PROPERTIES 


Aspen Plus offers two possibilities for accessing properties: Property Analysis and 
Property Estimation. Each is invoked from the drop-down, Run-type menu (Figure 2.1). 
The Property Analysis function can display pure component values such as the criti¬ 
cal compressibility factor, temperature-dependent properties such as the ideal gas heat 
capacity, and mixture properties from a variety of data banks and mixture properties 
from multicomponent functions such as equations of state. The Property Estimation 
capability can be used to estimate virtually the same values as are stored in the Aspen 
Plus database for user-defined components. References in this chapter to pure compo¬ 
nent properties, equations of state, activity coefficient equations, and property equations 
can all be found in the Aspen Plus Physical Property System documentation. 

2.1 PURE COMPONENT DATA BANKS 

All of the many data banks available in the Aspen Physical Property System can be 
identified by clicking the Help button at the top of Figure 2.1. The primary database is 
Pure22. To see the details of its content, searching Help for Pure22 will present various 
alternatives, one being to display the Pure22 Databank. The list below, taken from 
Aspen Plus’s documentation, describes the property categories for which parameters 
are stored. 

• Universal constants, such as critical temperature and critical pressure 

• Temperature and property of transition, such as boiling point and triple point 

• Reference-state properties, such as enthalpy and Gibbs free energy of formation 

• Coefficients for temperature-dependent thermodynamic properties, such as 
liquid-vapor pressure 
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Figure 2.1 Property analysis drop-down menu. 


• Coefficients for temperature-dependent transport properties, such as liquid 
viscosity 

• Safety properties, such as flash point and flammability limits 

• Functional group information for all Unifac models 

• Parameters for Soave-Redlich-Kwong and Peng-Robinson equations of state 

• Petroleum-related properties, such as API gravity, octane numbers, aromatic con¬ 
tent, hydrogen content, and sulfur content 

• Other model-specific parameters, such as the Rackett and Uniquac parameters 


The list of components stored can be found by clicking the link Pure Compo¬ 
nent Databanks and at that link, clicking the Pure Component Databank Parameters, 
which provides the parameter names, descriptions, units, and availability in other 
Aspen Plus data banks. An example of such data is shown in Figure 2.2. This dis¬ 
play is invoked by selecting the Run type Property Analysis, entering components (in 
this example, methanol and water), and selecting Tools, Retrieve Parameter Results, 
and OK. The display will be found under Results, Pure Component. Pure compo¬ 
nent temperature-dependent results are also available from the T-Dependent tab; for 
example, the selection of parameter Plxant-1 refers to the coefficients of the extended 
Antoine equation, as shown in Figure 2.3. The extended Antoine equation is 


C 2 


In p* = C x + 


t + c 3 


+ C 4 T + C 5 In T + C 6 T Cl C 8 < T < Cg 


( 2 . 1 ) 
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Figure 2.2 Pure component scalar properties. 
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Figure 2.3 Pure component temperature-dependent properties. 
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Figure 2.4 Temperature-dependent binary interaction parameters. 


where the element numbers in Figure 2.3 correspond to the parameters of equation (2.1). 

Selection of the Results Binary Interactions with the tab T-Dependent and parameter 
Wilson-1 displays the binary interaction parameters of the methanol-water system as 
shown in Figure 2.4. The elements refer to the Aspen Plus implementation of the 
Wilson equation for a binary pair: 

l n Y ij — a ij + ■— + c ij In r + dyT + (2.2) 

where i and j refer to the components and the element numbers correspond to the 
parameters of equation (2.2); for example a t j is element 1, and elements 9 and 10 refer 
to the applicable temperature range. Similar results appear when selecting a different 
activity coefficient equation. 

2.2 PROPERTY ANALYSIS 

Property Analysis enables the display of temperature-dependent properties at user- 
selected temperatures as well as phase equilibrium data for both two and three phases. 
Display of user-selected properties requires an analysis to be defined, including the vari¬ 
able range over which the properties are to be presented. Aspen Plus presents forms 
that need to be completed, similar to flowsheet run forms, followed by a properties 
specification input. Clicking the Next button initiates a display in which the property 
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Figure 2.5 Example of properties analysis specification PT-1. 


analysis to be generated is specified. For example, a vapor-liquid equilibrium prop¬ 
erties specification, labeled PT-1, is shown in Figure 2.5, where points along a flash 
curve are chosen in anticipation of preparing a T-x-y diagram. Aspen Plus requires 
the specification of composition, even though it is not needed for this example. Note 
the appearance of the adjusted variable range display shown in Figure 2.6. Prior to 
completing the analysis input, Aspen Plus prompts the user to define which proper¬ 
ties are required. These are defined in a Prop-Set, with an example labeled PS-1, and 
displayed in Figure 2.7, where the properties required are the bubble and dew points 
at 14.7 psi as a function of composition. It is possible for several property analyses 
and property sets to be executed with one Aspen Plus run, and a display is provided 
to associate a Prop-Set with an analysis. An example of such an association is shown 
in Figure 2.8. The results are the tabular display shown in Figure 2.9, which can be 
plotted with Aspen Plus’s Plot function, as shown in Figure 2.10. 

Another example of a Prop-Set, PT-2, which deals with an analysis of temperature- 
dependent properties, is shown in Figure 2.11. The tab Qualifiers is used to define the 
components, the phases, and other specifications for which properties are desired with 
the associated Prop-Set PS-2. Tabular results are shown in Figure 2.12. A property 
analysis example may be found at Chapter Two Examples/property analysis example. 

2.3 PROPERTY ESTIMATION 

There are many methods for estimating pure component properties such as the critical 
pressure. When using Aspen Plus to estimate such values, a user is offered a choice 
of methods, most of which are described by Poling et al. (2000), and in some cases, 
in an earlier edition of the same book. As an example, Joback’s method for estimating 
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Figure 2.6 Range definition for adjusted variable. 



Figure 2.7 Example of Prop-Set PS-1. 
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Figure 2.8 Association of PS-1 with PT-1. 



Figure 2.9 Results of PT-1. 
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Figure 2.10 Plot of PT-1 results. 



Figure 2.11 Variable selection for Prop-Set PS-2. 
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Figure 2.12 Tabular results for analysis PT-2. 


the critical temperature, is given by 


T C (K) = T b ■ 

0.584 + 0.965 

Y J N k(tck) 

_ 

r i 2 ' 

J2 N k(tck) 



_ k 


_ k 


where tck are structural group contribution values and Nt are the number of appear¬ 
ances of that group in the molecule. Note that the critical temperature is directly 
proportional to the boiling point modified by group contribution functions. 

Equations of state provide the means for estimating many other properties, such 
as the effect of pressure on enthalpy at constant pressure. Since the equation-of-state 
parameters are composed primarily of critical properties and in many cases the acentric 
factor (derived from vapor pressure), experimental values for the boiling point and 
vapor pressure data are very important. 

The property constant estimations system within Aspen Plus, sometimes referred to 
as PCES, is initiated by the same methods as in Chapter One, except that the Run type 
Property Estimation is entered. This function is used when a user-defined component 
is required. The component should not be entered into the component list, but rather, 
the new component wizard, shown in Figure 2.13, should be employed. Aspen Plus 
prompts for any available basic data and molecular structure, as shown in Figure 2.14 
followed by Figure 2.15 for additional data. It is imperative that the radio button to 
estimate properties by the Aspen Plus system be selected. Methyl vinyl ketone (MVK) 
is an example of a user-defined component. The molecule’s structure is 


CH 3 -CO-CH=CH 2 
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User Defined Component Wizard I? | X| 


Welcome to the User Defined Component Wizard, the quickest way to 
enter properties for User Defined Component.This wizard will lead you 
through the steps to enter the required physical properties for the 
User Defined Component based on its type. 


Component ID: |MVK 
Formula: |C4H6Q 


Type: 


Conventional 




The required properties for Conventional components include 

Molecular weight, normal boiling point, 
molecular structure, vapor pressure and 
ideal gas heat capacity. 


Click Newt> to continue or Finish to accept the component and exit. 


Cancel (Back Nent> Finish 


Figure 2.13 User defined component wizard. 


The molecule’s structure, defined in Figure 2.14, is entered by assigning a number to 
each nonhydrogen element and defining the bond type between them. In this case the 
carbons are numbered from the left 1, 2, 4, and 5, and the oxygen is numbered 2. 

For a multicomponent system, components stored in an Aspen Plus database or 
additional user-defined components may be added. Pushing the Next key presents 
Figure 2.16, which enables a choice of estimation options. If the choice to, estimate 
all missing parameters is made, it is possible that only the parachor will be estimated. 
It is preferable to choose the radio button and estimate only the parameters selected, 
in which case check boxes for selection of parameter types become available and 
each property desired can be selected via the parameter and/or property lists in the 
tabbed categories. Figure 2.17 shows an example in which all estimation methods 
have been selected. This may generate a warning message for some methods, due to 
their limitations on execution. The temperature-dependent properties input are similar to 
pure component properties except that an independent variable range must be specified. 

The tab Binary is used to estimate ciy pairs for activity coefficient equation pairs from 
estimates of infinite dilution activity coefficients by one of the Unifac methods available 
in Aspen Plus. It is also possible to estimate the ky parameter for the SRK equation of 
state. The component list must contain additional entries which specify the components 
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Conventional Component Basic.. ■ Molecular Structure-MVK 


al Functional Group | Formula | Structure | 


Component ID: MVK 
Enter Molecular structure 


Import Mol file 


r 

C Draw structure 

(* | Input molecule by its connectivity | Connectivity form was 


Enter available physical properties 
Molecular weight: 

Normal boiling point: 

Specific gravity at 60 deg F: 
Standard enthalpy of formation: 
Standard Gibbs energy of formation: 


;ra 

| Btu/lbmol 


| Btu/lbmol 


Click Next> to continue or Finish to save the above properties and exit. 


YGei 



Atoml 

Number Type 

Atom2 

Number Type 

Bond type 



C 

2 

C 

Single bond 


2 

C 

3 

0 

Double bond 


2 

c 

4 

C 

Single bond 


4 

c 

5 

c 

Double bond 

* 







Atom number ■ atom type correspondence 










Atom number 


2 

3 

4 

5 


Atom type 

C 

C 

0 

C 

C 


:\...pplication Data\AspenTech NUM 


|Required Input Incomplete 


| <Back | Next> | Finish | 


Figure 2.14 User defined component data. 


Conventional Component Addit... 


Click buttons 1 to 5 to enter additional properties or data, 
m Liquid density 
2 | Vapor pressure data 

Extended Antoine vapor pressure coefficients 

4 | Ideal gas heat capacity data 

5 | Ideal gas heat capacity polynomial coefficients 


C Evaluate using TDE 


Evaluate using TDE 


f* Estimate using Aspen Plus property estimate system 


Cancel 


<Back 


Next> 


Finish 


Figure 2.15 Additional user defined component wizard. 
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Figure 2.16 Estimation setup. 



Figure 2.17 Pure component parameter selection. 

with which the user-defined component will be interacting. Parameter estimation can 
also be done using experimental data. An example is shown in Figure 2.18. 

Occasionally, it may be necessary to supplement the structural data supplied in 
Figure 2.14 with group contribution data for one of the methods available. An example 
of Unifac data for the MVK molecule is given in Figure 2.19. The identification of 
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Figure 2.18 Binary parameter estimation setup. 



Figure 2.19 Unifac group selection. 
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the group numbers is provided at the bottom of the display, but they cannot be seen 
unless the display is full size. The setup and results for this example may be found at 
Examples/mvkwaterestimation. 

2.4 WORKSHOPS 
Displaying Properties 

Open an Aspen Plus run selecting the blank simulation template with the Run type 
Property Analysis. Use the Next key to proceed with the data input. The components 
to be used are mesityl oxide and toluene. Select Unifac and Ideal to represent the vapor 
and liquid phases, respectively. 

Workshop 2.1 Select the Generic option for your property type and enter 50 mol/hr 
for each component when required. Generate a T-x-y diagram at 14.7 psi. Use a 
Property Set to choose which properties are required. Plot the T -x -y diagram. Display 
all pure component parameters by using the Tools menu. Save the values of the boiling 
point, critical temperature, and critical pressure. 

Estimating Properties 

Open an Aspen Plus run by selecting the blank simulation template with the Run type 
Property Analysis. Use the Next key to proceed with the data input. The components 
to be used are mso-e (mesityl oxide) and toluene. Treat mso-e as a user-defined 
component even though it is in the data bank. Provide the molecular structure and 
molecular weight only. 

Workshop 2.2 Select Setup and estimation options, and estimate only the parameters 
selected and pure component scalar properties. Under the Pure Component tab, choose 
boiling point, critical temperature, and critical pressure for the component mso-e using 
all available methods. Interpret the error messages. 

Workshop 2.3 Remove the property estimation methods that caused the error mes¬ 
sage, but ignore the warning message. Reinitialize and rerun. 

Workshop 2.4 Repeat Workshop 2.3 but provide as input the true boiling point 
of mesityl oxide: 129.76°C. Be sure to clear all entries under the setup menu entry 
PCES. Clear the request to estimate the boiling point. Reinitialize and rerun. Compare 
the results of Workshops 2.3 and 2.4. 

Workshop Notes 

Workshop 2.1 Unifac is used to estimate the vapor-liquid equilibrium for this system 
to demonstrate the possibility of obtaining such data in the absence of experimental 
values. If possible, experimental verification is extremely important. 

Workshop 2.2 Although all methods for estimating the three properties appear avail¬ 
able, some of the methods have limitations. For example, Ogata’s method for estimating 
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boiling points is limited to one radical per molecule. Some of these difficulties inhibit 
execution of the Property Estimation run. 

Workshop 2.3 Upon removal of those methods with problems from the property 
specifications, the execution is satisfactory despite the warning message concerning 
the boiling point and vapor pressure. Note that the two methods for estimating the 
boiling point are about 10 K apart. Notice also by comparison to the database values 
presented in Workshop 2.1 that the estimated boiling points by both methods are off 
by 5 to 15°F. 

Workshop 2.4 Comparing the critical temperatures predicted by Workshops 2.3 and 
2.4, where the experimental boiling point is included in the data, the Workshop 2.3 
data average about 589 K, and the Workshop 2.4 data, 597 K. The critical pressure 
does not depend on the boiling point. 


REFERENCES 

Aspen Plus version 7, Physical Properties System documentation. 

Poling, B. E., Prausnitz J. M., and O’Connell, J. P., The Properties of Gases and Liquids, 5th ed., McGraw- 
Hill, New York, 2000, Chaps. 2 and 3. 


CHAPTER THREE 


THE SIMPLE BLOCKS 


Aspen Plus provides several blocks that are suitable for, but not limited to, material- 
balance-only calculations: the stream splitter, the mixer, and the simple separator 
blocks. Additionally, the stoichiometric reactor, covered in Chapter Ten, is also very 
useful for material-balance-only calculations. These blocks are especially important 
since when they are combined to form a simple process model, the results are useful 
as an initial estimate for a rigorous model of a process. In some cases, such calculations 
are used as the basis for design. 

Aspen Plus also provides two blocks with general utility in various aspects of 
simulation. The stream duplicator block enables, for example, the execution of duplicate 
processes with different thermodynamics and parameters but the same feeds. The stream 
multiplier enables all elements of an input stream to be modified by a calculated 
multiplier subject to the specification of a block output, such as the composition of a 
component in an output stream. 

3.1 MIXER/SPLITTER BLOCKS 

The blocks Mixer, Fsplit, and S split are available from the Aspen Plus model library 
shown on the bottom of the Aspen Plus displays under the Mixers/Splitters tab, as 
shown in Figure 3.1. 

3.1.1 Mixer Block 

The Mixer block combines an arbitrary number of input streams and produces a single 
stream by a simple material balance. If the “no heat balance” option is chosen, no 
other input is required. An example is shown in Figure 3.2. If F t and hi represent the 
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f* 

-CM 


Mixers/Splitteis | Separators | Heat Exchangers | Columns | Reactors | Pressure Changers | Manipulators | Solids | User Models | Com_l 


Mixer FSolit SSolit 


Material 

STREAMS 1 Mixer_FSplit S Split 


Figure 3.1 Mixer s/splitters. 



Figure 3.2 Mixer block. 


flow rate and enthalpy/mole of stream i , the material balance and energy balance are 

F 4 = F l + F 1 + F 3 (3.1) 

h\F\ = h\F\ + /?2^2 “t - (3.2) 

When energy balances are selected, it is necessary to specify the outlet stream’s 
pressure or the pressure drop across the block. If no pressure is specified, the block 
uses the lowest pressure of the feed streams. Additionally, the allowable phases of the 
output stream must also be specified. Since the output stream is of known composition 
and pressure and its molar enthalpy is calculated using equation (3.2), its state is 
defined, and therefore an adiabatic flash yields its temperature. 

3.1.2 Fsplit Block 

The Fsplit block is designed to divide a single stream or a mixture of streams of 
known state (i.e., flow rate, composition temperature, and pressure) into an arbitrary 
number of product streams with identical states. Note that since the state is known, 
the enthalpy/mole is calculable. A choice of product stream split fraction (based on 
the sum of the flow rates of the feed streams), mass or molar flow rate, volumetric 
flow rate, or flow rate of a component may be made. If there are n product streams, 
specifications for any n — 1 streams must be made. Figure 3.3 shows an example of 
an Fsplit block. 

When the model is executed, all the feed streams are mixed and the combined flow 
rate, composition, and molar enthalpy is computed. It is important to note that any 
of a block’s product specifications can be converted to product stream split fractions. 
This provides the basis for the describing equations of the model. For example, if the 
molar flow rate of component i in stream j, p \, is specified, the split fraction a, can 
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Figure 3.3 Fsplit block. 


be calculated as 


a,- = 


Pj_ 

fi 


(3.3) 


where f is the combined feeds flow of component i . Componential material balances 
are as follows. Let F and h represent the total flow rate and the molar enthalpy of 
the combined feed streams, and let Pj and a, represent the flow rate and split fraction 
of the product streams. Since there are n product streams n — 1 split-fraction material 
balance equations, each described by the following equation, may be written for the 
product streams: 


Pi=oaF (3.4) 

with the flow rate of the last product stream calculated by the overall material 
balance 


n —1 


Pn=F-J2 P k 


k= 1 


(3.5) 


The combined feed molar enthalpy, h, is assigned to each product stream. The 
pressure of all the product streams is identical and is either specified on the flash 
tab or if not specified, the lowest input stream pressure is used and thus all product 
streams are at the same state. If the energy balance option is used for the simulation, 
the temperature of a product stream is calculated by an adiabatic flash. 

The Ssplit block available from the model library, except for the input format, has 
the same functionality as Fsplit. 

3.2 SIMPLE SEPARATOR BLOCKS 

The simple separator blocks Sep and Sep2 are available from the Aspen Plus model 
library shown on the bottom of the Aspen Plus displays under the Separators tab, as 
shown in Figure 3.4. 
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1 * 

-cm 


Mixers/Splitters Separalois | Heat Exchangers | Columns | Reactors | Pressure Changers | Manipulators | Solids | User Models | Com < H 


(0-0 -O- 

Flash2 Flash3 Decanter 


Material 

STREAMS 1 Flash2 Flash3 Decanter Sep Sep2 


Figure 3.4 Separators. 


3.2.1 Sep Block 

When the Sep model is executed, all the feed streams are mixed and the combined flow 
rate, composition, and molar enthalpy is computed. The block permits the assignment 
of the flow rates or componential split fraction of each component in each of the n — 1 
product streams, based on the combined feeds. The composition of the nth stream is 
calculated by the overall material balance. An example of a Sep block is shown in 
Figure 3.5. 

As part of the block’s input specifications, there is a tab, Outlet Flash, which enables 
a user to specify the state of each outlet stream. A flash of each stream establishes the 
molar enthalpies and facilitates calculation of the block energy balance. 

Since Sep is a typical sequential modular block, for an n -component, m -product 
stream configuration, with the temperature and pressure of the feeds and products 
specified, a degree-of-freedom analysis shows that for a material-balance-only appli¬ 
cation there are nm unknowns. If n(m — 1 ) products are associated with a specification 
such as the fraction of the sum of the feeds of a component, to exit in a particular 
stream, n unknowns remain to be calculated by a material balance; thus, the com¬ 
position is calculable for all streams. For an energy balance calculation there is one 
additional unknown per product stream: the molar enthalpy. Since the composition 
temperature and pressure of each stream is known at the conclusion of the material 
balance calculations, an adiabatic flash produces the molar enthalpy of all streams, and 
an overall energy balance can be performed. 

If f k represents the quantity of component i in feed stream k , the amount of 
component i entering the block, is given by 

/ 

Fi=Y,fi (3- 6 ) 

k=\ 



Figure 3.5 Sep block. 
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where / is the number of feed streams. If the split fraction of a particular component 
i , in stream j, a \, is given, the quantity of component i leaving the block in stream j , 
Sj , is given by 


Sj =ocjFi (3.7) 

Alternatively, the quantity of component i leaving the block in stream j, Sj, can be 
specified directly. These variables in the equations do not depend on each other, and 
thus the solution, where not uniquely specified as in equation (3.7), is direct. 

3.2.2 Sep2 Block 

The Sep2 block is designed to model a simple two-product separator such as a dis¬ 
tillation or extraction process. Like the Sep model, when a Sep2 block is executed, 
all the feed streams are mixed and the combined flow rate, composition, and molar 
enthalpy is computed. For n product streams, the block permits assignment of the 
flow rates or componential split fraction of each component in each product stream of 
n — 1 streams, based on the combined feeds. The Sep2 model, however, provides more 
flexibility in specification than that of the Sep block. The additional specs allowed are 
overall stream split, that is, the ratio of flow of a product stream to the sum of the flows 
of the feed streams, and the mole or mass fraction of a component in a product stream. 
The composition of the nth stream is calculated by the overall material balance. An 
example of a Sep2 block is shown in Figure 3.6. 

Unlike the Sep block, there is interdependency between some of the block specifica¬ 
tions. Equations (3.6) and (3.7) describe the specifications of the Sep block that apply 
to the Sep2 block, but the other Sep2 specifications involve interaction between the 



Figure 3.6 Sep2 block. 
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variables. The mole fraction specification is an example. If x\ represents the specified 
mole fraction of component i in product stream j , when / is the number of components, 
then 


si 


(3.8) 


EL, Si 

applies. An equation similar to (3.8) which describes the total flow specification is 


EL fj 

EL, F m 


(3.9) 


where P j is the overall split fraction of stream j to the sum of flows of all p feeds, F m 
represents the sum of the flows of all components in feed stream m, and l represents 
the number of components. 

Note that the flows of the components may not be known a priori as in the Sep 
block; therefore, the solution of a set of simultaneous equations in which the unknowns 
are the individual component flows of product streams is required. 

Like the Sep block, there is a tab, Outlet Flash, which enables a user to specify the 
state of each outlet stream. A flash of each stream establishes the molar enthalpies and 
facilitates calculation of the block energy balance. 


3.3 SOME MANIPULATOR BLOCKS 

Manipulators may be found in the Aspen Plus library under the tab Manipulators 
(see Figure 3.7). The only manipulators considered here are the stream duplicator and 
stream multiplier blocks. 

3.3.1 Dupl Block 

An example illustrating the solution of the same problem with both the Sep and Sep2 
blocks is shown in Figure 3.8. Upon execution, streams 2 and 3 are identical to stream 1. 

A feed of 100 lbmol/hr of methanol, 200 lbmol/hr of water, and 300 lbmol/hr of 
ethanol is to be separated into an overhead product of 100 lbmol/hr with composition 
0.95 mole fraction methanol and 0.04 mole fraction ethanol. The column, operating at 
14.7 psi, is to be simulated twice, once with a Sep block and once with a Sep2 block. 
A Duplicator block is to be used to feed both columns. Figures 3.9, 3.10, and 3.11 
show the specifications for both separator blocks. 

Figure 3.12 shows that identical results using the file Chapter Three Examples/ 
DuplicatorExample are obtained; however, the use of Sep would not be straightforward 
if it were in the middle of a process where the feed was not known a priori. 


r*r 

-O* 

Material 

STREAMS 





Figure 3.7 Manipulators. 
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Figure 3.8 Duplicator block example. 


^Specifications| Feed Flash j Outlet Flash | Utility ] 


Outlet stream conditions 


Outlet stream: [T 
Substream: 


~TZI 

| MIXED 



Component ID 

Specification 

Basis 

Value 

Units 


METHANOL 

Flow 

Mole 

95 

Ibmol/hr 


ETHANOL 

Flow 

Mole 17 

4 

Ibmol/hr 


WATER 

Flow 

Mole 

1 

Ibmol/hr 


Figure 3.9 Sep specifications. 
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'/Specifications| Feed Flash | Outlet Flash | Utility | 


Substream: | MIXED ▼ | Outlet stream: |6 

Stream spec: plow _^jjl00 |Mole ▼||lbmol/hr ▼ | 


Component ID 

1st Spec 

2nd Spec 



METHANOL 

Split fraction 




Mole frac 

0.95 


ETHANOL 

Split fraction 




Mole frac 



WATER 

Split fraction 




Mole frac 





Figure 3.10 Sep2 specifications. 


^Specifications! Feed Flash j Outlet Flash | Utility | 


Substream: | MIXED Outlet stream: P-3 

Stream spec: J Split fraction ▼ \\~ 



Component ID 

1st Spec 

2nd Spec 


METHANOL 

Split fraction 




Mole frac 



ETHANOL 

Split fraction 


0.98 


Mole frac 



WATER 

Split fraction 




Mole frac 



Figure 3.11 Sep2 specifications. 


3.3.2 Mult Block 

The Multiplier block contains a user-specified parameter which multiplies all the flow 
rates in the block’s feed by a specified value or by a value calculated from another part 
of the process. An example of the use of the multiplier block is shown in Figure 3.13. 
In this example, stream 1, the feed to the Mult block, is specified, but stream 2, the 
feed to the rest of the process, is unknown. In this case a process quantity such as 
the flow rate of a component in a stream, for example in stream 8, will be targeted to 
achieve a user-specified value. A design specification, described in detail in Chapter 
Five, is used to manipulate the value of the multiplier parameter, operating on stream 
1, so that stream 2 attains a value such that the desired flow rate of the component in 
stream 8 is obtained. 
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Display: Format j FULL 
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Substream: MIXED 
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Mole Flow Ibmol/hr 











METHANOL 

100.0000 

100.0000 

100.0000 

95.00000 

5.000000 

95.00000 

5.000000 




ETHANOL 

200.0000 

200.0000 

200.0000 

4.000000 

196.0000 

4.000000 

196.0000 




WATER 

300.0000 

300.0000 

300.0000 

1.000000 

299.0000 

1.000000 

299.0000 




Mole Frac 











METHANOL 

.1666667 

.1666667 

.1666667 

.9500000 

.0100000 

.9500000 

.0100000 




ETHANOL 

.3333333 

.3333333 

.3333333 

.0400000 

.3920000 

.0400000 

.3920000 




WATER 

.5000000 

.5000000 

.5000000 

.0100000 

.5980000 

.0100000 

.5980000 




T otal Flow Ibmol/hr 

600.0000 

600.0000 

600.0000 

100.0000 

500.0000 

100.0000 

500.0000 

J 

i 

1 -lJ I 3 
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Figure 3.12 Duplicator example results. 



O 


Figure 3.13 Mult block. 


3.4 WORKSHOPS 

Workshop 3.1a Three streams are to be mixed. Their specifications are given in 
Table 3.1. Use the option “material balance only.” The property option to be used is 
Chao-Seader. Solve the flowsheet. The components are propane, -butane, n -pentane 
and n -hexane. Stream specifications are given in Table 3.1. 

Workshop 3.1b Repeat Workshop 3.1a using material and energy balances. What 
do you observe? Is your solution correct? 

Workshop 3.2a Three streams are to be fed to a Fsplit block used to emulate 
a distillation column with a sidestream. The process specifications are given below. 
Use the option “material balance and energy balances.” There will be three product 
streams: D, SS, and B. The property option to be used is Chao-Seader. The component 
list and stream specifications are the same as in Workshop 3.1a. Make the following 
specifications: 
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TABLE 3.1 Stream Specifications for Workshop 3.1a 


Stream/ 

Component 

Temperature 

(°F) 

Pressure 

(psi) 

Flow Rate 
(lbmol/hr) 

Vapor 

Fraction 

1 

200 

300 



C3 



20 


NC4 



30 


NC5 



30 


NC6 



20 


2 

250 

400 



C3 



30 


NC4 



30 


NC5 



20 


NC6 



20 


3 

200 



0.5 

C3 



50 


NC4 



0 


NC5 



30 


NC6 



20 



1. 50% of the combined feeds are to be directed to stream D. 

2. 50 lbmol/hr of NC4 will exit in stream SS. 

Solve the flowsheet. Is something wrong? What? Why? 

Workshop 3.2b Change Workshop 3.2a such that 20 lbmol/hr of NC4 exits in stream 
SS. Solve the flowsheet. What do you observe? 

Workshop 3.3a A two-feed, two-product distillation column operating at 1 atmo¬ 
sphere is to be simulated using the Sep2 block. The system components are methanol, 
water, and normal hexane. You may assume that the system thermodynamics can be 
represented by an ideal vapor phase and the Unifac correlation for the liquid phase. 
The componential flow rates of the first feed, a saturated liquid at 14.696 psi, are 50 
lbmol/hr of methanol and 100 lbmol/hr of water. The second feed flow rate is 150 
lbmol/hr of pure n -hexane at 120°F and 14.696 psi. Both feeds are at 1 atmosphere 
and are saturated liquids. The column performance specifications are as follows: 

1. One-half of the methanol fed is to be recovered in the overhead product. 

2. 99% of the water fed is recovered in the overhead product. 

3. 99% of the n -hexane fed is recovered in the bottom product. 

Set up the problem and solve. 

Workshop 3.3b Specifications 2 and 3 above are to be changed as follows: 

2. The mole fraction of water in the overhead product is 0.788. 

3. The mole fraction of n -hexane in the bottom product is 0.85. 


Modify Workshop 3.3a and solve. 
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Workshop 3.3c Specification 2 is to be changed as follows: 

2. The mole fraction of water in the overhead product is 0.8. 
Modify Workshop 3.3b and solve by hand. What do you observe? 
Workshop 3.3d Specification 2 is to be changed as follows: 

2. The mole fraction of water in the overhead product is 0.75. 
Modify Workshop 3.3c and solve. What do you observe? 


Workshop Notes 

Workshop 3.1b Care must be taken with the specifications of the allowable phases 
in the product. If there is a possibility that the resulting state is a mixture of phases, it 
must be specified on the Flash Options/Valid phases data-entry location. 


Workshop 3.2a The specifications stated in the problem have not been achieved. It 
is recommended that the problem be solved manually. This will show that a negative 
flow of one of the products results. This problem is specified improperly. 


Workshop 3.2c A satisfactory solution is 

Workshop 3.3a A satisfactory solution is 

Workshop 3.3b A satisfactory solution is 

Workshop 3.3c A satisfactory solution is 
below using the following nomenclature: 


obtained. 

obtained. 

obtained. 

not obtained. A manual material is given 


/i, f 2 , are the combined feed flow rates of methanol, water, and n -hexane, 
respectively. 

d \, d 2 , d 2 are the componential distillates. 
b\, b 2 , b 2 are the componential bottoms. 


The material balance equations are then as follows: 


d\ T b\ — 50 


d 2 + b 2 = 100 


ds + b 2 

h 

b\+b 2 + b 3 
d 2 


= 150 
= 0.85 

= 0.8 


d\ + d 2 + d 2 
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The first specification results in d\ = b\ = 25. Simultaneous solution of the remain¬ 
ing four equations and substituting the known values for d\ and b\ results in the follow¬ 
ing results: b^ = 1.59; b 3 = 151.4; d 2 = 98.46; d?> = —0.38. Since there is a negative 
flow for ^ 3 , ft-hexane in the distillate, this is a physically unrealistic solution; therefore, 
the specifications are improper. 

Workshop 3.3d A satisfactory solution is obtained as a result of the change in the 
overhead specification for water. 


CHAPTER FOUR 


PROCESSES WITH RECYCLE 


Virtually all Aspen Plus blocks are designed to operate as part of a sequential-modular 
(SM) simulator; that is, given the block’s inputs and operating parameters, the block 
will calculate internal conditions and its outputs. In most cases the block’s operating 
pressure and product stream pressures are specified in some way. This is illustrated by 
the example of a Flash2 block shown in Figure 4.1. 

Streams 1 and 2 are either specified or calculated by a block executed prior to block 
Bl. Stream 5 is a special case of free water, usually zero. The block’s products, each 
of streams 3 and 4, have n c + 2 variables, where n c is the number of components. 
These are the flows for each of the components temperature and pressure, for a total of 
2n c + 4 variables. Counting equations for a Flash2 block, there are n c material balances, 
n c equilibrium equations, and two equations that satisfy the equilibrium requirement 
that the temperature and pressure of the equilibrium streams, 3 and 4, are equal, for 
a total of 2n c + 2 equations with two degrees of freedom. The specification of any 
of the possible two flash conditions—for example, the heat added and the vapor/feed 
ratio—satisfies the block’s requirements and permits calculation of the state of both 
product streams. A process composed solely of SM models simulates the performance 
of a process. 

A variation on such a model poses the question: What temperature should the flash 
operate at to give a desired product composition, within the constraints of thermody¬ 
namic reality? As an example, if flashing a mixture of hydrogen and rocks, the basic 
laws of nature will prohibit rocks from being in the overhead product. This mode 
of simulation encompasses a design function. To accomplish this, one may employ a 
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Figure 4.1 Flash2 block. 


SM model by releasing a specification on the flash conditions and, instead, specify¬ 
ing the composition desired for one component in a product. This is termed a design 
specification. Details are described in Chapter Five. 

4.1 BLOCKS WITH RECYCLE 

When blocks are connected for a process with recycle such as Figure 4.2, traditional 
chemical engineering practice identifies stream 6 as a recycle stream. The technique 
for solution involves a guess of the component flows, temperature, and pressure of the 
stream and uses the method of direct iteration. The values of the variables in stream 
6 are estimated and the blocks are executed in the sequence React, Flash, and Split, 
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and produce a calculated set of values for stream 6. The values of stream 6 guessed 
and calculated are compared, and if they are not equal within a set tolerance, the cal¬ 
culated values of stream 6 are used as the set of values for the next iteration in the 
calculations. The method of direct iteration may require many iterations to converge 
and does not guarantee a solution. If a particularly difficult process is encountered, it is 
sometimes beneficial to select the convergence option, direct, and limit calculations to 
one or two iterations. The maximum number of iterations can be selected under Con- 
vOptions/Methods under the tab corresponding to the convergence method, and chosen 
as Maximum flowsheet evaluations. This frequently helps in analysis and debugging 
of models. 

The choice of stream 6 as a recycle stream is intuitive, but it is not the only approach. 
If one were to choose stream 2 as a “forward” recycle, one could still complete the 
calculations by using the sequence Flash, Split, and React to solve the process. Any 
stream selected as a “forward or backward recycle” is termed a tear stream. Aspen 
Plus has a facility to enable the user to select tear streams and to choose the sequence 
of calculations. These options can be found under Convergence, Tear, Convergence, 
ConvOrder, and Convergence, Sequence. Whether Aspen Plus or the user selects the 
tear streams, it is necessary to initiate the calculations with a set of starting values. 
Aspen Plus provides zeros, which in many situations will lead to a converged solution. 
If it does not, the user must provide reasonable estimates by entering values on the 
stream input form. Prior to rerunning the model, the calculations should be reinitialized 
by selecting the Run menu and clicking on Reinitialize. 

As an illustration of the method of direct iteration, consider solution of the 
equation 


f(x)=x 2 -c (4.1) 

Add the quantity v to both sides of the equation to form the function Fix), 

F(x) = x 2 — c + x = f{x) + v (4.2) 

and look for a solution where 


F{x) = x (4.3) 

which occurs when fix) = 0. Repeated iterations yield the solution desired. The 
derivative F\x) can be evaluated, numerically, from two successive iterations. The 
process converges under the conditions 

0 < F\x) < 1 (4.4) 

-1 < F\x) < 0 (4.5) 

Aspen Plus provides several other convergence options which may be selected by 
opening Data, Convergence, ConvOptions, Defaults, and the tab Default Methods. 
The default method (Wegstein, 1958) is usually suitable for most applications. Weg- 
stein’s method provides substantial robustness by making use of the first derivative 
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and modifying the next trial value of x 2 as follows: 

T = 1 - F'Otj) 
x 2 = (1 - r)vi + TF(x i) 


(4.6) 

(4.7) 


The direct iteration and Wegstein’s methods treat each of the state elements of a 
stream as if it were independent of the others, and hence for a single, n -component tear 
stream there are n + 1 equations in either method. Some of these ideas are illustrated 
in Figure 4.3, due to Rubin (1964), which has been modified minimally to facilitate 
the use of Aspen Plus’s blocks. Using the traditional approach, the flowsheet appears 
to have five tear streams, 1, 3, 4, 6, and 7, using the calculation sequence Mix, B2, B3, 
B4, SP1, B5, and SP2. Selection of these streams as tears is illustrated as follows. The 
selection Convergence under the setup drop-down menu opens the object manager, and 
when New is selected, Figure 4.4 appears. It is then possible to enter the desired tear 
streams. Aspen Plus detects whether the sequence is feasible and issues error messages 
if it is not. 

Through visual trial and error the reader should be able to discover the sequence B3, 
B4, SP1, B5, SP2, Mix, and B2 with tear streams 2 and 5, thus reducing the number 
of tear streams from five to two and potentially reducing the calculations considerably. 
As part of Rubin’s paper he developed an algorithm that can produce the minimum 
number of tears required for an arbitrary flowsheet. Such an algorithm is a part of Aspen 
Plus, and when this problem was solved with Aspen Plus, the optimal sequence was 
identified. This problem is available for study at Chapter Four Examples/Rubin.bkp. 

The well-known Newton method is an alternative method for solving process flow¬ 
sheets. Details of the method are given in Section 6.5. A significant difference between 
Newton and Wegstein is the formulation of equations that are to be solved. Each 
equation still has the form of equation (4.2), except that with Newton’s method each 
of the 2 n c + 1 equations is a function of the 2 n c + 1 unknowns associated with that tear 
stream. For example, if there are n t tear streams, n t (2n c + 1) unknowns and equations 
are to be solved. The Newton method requires partial derivatives of each equation with 
respect to each variable for each iteration. To create these numerically, after Aspen Plus 
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Figure 4.3 Rubin’s flowsheet modified. 

























































HEURISTICS 53 


EEI-Gfl Setup 

© Specifications 
© Simulation Options 
| © Stream Class 
EC1 Substreams 
| © Costing Options 
B D Stream Price 
El Units-Sets 
| © Custom Units 
© Report Options 
EE-Gfl Components 
E Properties 
ffi-Sl Flowsheet 
E Streams 
ffl-Sl Blocks 
Cl Utilities 
( D Reactions 
EHSfl Convergence 
B Gfl Conv Options 
E EO Conv Options 
O Tear 

B-SI Convergence 
! B <31 C-l 

|© Input 
□ Results 


^Tear Streams Calculator Tears | Parameters J 


T ear streams 



Stream 

T olerance 

Trace 

State variables 

Component group 

► 

1 

0.0001 


Pressure St enthalpy 



3 

0.0001 


Pressure St enthalpy 



4 

0.0001 


Pressure St enthalpy 



8 

0.0001 


Pressure St enthalpy 



7 

0.0001 


Pressure St enthalpy 


* 







Figure 4.4 Specifying tear streams. 


has an unconverged trial solution, each variable is in turn perturbed by a small amount; 
thus, n t (2n + 1) passes through the flowsheet are made. Convergence is quadratic. 

The Broyden (1965) method is a variant of Newton’s which estimates the values of 
the partial derivatives from the current values and may require fewer passes through 
the flowsheet. Convergence is linear. 

4.2 HEURISTICS 

1. Prior to attacking a material and energy balance problem, attempt a solution with 
simple blocks without energy balances. 

2. Solve easy variations of a problem prior to attempting the actual problem (i.e., 
do not use tight specifications). 

3. If errors are found in the solution and corrected, for the next run be sure to 
reinitialize the flowsheet; otherwise, Aspen Plus will continue the solution using 
the current, incorrect, values when it executes, which may prevent convergence. 

4. For a large flowsheet or one with a complex configuration in which the tear 
streams are dependent on each other: 

a. If it fails to converge try solving with Aspen Plus’s defaults. 

b. Break the flowsheet into smaller segments and converge them separately prior 
to integrating them into the complete configuration. 
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4.3 WORKSHOPS 

Workshop 4.1 The separation of ethanol from water by heterogeneous azeotropic 
distillation is carried out by two distillation columns, each with its own overhead 
condensor but sharing a common decanter, as shown in Figure 4.5. The first col¬ 
umn has two feeds, one of fresh material and one containing the water-rich phase 
recycled from the decanter. The bottoms product of this column is primarily water. 
The fresh ethanol-rich feed is composed of 1000 lb/hr of ethanol and 9000 lb/hr 
of water. The overhead product of this column is condensed and sent to the shared 
decanter. 

The decanter has a fresh cyclohexane feed of 1 lb/hr, which represents the makeup 
of cyclohexane lost in the products of both columns. The feed to the second column is 
the organic product of the decanter, which contains primarily ethanol and cyclohexane. 
The bottoms product of this column is primarily ethanol. The overheads product is 
condensed and sent to the shared decanter. All fresh feeds are saturated liquids. All 
operations are at 1 atmosphere pressure with negligible pressure drops within the 
equipment. The two distillation columns are to be simulated with Sep2 blocks and 
the decanter with a Sep block. Solve using the material-balance-only option and use 
Wegstein’s default method for convergence. 

The performance specifications of the unit operations are given in Table 4.1. For 
both columns the fraction of the feed of each component leaving with the respective 
bottoms product is given. For the decanter, the fraction of the feed of each component 
leaving with the organic phase is given. 

a. Why may one assume that the system can be represented by ideal thermodynam¬ 
ics for the vapor and liquid phases? 

b. Has Aspen Plus chosen the best tear streams? 



Figure 4.5 Heterogeneous azeotropic distillation. 
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TABLE 4.1 

Process Performance Data 



Component 

Fraction of Feed to Bottoms 

Fraction of Feed 
to Organic Phase: 
Decanter 

Column 1 

Column 2 

Ethanol 

0.01 

0.97 

0.98 

Water 

0.97 

0.0001 

0.01 

Cyclohexane 

0.09 

0.0001 

0.99 


Workshop 4.2 Increase the number of Wegstein iterations to 100. Reinitialize the 
problem, and solve. What do you observe? 

Workshop 4.3 Select C and F as the recycle streams. Provide initial values and 
reinitialize. 

a. Solve using Wegstein’s method. What do you observe? 

b. Solve using Newton’s method. What do you observe? 

Workshop 4.4 Figure 4.6 shows a continuous process for the recovery of solvents 
from a mixture that includes a soluble solid. Upon concentration, the solid precipitates 
and can be removed by filtration. You are to perform preliminary calculations without 
energy balances using the process description given below. 

The primary feed to the process, FD1, is composed of 100 lb/hr of ethanol (ETOH), 
750 lb/hr of water (WATER), 100 lb/hr of toluene (TOL), and 50 lb/hr of soluble 
solids (SOLIDS) at 14.696 psi. FD1 is sent to a column (CRUDSTL) along with a 
second feed, FD2, composed of 250 lb/hr of saturated steam (H 2 O) at 30 psig. The still 
produces two products: SOLVNTS and CRSOLIDS. 95% of the ethanol fed, 99% of 
the toluene fed, and 50% of the water fed end up in the stream SOLVNTS, and all of 
the solids end up in the stream CRSOLIDS. All of the solids in this stream precipitate 
out of solution. 

The stream CRSOLIDS is the feed to a filter that produces two streams: AQSLVNT 
and SOLIDS. All of the solids fed, 1.0% of the ethanol fed, 1.0% of the toluene fed, 
and 0.5% of the water fed end up in the stream SOLIDS. 

The stream SOLVNTS is fed to a column (TOLSTILL) whose role is to sepa¬ 
rate ethanol and toluene. The overheads (CRUDETOH) from this still are a stream 
that is primarily ethanol and water which is purified in the azeotropic distillation 
columns ESTILLHI and ESTILLLO. The bottoms (CRUDTOL) are sent to a decanter 
(DECANT). The performance characteristics of TOLSTILL are as follows: 25% of the 
ethanol fed, 45% of the water fed, and 99% of the toluene fed end up in the stream 
CRUDTOL. 

The decanter (DECANT) produces a relatively pure toluene product (TOL) and an 
aqueous product (AQTOL), which is recycled to the column TOLSTILL. The perfor¬ 
mance characteristics of the decanter are as follows: 50% of the ethanol fed, 99.9% of 
the water fed, and 0.1% of the toluene fed end up in the stream AQTOL. 

The stream CRUDETOH feeds the column ESTILLHI. The column produces two 
products: ETAZH, whose composition is to be 88 mol % ethanol, and ETAQ, whose 
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Figure 4.6 Distillation train. 


composition is to be 99% water. All of the toluene fed ends up in the stream ETAQ. 
The stream ETAZH feeds the column ESTILLLO. This column produces the relatively 
pure product ETH, which is to contain 99 mol % ethanol, and the product ETAZL, 
which is to contain 82 mol % ethanol. The stream ETAZL is recycled to column 
ESTILLHI. 

A steam stripper (STMSTRIP) is provided to remove solvents from the aqueous 
products of the process, ETAQ and AQSLVNT. The stream FD3 consists of 50 lb/hr 
of 30-psig steam that is used as the energy source for the stripping the solvents. The 
stripper produces the stream RSLVNT, which contains the stripped solvents and is 
recycled to the column CRUDSTL. The stream WATER is discharged to the envi¬ 
ronment. The performance characteristics of the stripper are as follows: 0.5% of the 
ethanol fed, 90% of the water fed, and 0.01% of the toluene fed end up in the stream 
WATER. 

a. There is no component SOLIDS. Devise a method that will permit this component 
to be used in the simulation. 

b. To calculate the material balance for the complete process flow sheet, it may 
be advisable to break up the process into manageable segments that are to be 
combined when the complete material balance is solved. Select the segments and 
solve each using an estimate for the required feeds when necessary. 
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c. Solve the complete process. Take care to place starting values for the tear streams 
estimated from the segments. 

Workshop 4.5 Search Aspen Plus Help for adding a stream table to your results. 
Replace all simple boxes with suitable icons. Place a stream table on your flowsheet 
graphic. Place suitable text labels on your flowsheet. 


Workshop Notes 

Workshop 4.1 

a. Since this is a material-balance-only approach, no rigorous blocks are used, and 
all internal calculations are by simple blocks, the type of thermodynamics chosen 
is arbitrary since it is not used in the simulation. If this problem were solved with 
rigorous blocks such as Radfrac and Decant, rigorous thermodynamics would be 
required. 

b. The Aspen Plus execution does not converge. You may try the following: 

• Increase the number of iterations. Check the overall material balance. 

• Check the performance of each block. 

• Check that no unrealistic specifications were used. 

• Change the tear streams. 

• Change the convergence method used. 

Workshop 4.2 Four-2.bkp uses 100 Wegstein iterations and it, too, fails to converge. 
This looks like a situation where Wegstein may not work or the tear streams Aspen 
chose were not suitable. Examining this process from a chemical engineering point of 
view, the composition of the streams that are easiest to estimate are those leaving the 
decanter: streams C and F. 

Workshop 4.3 Streams C and F may be initialized with unconverged results from 
Four-2.bkp. As an alternative, the compositions of streams C and F can be estimated 
from a phase diagram of the ethanol-water-cyclohexane system. 

a. Four-3.bkp using Wegstein’s method again fails to converge. This problem is not 
suitable for Wegstein convergence, even with 100 iterations. 

b. Four-3n.bkp using Newton’s method converges easily with only a few iterations. 

Workshop 4.4 

a. Since Aspen Plus deals mostly with solids in a special way (an advanced appli¬ 
cation) and the states considered are either vapor or liquid, there is no specific 
provision for dealing with solids dissolved or suspended in a solution. Further¬ 
more, for a situation in which there are no specific database components for 
unidentified by-products of reactions, one must assure that these are accounted 
for on a mass basis, and that their behavior in rigorous blocks is accounted for 
in some way. For example, it is safe to assume that the unknown vapor pressure 
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of a solid by-product is much lower than that of the other components; there¬ 
fore, none of that species will appear in a vapor phase. A simple method to 
accomplish this is to select a very heavy component, such as eicosane, C 20 H 42 , 
from the database and change its name, for example, to HEAVY, and ignore the 
properties that are calculated when that species is part of a mixture, 
b. Select the sequences 

• CRUDSTL, FILTER, STMSTRP 

• TOLSTILL, DECANT 

• Bl, B2 

The complete solution may be found at Four-4.bkp. Note that if this were to be 
solved rigorously, many problems would be raised by the complex thermodynamics 
required. This problem is addressed rigorously in Chapter Fourteen. 

Workshop 4.5 A solution to the flowsheet is required. Under the View menu, the 
PFD mode must be checked. 
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CHAPTER FIVE 


FLOWSHEETING AND MODEL 
ANALYSIS TOOLS 


Two of the most useful tools in Aspen Plus are the Sensitivity study and the Design 
specification, both located on the Data Browser. These tools work hand-in-hand and 
are used when a design specification is used to specify a block output variable. For 
example, the file Chapter Five Examples/mixer5.bkp presents a mixer block with three 
feeds. When all input streams are specified properly, the block will calculate the com¬ 
position and state of the output stream. Figures 5.1 and 5.2 show the graphic setup 
and the solution. Note that the product, stream 4, is a straightforward summation on a 
componential basis, of all the feeds, excluding an energy balance. Suppose, however, 
that in the example one requires that the composition of a specific component, say 
ethanol, is to be 0.5 mole fraction. It will be necessary to increase the flow of stream 
2. This is easy to do with a hand calculation, but if this block were in the middle of 
a process flow diagram, it might not be a trivial task. The Sensitivity study permits 
the specification of a trial range of flow rates for stream 2, which may be plotted 
against the resulting ethanol composition. Once an approximate solution is obtained, 
by plotting the results and checking the location of the target composition, the Design 
specification can be employed to provide a means for Aspen Plus to solve for the 
exact value of the flow rate of stream 2, which satisfies the specification desired. In 
setting up both the Sensitivity study and the Design specification, it may be necessary 
to employ the Fortran language capability embedded in several blocks and available 
on a stand-alone basis in the Calculator block. 

5.1 INTRODUCTION TO FORTRAN IN ASPEN PLUS 

Fortran (Formula Translation) is a very old but rich language which has traditionally 
been used for scientific applications. Details of the language are very complex and can 
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Figure 5.1 Mixer5 example. 



|i 3 

I 2 Jll 

I 3 d 

I 4 d 

Substream: MIXED 





Mole Flow Ibmol/hr 





METHANOL 

100.0000 

0.0 

0.0 

100.0000 

ETHANOL 

0.0 

100.0000 

0.0 

100.0000 

WATER 

0.0 

0.0 

100.0000 

100.0000 

Mole Frac 





METHANOL 

1.000000 

0.0 

0.0 

.3333333 

ETHANOL 

0.0 

1.000000 

0.0 

.3333333 

WATER 

0.0 

0.0 

1.000000 

.3333333 

T otal Flow Ibmol/hr 

100.0000 

100.0000 

100.0000 

300.0000 

T otal Flow Ib/hr 

3204.216 

4606.904 

1801.528 

8612.648 

T otal Flow 





T emperature F 

70.00000 

70.00000 

70.00000 


Pressure psia 

14.70000 

14.70000 

14.70000 

14.70000 


Figure 5.2 Mixer5 solution. 


be found in many books. Aspen Plus employs both interpreted Fortran (line-by-line 
execution) for the blocks in this chapter and compiled Fortran (dynamic link library 
creation) for the creation of user models. An overview of the language criteria can be 
found in Aspen Plus Help and by searching About the Interpreter. Aspen Plus user 
model documentation recommends use of the Intel Fortran compiler. 

5.2 BASIC INTERPRETED FORTRAN CAPABILITIES 

For the purposes of this chapter’s functions, what follows is the primary function¬ 
ality that is available when using interpreted Fortran. There is no intent here to 
present a thorough documentation of the Fortran programming language, merely what is 
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sufficient to be able to work with the subject matter of this chapter and some material 
in later chapters. 

5.2.1 Primary Fortran Operators 

The primary Fortran operators are given in Table 5.1. Prior to execution of any state¬ 
ment, all variables appearing on the right-hand side of an equal sign must be defined: 
that is, contain either a numerical value or a value calculated with previous legal 
statements. 

A list of all the functions available is given in the Aspen Plus documentation About 
the Interpreter. 

5.2.2 Precedence of Calculations 

The following rules apply when a Fortran statement is evaluated: 

• Fortran statements may be organized in groups of operations within parentheses. 

• The contents of each group is evaluated prior to evaluation of the complete state¬ 
ment. 

• Function calls precede group evaluation. 

• Exponentiation precedes group evaluation. 

• If part of a group, divisors are evaluated prior to numerators. 

For example, the precedence of calculations for the Fortran statement 

(f * dlog(b) + c) 

a =- (5.1) 

(a + b) + (c — d * e ** f) v ' 


is as follows: 


TABLE 5.1 Fortran Operators 


Operator 

Definition 

Example 

Result 

= 

replacement 

A = B 

B 

+ 

addition 

A = A + A 

2A 

- 

subtraction 

A = A — A 

0 

/ 

division 

A = A/A 

1 

* 

multiplication 

A = A* A 

A 2 

** 

exponentiation 

A = A**B 

A b 

.LT. 

.GT. 

.EQ. 

..AND. 

logical less than 
logical greater than 
logical equal 
logical and 
function call 
subroutine call 

A = Dlog(B) 

CALL XYZ 

Ln(B) 

XYZ executes 
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1 e ** f 

2 . d*e**f 

3. (c-d*e**f) 

4. (a+b) 

5. (a+b)+(c-d*e**f) 

6 . dlog(b) 

7. f* dlog(b) 

8 . (f*dlog(b)+c) 

9. (f*dlog(b)+c)/(a+b)+(c-d*e**f) 

5.2.3 Statement Format 

Intel Fortran, as implemented in Aspen Plus’s interpreter, is not quite standard (there 
are differences of opinion as to what constitutes standard Fortran). The main formatting 
rules, taken from Aspen Plus’s documentation, are given below. 

• Column 1 can contain a C. All other entries on the same line are comments. 

• Column 2 is left blank. 

• Columns 3 to 5 are reserved for statement labels. Under certain logical conditions, 
program execution can be transferred to a labeled statement. 

• Executable statements begin at column 7 and beyond. 

• No variables may begin with the characters IZ or ZZ. 

• Integer variables begin with the letters I through N. 

• Real variables begin with the letters A through H or O through Z. 

• Variables names are limited in length to seven characters. 

• Lowercase characters are permitted. 


5.2.4 Program Logic Control 

Logical statements may take a variety of forms. An example of one of the IF variations, 
the most common, is as follows: 

IF ((expression a). EQ. (expression b)) (expression c) 
expression d 

Expression c will execute if the logical grouping is true; otherwise, expression d exe¬ 
cutes. In all statements of this type, the number of left parentheses must equal the 
number of right parentheses. The logical operators are all interchangeable and can be 
nested within each expression to meet program requirements. 

Another common logic control statement is the direct Go To: 

Go To statement label 

All statements following the Go To are skipped until the labeled statement is found: 
then execution proceeds normally. 
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5.3 SENSITIVITY FUNCTION 

The features of the sensitivity function will be displayed by the continuation of example 
mixer5.bkp through example mixer5s.bkp. The Sensitivity function is located under 
the setup menu Model Analysis Tools and when selected initiates the Sensitivity Data 
Browser. When New is selected, a Sensitivity ID such as S-l is assigned. This initiates 
a display that facilitates the association of a Fortran-accessible variable with a variable 
within the process. It is imperative that this process be started by selecting the New 
button, which generates the displays shown in Figure 5.3. The variable xetoh, which is 
user selected, represents the mole fraction of ethanol in output stream 4 of the example. 
Figure 5.4 shows the assignment of the process variable to the Fortran variable. For 
this example, a Fortran expression will be created to calculate the fraction of alcohols 
in stream 4. For this purpose it will be necessary to create three additional Fortran- 
accessible variables (fmoh, fetoh, and f4), the flow rates of methanol and ethanol in 
stream 4, and the total flow of stream 4. 

Selection of the tab Vary in the Sensitivity input Data Browser brings up Figure 5.5, 
which defines the range of trial values (the independent variable) that stream 2 should 
use in seeking for a value that would produce an approximate solution to the desired 
target of ethanol mole fraction in stream 4 of 0.5. 

Selection of the tab Tabulate permits the assignment of any defined Fortran variables 
or legitimate Fortran expression to the column of a table on which the results of the 
Sensitivity study wifi be presented. Note that column 3 contains an expression for the 
calculation of the fraction of alcohols leaving the mixer. Tabular and plotted results are 
displayed in Figure 5.6. One can see that the flow rate required to produce a 0.5 mole 
fraction of ethanol is bracketed between a stream 2 flow of 100 and 300 mol/hr, 
and from the plot the solution is approximately 225 mol/hr. Note that the mole 
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Figure 5.3 Creation of a Fortran-accesible variable. 
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Figure 5.4 Association of flowsheet variable with Fortran variable. 
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Figure 5.5 Defining the range of a trial variable. 
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Figure 5.6 Results of sensitivity study Mixer5s. 


fraction of alcohols in stream 4 is tabulated. This example may be found at Examples/ 
mixer5s. 

5.4 DESIGN SPECIFICATION 

The features of the design function will be displayed by the modification of the example 
mixer5s.bkp and is given at mixer5d.bkp. The Design-Spec is located under the setup 
menu Flowsheeting Options, and when selected initiates the Design-Spec Data Browser. 
Selection of the New button assigns an identifier to a design spec such as DS-1. The 
key to implementing the design specification is to manipulate the feed, stream 2, so that 
it can take on the required flow rate to achieve the required stream 4 specification. One 
possibility is shown in Figure 5.7. The inclusion of stream 5, which has a larger flow 
than what is required to meet the stream 4 specification (the flow rate is determined 
from the sensitivity study), is fed to a multiply block. The design spec function will 
manipulate the multiplier value, which will place the appropriate value into stream 2 to 
satisfy the design specification. The implementation requires the definition of a Fortran 
variable, in this case, XETOH, just like the setup in the Sensitivity study; however, 
there are two other tabs in the Design-Spec Data Browser, the target value and an 
acceptable tolerance, which must be selected and implemented. Figure 5.8 shows the 
setup of the function to be solved: 

/(factor) = xetoh — 0.5 (5.2) 

where factor is the multiplying factor of the block Ml. The function is solved by 
Aspen Plus’s built-in secant method equation solver within a user-specified tolerance, 
in this case 0.001. Results are shown in Figure 5.9. 
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Figure 5.7 Flowsheet for Mixer5c. 
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Figure 5.8 Mixer5d design specification. 
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Figure 5.9 Mixer5d design specification satisfied. 
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5.5 CALCULATOR FUNCTION 

The calculator may be found under the Setup menu Flowsheeting options. Like the 
Design and Sensitivity functions, the user associates flowsheet variables with user- 
defined Fortran variables. The basic calculator display is shown in Figure 5.10. The 
association of a Fortran variable occurs after the New button is clicked, but requires 
that the line upon which the new variable is entered be selected and the Edit button 
clicked. The entries for the creation of the variable FL02 are shown in Figure 5.11. 
Two very important selections are Import, which copies current values of variables 
from the flowsheet into the Fortran variables, and Export, which copies values of 
the Fortran variables into current variables of the flowsheet. An example of Fortran 
coding that accomplishes the same result as in the example Mix5d.bkp is given in 
Figure 5.12, which uses the definitions in Figure 5.10. Figure 5.13 illustrates the speci¬ 
fications, which permit placement of the calculation virtually anywhere in the sequence 
of flowsheet calculations. 

The results of the calculations are shown in Figure 5.14. The results were accom¬ 
plished with the first Fortran statement, a material balance to calculate the amount of 
ethanol that must be added to stream 4 to make the 0.5 mole fraction specification. 
The second statement exports the sum of the original value of ethanol in stream 4 and 
the results of the first statement. When these operations are completed, Aspen Plus 
issues a series of errors and warnings which are documented in the control panel as 
shown in Figure 5.15. This example may be found at Examples/mixer5c. Attention 
should be paid to the information in the control panel because it will caution the user 
to any irregularity. As an example, note the warning that the mixer block is not in 
material balance. This occurred because the mixer was executed prior to the calculation 
function. 
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Figure 5.10 Basic calculator display with Fortran variables defined. 
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Figure 5.11 Association of flowsheet variable with Fortran variable. 
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Figure 5.12 Fortran coding for Mixer5c equivalent to Mixer5d. 
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Figure 5.13 Sequence control of Fortran code. 
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Figure 5.14 Results of Mixer5c. 


5.6 TRANSFER FUNCTION 

The error “mixer in the example mixer5c.bkp is not in material balance” occurs because 
stream 2 has not been updated with the contents of stream 4, as shown in Figure 5.14. 
Aspen Plus provides the Transfer function, which enables the copying of data from 
stream to stream. When the Transfer function is invoked, an ID is assigned and the 
display shown in Figure 5.16 is displayed. Selection of tabs enables the selection of 
streams to copy from and to as well as which elements of those streams to select. The 
tab sequence enables the user to choose the order within the flowsheet calculation to 
run the function. Figure 5.17 shows the results of the transfer. 
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Figure 5.15 Warnings and errors in Mixer5c results. 
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Figure 5.16 Primary transfer function display. 
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Figure 5.17 Results of Mixer5t. 


5.7 WORKSHOPS 

Workshop 5.1 Copy the file four3n.bkp, which will be the basis for using Aspen 
Plus’s Sensitivity Studies capability to determine the effect of the decanter’s perfor¬ 
mance with regard to the fraction of cyclohexane leaving the decanter in its aqueous 
product, stream C. For the base case this specification is 0.02. The Sensitivity Studies 
should vary this value between 0.01 and 0.10 in increments of 0.01. Select as tear 
streams C and F. Make a reasonable estimate of the componential flows and solve the 
flowsheet using the Newton convergence option. This becomes the base case solution. 
This file is to be modified using the Sensitivity function. Select the following as the 
dependent variables: 

1. The total flow of stream G, the cyclohexane-rich distillate recycled to the 
decanter. Call it FLOWG. 

2. The flow of cyclohexane in stream G. Call it FLOWGC. 

Plot the results. If the pipeline in which G is to flow is limited to 300 lb/hr, what is 
the fraction of the cyclohexane entering the decanter that leaves in stream C? 

Workshop 5.2 Using Workshop 5.1 as a basis for solving Workshop 5.1 by means 
of the Design Specification capability of Aspen Plus, delete the Sensitivity function and 
replace it with Design-Spec. Vary the decanter performance with regard to cyclohexane 
in stream C and use it as brackets for the iterations 0.01 and 0.10. Set up the design 
specification such that the total flow rate of stream G is 300 lb/hr with a tolerance of 
0.1 lb/hr. What is the fraction of the cyclohexane entering the decanter that leaves in 
stream C? 
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Workshop 5.3 Using Workshop 5.1 as a basis, create a Sensitivity study to determine 
the effect of feed composition on the concentration of ethanol in stream B and water 
in stream H. To do this it is necessary to change the quantity of ethanol in stream 
A while reducing the amount of water in stream A by the same amount for each 
case. Aspen Plus’s Sensitivity function will permit variation in the ethanol quantity for 
each case, but it will be necessary to write Fortran that executes prior to each time 
the Distillation I block executes. This capability will be found under “Flowsheeting 
Options—Calculator.” Use the following values of ethanol in A, in units of lb/hr, for 
your Sensitivity studies: 1000, 1100, 1500, 2000, and 5000. 

Workshop 5.4 Part of a plant-wide environmental control facility (ECF) is a system 
for diverting aqueous waste fluids from a discharge point such as a river in the event 
that the concentration of volatile organic compounds (VOCs) exceeds a certain level. 
This logic is to be modeled by employing a splitter with two product streams. If the 
composition of the VOCs exceed 20 parts per billion, the stream is to be diverted to the 
second outlet of the splitter. Otherwise, the stream is to be diverted to the first outlet. 
Execute a calculation block in front of the splitter to program the logic in Fortran. Test 
with the two potential feeds listed in Table 5.2. 


TABLE 5.2 ECF Feeds 


Component 

voc 

Feed 1 lb/hr 

Feed 2 lb/hr 

Methanol 

Yes 

1 

10 

Acetone 

Yes 

2 

20 

Toluene 

Yes 

3 

30 

Methylene chloride 

Yes 

4 

40 

Water 

No 

1.0 x 10 9 

1.0 x 10 9 


Workshop Notes 

Workshop 5.3 The Fortran statement required may be found in the file Workshops/ 
Five-3.txt. 

Workshop 5.4 The Fortran statements required may be found in the file Work- 
shops7Five-4a.txt. 
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CHAPTER SIX 


THE DATA REGRESSION SYSTEM 


Aspen Plus’s regression system can be used to fit various types of pure component 
physical property data such as vapor pressure data, or a parameter for a new component 
in the Unifac database, but its primary use is to customize thermodynamic models 
dealing with vapor-liquid equilibrium (VLE) and liquid-liquid equilibrium (LLE). 
Most of the equations of state and activity coefficient models commonly used are 
described in detail by Poling et al. (2000), as are listings of textbooks and data sources 
relating to fluid-phase equilibrium. Additional descriptions of models and variants may 
be found in Aspen Plus’s documentation. 

When an ideal system is modeled, the equations and/or correlations that are used to 
calculate the physical and thermodynamic properties of both phases are well defined 
and do not require augmentation with experimental data; however, the opposite is true 
when a real system is considered. In this regard there are two possibilities: 

1. When both phases of a system can be represented by an equation of state, the 
fundamental equation for vapor-liquid equilibrium for component i is given by 

=xi§[ (6.i) 

where, y and x are mole fractions of the vapor and liquid phases, respectively, 
and (\> is a fugacity coefficient which depends on temperature, pressure, and 
composition and is calculated from an equation of state. 

2. When the vapor phase can be represented by an equation of state and the liquid 
phase by an activity coefficient equation, the fundamental equation for vapor 
equilibrium is 


y t <\>iP = wifi 


( 6 . 2 ) 
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where P is pressure, c|> is calculated from an equation of state, y is an activity 
coefficient calculated from an activity coefficient equation, and f t ° is the standard 
state fugacity, which is given by 


f? = exp (yj at ^dp'j (6.3) 

The exponential term above, known as the Poynting correction , is used for very 
light components and is usually negligible. In most cases the vapor pressure, 
P sat , is used as the componential liquid-phase fugacity. 

The fundamental equation for liquid-liquid equilibrium equates two instances 
of the right-hand side of equation (6.2) written once for each liquid phase. 

Aspen Plus supplies property methods identified by a specific vapor-phase fugacity 
equation and a specific activity coefficient equation which defines the use of equation 
(6.1) or (6.2). For example, the method NRTL-2, which represents the nonrandom 
two-liquid activity coefficient equation, also known as the Renon equation, is used in 
conjunction with the ideal gas equation and employs equation (6.2) in calculations. 
The “2” identifies the data set from which the NRTL parameters were derived. Aspen 
Plus has the possibility of two sets of activity coefficient parameters for most activity 
coefficient equations with an indicated range of applicability (i.e., temperature and 
pressure). Henry’s law is used in some property methods for light components with 
Henry’s law parameters provided in the Aspen Plus database. Aspen Plus supplies 
more than 50 property methods. 


6.1 PARAMETERS OF EQUATIONS OF STATE 


The most commonly used nonideal equations of state available in Aspen Plus all have 
the same basic structure: third order. The general formulation from which all modern 
forms can be derived is the five-parameter form 


RT O(V-ri) 

V-b ~ (V -6)(V 2 + 8V + s) 


(6.4) 


The parameters a and b are functions of the critical properties and the acentric factor 
oo. The symbols 0, T], 8, and e have values that depend on the equation of state (chosen 
details may be found in Poling, 2000). For example, the values for pure components 
of the Peng-Robinson equation are 


0 = aa(T r )* (6.5) 

8 = 2b (6.6) 

e = - b 2 (6.7) 

a (Jr)* = [1 + (0.37464 + 1.5422 co - 0.2699 co 2 )(l - T r 0 - 5 )] 2 (6.8) 


For mixtures, the most commonly used equations of state are ideal, Redlich-Kwong, 
Soave-Redlich-Kwong with variants, and Peng-Robinson with variants. Most use 
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the mixing rules given by 


a = y j Y XjXj^djcij (1 — ki j) 
i i 

(6.9) 

W 

II 

(6.10) 


i 


Modifications of the mixing rules for some variants may be found in Aspen Plus 
documentation. The introduction of the k LJ parameter, one for each binary pair in a 
mixture, improves the ability of the equations of state to fit experimental data. These 
values, which are available from a variety of sources (e.g., Knapp et al., 1982), are 
included in the Aspen Plus database and may be used within the temperature and 
pressure ranges given in Aspen Plus’s documentation. There is rarely a need to fit the 
kij parameter. As an example, Figure 6.1 shows the kij parameters for the six binary 
pairs in the system nitrogen-methane-propane-ft-pentane. Note that the upper and 
lower bounds of the working temperature ranges are given. 

If the need arises to correlate a binary system not in the database, a set of PVT data 
with system vapor and liquid compositions is required. Aspen Plus implements the kij 
parameter with two terms: 

hj = aij + bijT ( 6 . 11 ) 

with the introduction of the b t j term for temperature dependency. Some variants have 
additional terms which may be selected during the setup of the regression. 



Figure 6.1 Soave-Redlich-Kwong Parameters. 
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6.2 PARAMETERS OF ACTIVITY COEFFICIENT EQUATIONS 

Activity coefficient equations are derived from the general formulation of excess free 
energy, 


RT lny* = 


dG E 

drii 


T,P,nj^i 


( 6 . 12 ) 


where is the number of moles of component i, and ri j^i indicates any component 
except i. Applying equation (6.12) to a model of g E , the total excess free energy per 
mole of mixture produces the activity coefficient equation. This is illustrated by the 
application of equation (6.12) to 



= —x\ \n(x\ + A\ 2 x) — X 2 ln (^2 + A 21 JC 1 ) 


(6.13) 


the two-component formulation of g E for the Wilson equation, and results in 


RT In Yi = —x\ ln(vi + A 12 *) — *2 


A12 


X x + A 12^2 


A 2 1 \ 

A2I-T1 + X2) 


(6.14) 


where the activity coefficient for component 1 is yi, and A 12 and A 21 are parameters 
which are determined by regressing experimental data. When used in multicomponent 
calculations, binary parameters representing all permutations of the components are 
used in the multicomponent form of the Wilson equation: 




+ '-£ 




(6.15) 


The Aspen Plus implementation of the Wilson parameter A z -y includes the user’s choice 
of terms for temperature dependency and is shown in 

l n Aij = ciij + + Qj In T dijT + (6.16) 

where in Aspen Plus’s terminology, A replaces A. 

When regressing data the user may select which parameters of equation (6.16) are 
to be regressed from the Properties/Regression/R-1/Input form parameter tab, shown 
in Figure 6.2. On each column of the central part of the display reading downward, 
the parameter is identified as binary, the name of the activity model used, the element 
number [which corresponds to the coefficients in equation (6.16), e.g., a t j is 1], and 
the components are identified in order; for example, ethanol is 1 and water is 2, so 
column 1 refers to parameter <212 and column 2 refers to parameter < 221 • The remainder 
of the input permits the user to omit, fix the value of, regress the parameter and assign 
an initial value, define the bounds of the regression, and set a scale factor such that all 
parameters have comparably scaled contributions to the regression. The last entry on 
the column is used for regression of equations of state. 
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Figure 6.2 Parameter Selections for Regression. 


The most commonly employed activity coefficient equations are the van Laar, 
Wilson, NRTL, Uniquac, and their variants. The input for all equations and the iden¬ 
tification of the parameters is virtually the same as in the example above. 

A very large collection of experimental data for many binary systems, from many 
sources, fit to the activity coefficient models is incorporated within Aspen Plus. An 
example of Wilson-2 parameters for an ethanol-water system is shown in Figure 6.3. 
These are available for user development of Aspen Plus simulations. Note that the tem¬ 
perature range of applicability is shown. These may not, however, satisfy requirements 
for many applications. For example, temperatures and pressures may not be appropri¬ 
ate to the modeling needs of specific applications. An important task for a user is to 
locate suitable data and customize (i.e., fit) the data to a model that lends itself to the 
application. If the required binary data are not available, the generalized correlation, 
Unifac, can be used to estimate the VLE data, which may then be fit to any of the 
activity coefficient equations. In such a situation the accuracy of a simulation may be 
compromised, and caution is advised. 

6.3 BASIC IDEAS OF REGRESSION 

Frequently, engineers and scientists are faced with sets of data that are proposed to fit 
a linear equation. The typical first pass at fitting the data is an eyeball plot on linear 
coordinates. This plot will usually differ when several persons attempt to plot the same 
data, and for the model 


y = ax + b 


( 6 . 17 ) 
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Figure 6.3 Wilson equation parameters. 


the values of the slope, a, and intercept, b , will differ. This problem is solved using 
the well-known method of least squares, which works in the folllowing manner. An 
objective function, 


rcdata 

T = X ^ calc “ flXdata “ ( 6 - 1 8) 

i 

is a measure of the goodness of fit and is defined as the sum of squares of errors 
between the individual values of the dependent values of each of the y data points 
observed and the values of y calculated from the linear model, with trial values of a 
and b. The minimum of this function can be found by setting 9\J r/da and 9\J r/db equal 
to zero. This generates a simultaneous set of linear equations which can be solved for 
the values of a and b that best fit the data. 

The problem of fitting data to the parameters of a nonlinear activity coefficient model 
is similar to the above. In Aspen Plus an objective function must be selected from 
several alternatives, given in Table 6.1. The maximum likelihood objective function is 
minimized directly subject to the constraints imposed by the applicable thermodynamic 
relationships using nonlinear programming methods such as those described by Edgar 
et al. (2001). All the remaining objective functions may be minimized by solving the 
nonlinear equations that arise by setting the partial derivatives of the objective function 
relative to the model parameters equal to zero. 
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TABLE 6.1 Objective Functions, Aspen Plus Version 7.0 


Objective Function 

Formed By: 

Maximum likelihood 

All measured variables 

Ordinary least squares 

Pressure and vapor compositions for isothermal VLE data 
Temperature and vapor compositions for isobaric VLE data 

Barker’s method 

Pressure only 

Modified Barker’s method 

Pressure and vapor compositions 

Activity coefficients 

Activity coefficients 

Equilibrium constants 

Equilibrium constants (K values) 

Relative volatility 

Relative volatility relative to first component 


The choice of objective function depends on the type are data that are available. For 
vapor-liquid equilibrium the data may be of the type T-P-x-y or T-P-x, where T 
is temperature, P is pressure, and x and y are liquid- and vapor-phase mole fractions, 
respectively. T-P-x-y data can be used with all objective functions except Barker’s 
method, which is used with T-P-x data only. Selection of an objective function may 
be found at Propereties/Regression by choosing New regression and the tab Algorithm. 

An example of an objective function suitable for use with T-P-x-y data is Activity 
Coefficients: 


miata ft com p / 0 b s calc ^ 

* = E E f ~ Yi 


v °b S / 

7 — 1 *‘ = 1 V Yi / 

where the experimental activity coefficients are calculated using 


bs = y^i(y{,y J 2 Pj,T j )P j 

Yi 


x i Py (T°j) 


(6.19) 


( 6 . 20 ) 


the experimental values and trial values of the activity coefficient parameters. Note the 
scaling of the activity coefficients in the objective function by the denominator, y° bs : 
When VLE data are of the type T-P-x the Barker objective function 

data 

i|/ = J2 ( p j - Pf c f (6-21) 

7 — 1 

is indicated. Here PJ alc is calculated using the trial values of the activity coefficient 
equation parameters and the experimental values 

pj alc = y[x{Pv{ + y[x[Pv[ (6.22) 

Fitting liquid-liquid equilibrium data (i.e., T -x 1 -* 11 , where x 1 and x 11 refer to the 
composition of the two liquid phases in equilibrium) can be done with the least squares 
or the maximum likelihood objective function. 
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6.4 MATHEMATICS OF REGRESSION 

Two approaches for regressing data are available in Aspen Plus, minimizing an objec¬ 
tive function by solving simultaneous nonlinear equations similar to linear regression 
and minimizing an objective function subject to the constraints of equilibrium. The 
partial derivatives of the objective functions with respect to the parameters are calcu¬ 
lated numerically by perturbing each of the parameters in turn by a small value, A a, 
a small fraction of the parameter, while keeping the remaining parameters constant. 
Using the perturbed and current values of the objective function, partial derivatives are 
calculated using 


+ ^<2 
da A a 


(6.23) 


6.4.1 Newton-Raphson Method for Solution of Nonlinear Equations 

The generalized Newton-Raphson method involves the arrangement of a set of n 
equations to be solved in the format 

fi(xi,x 2 , ...,x n )=0 


fn(*l,X 2, ...,x n )=0 

Each of the functions is expressed as a first-order Taylor’s series about a point to 
represent the functions at a point X^ +1 in the vicinity of X*: 


fi(X k+1 ) = fi(X k ) + ^ 

OX 1 


/ 2 (X* +1 ) = / 2 (X*) + 


dX] 


(• 4 +l -xb + JT 

X k dx 2 

(*f +1 -**> + §£ 

X k dx 2 


Xk 


Xk 


(x* +l -4) + ... + y± 

ox n 


x k 


K + -K) 


x k 


(x™-x k n ) 


f n (X k+l ) = /„(X*) + 


k, , df n 


dxi 


(x ^-4)+ d JL 

x k dx 2 


Xk 


(V2 + 1 — X%) + * * * + 


% 

dx n 


(4 +l -4) 


Xk 


Hypothesizing that all functions are zero at the point X^ +1 , evaluating all derivatives 
at the point X*, and changing the nomenclature such that 


A Xi = xf +1 - x- 


produces the following linear set of equations, which can be solved for the vector of 
Av. 
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The algorithm for solving a set of nonlinear simultaneous equations by the 
Newton-Raphson method is given below. If the equations to be solved are not 
analytical, it is necessary to calculate derivatives numerically. 

1. Guess X k . 

2. Evaluate all fi(X k ). 

3. If all \fi(X k )\ < s, the problem is solved. 

4. Calculate all partial derivatives at X k . 

5. Solve the matrix-vector equations 


PAX = -F for AX AX = P _1 (-F) 


where P is the matrix of partial derivatives, AX is the vector of x moves, and F 
is the vector of function values. 

6. For all x z -, calculate 


xf +1 = x\ + A Xi 


7. Return to step 2, substituting X i+I for X k . 

6.4.2 Direct Optimization of an Objective Function 

The default regression method in Aspen Plus is maximum likelihood. In this method 
the objective function is given by 



where the superscripts e and m refer to estimated and measured values. For each point 
in the set of data the following constraints apply. For simplicity, Poynting correction 
is ignored; then for each component one may write 


yn t>r p - Vi*/? = o 


(6.25) 

(6.26) 
(6.27) 


Y i = Y(T, x ) 

4>r =4> V (T P,y) 
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and 



1 =0 


(6.28) 




1 =0 


(6.29) 


Aspen Plus documentation does not provide information on the details of the imple¬ 
mentation of the minimization procedure; however, the papers of Anderson et al. (1976) 
and Fabries and Renon (1976) describe possibilities. 

6.5 PRACTICAL ASPECTS OF REGRESSION OF VLE OR LLE DATA 

The first question that arises when activity coefficient equations are employed is 
whether to use a two- or four-parameter approach. With the NRTL equation it would 
be either three or five. Aspen Plus recommends that the NRTL c parameter be fixed 
according to the following rules: 

• 0.2 for saturated hydrocarbons with polar nonassociated liquids and systems that 
exhibit liquid-liquid immiscibility 

• 0.3 for nonpolar substances; nonpolar with polar nonassociated liquids; small 
deviations from ideality 

• 0.47 for strongly self-associated substances with nonpolar substances 

When using two parameters, a question arises as to the selection of the constant a t j 
or the temperature-dependent bij version of the parameters. Prausnitz et al. (1999) 
discuss this point and no clear advice is available; however, experience has shown 
that the difference between the two approaches in quality of fit is very small. Often, 
the difference between the two- and four-parameter approach is also very small. When 
regressing four parameters simultaneously, steps taken during the regression are related 
to the current values. It is important to note that relative to the bij values, the a x j 
values vary sometimes by two orders of magnitude. Aspen Plus provides the ability 
to assign a scaling factor to any parameter so that step sizes during the regression are 
not dominated by a particular parameter. 

In most cases of nonlinear data regression, more than one set of values that satisfy 
the equations can be found. It is imperative that results be checked to assure that the 
results produce calculated values of the independent values that fit the data consistent 
with the requirements of the simulation for which the results will be employed. In 
many cases the values calculated by the regression will depend on the starting values 
of the parameters. In some cases a regression that employs Aspen Plus’s default values 
will suffice. 

6.5.1 Regression of VLE Data 

A useful technique for providing starting values for vapor-liquid equilibrium regression 
using T-P -x -y data is to estimate the infinite dilution activity coefficients by calcu¬ 
lating activity coefficients using equation (6.18), plotting versus x\ and extrapolating to 
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TABLE 6.2 

Infinite Dilution Activity Coefficients for Various Equations 

Equation 

In yf and In yf 

Van Laar 

Wilson 

NRTL 

Uniquac 

Inyf = Ai 2 ; lny£° = A 2 \ 

lny'j" 0 = 1 — In A 12 — A 21 ; lny“ = 1 — In A 21 — A 12 

lny'j’ 0 = T 21 + Ti 2 exp(— 0112 T 12 ); lny|° = t 12 + t 2 i exp(-ai 2 t 2 i) 

Inyj*° = ln(n/r 2 ) + [5\n(qir 2 /qm) - lnt 2 i + 1 - ti 2 ] + h - ( n/r 2 )I 2 
lny-f = In(r 2 /n) + [5In {q 2 n/qir 2 ) - lnx n + 1 - x 2l ] + I 2 - {r 2 /n)h 


the 0.0 and 1.0. values of x\. These are related to the parameters by various algebraic 
expressions, as shown in Table 6.2. 

The adjustable parameters are Aij for van Laar, A /y for Wilson, x /y and 042 for 
NRTL, and t lj for Uniquac. The r, q, and I parameters are pure component values 
which can be found in the book by Poling et al. (2000). 

An alternative approach is to recognize that a single data point yields two equilibrium 
equations with the parameters < 2 /y as unknowns. Solution of these equations yields good 
estimates of the a t j for a reasonable data point. 

Both of these approaches can be implemented easily using the Aspen Plus regression 
system. The former approach would involve using the entire data set for the regression 
with the initial values set to zero. If this produces good results, nothing more need be 
done; if not, the experimental activity coefficients are listed under regression/results/tab 
profiles and can be displayed with Aspen Plus’s plotting capability and extrapolated 
by eyeball. 

The latter approach can also be implemented with Aspen Plus by regressing a single 
data point. Subsequently, the regression with the full data set can be done with the 
estimated values of the pair as starting values. If one decides to regress the bij 
pair, an initial value can be obtained by noting that = b^ /T in the absence of a 
four-parameter fit. 

Aspen Plus provides a very useful Plotting Wizard which has many built-in plots 
related to regression of equilibrium data (see Figure 6.4, as well as the ability to 
select which variables to plot from a variety of Aspen Plus data sources. The Wizard 
is invoked by selecting Plot from the main display menu. All results shown may be 
obtained by selecting Properties/Regression/Rl with tab Profiles and may be plotted 
by selecting the column to be plotted and associating it with the desired axis available 
in the Plot menu. Figure 6.5 shows an example of a gamma(l) and gamma(2) versus 
x\ plot. 

As an example, Table 6.3 gives the T-P-x-y data for the system cyclohexane 
(CH)-isopropyl alcohol (IPA) of Nagata (1963). A plot of the activity coefficients of 
both components versus the mole fraction of component 1 (CH) is shown in Figure 6.5. 
Extrapolation of the curves produces infinite dilution activity coefficients of about 4.8 
for CH and 9.5 for IPA. 

Applying the infinite dilution form of the van Laar equation to the estimated infinite 
dilution coefficients yields 1.57 and 2.25 for a \2 and < 221 , respectively. Calculating b \2 
and Z? 2 i and noting that the original van Laar parameter A 12 is modified to 


A12 — ^12 + — 
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Figure 6.4 Aspen Plus plot wizard. 


TABLE 6.3 Vapor-Liquid Equilibrium 
Cyclohexane (1)-Isopropyl Alcohol (2) 


Temperature 

(°C) 

Pressure 

(mmHg) 

*1 

y\ 

80.29 

760 

0.021 

0.093 

78.1 

760 

0.047 

0.18 

78.85 

760 

0.091 

0.279 

73.8 

760 

0.147 

0.373 

72.13 

760 

0.217 

0.442 

70.88 

760 

0.279 

0.481 

70.13 

760 

0.373 

0.519 

69.99 

760 

0.478 

0.561 

69.56 

760 

0.602 

0.597 

69.79 

760 

0.748 

0.649 

70.99 

760 

0.882 

0.704 

74.61 

760 

0.972 

0.815 


then either A 12 = an or A 12 = bn/T when two parameters are used. Then for constant 
A 12 , an initial estimate for bn = a nT and similarly for b 2 \- For this example, 


b n = 1.57(273.1 + 80.3) = 530.1 
b 21 = 2.25(273.1 + 74.6) = 782.3 
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Figure 6.5 Activity coefficient plot. 


This is illustrated from Aspen Plus regression runs, which give the following results. 
When bij = 0: 

d\2 = 1.435 
<221 = 2.298 

When ciij = 0: 

b n = 498.8 
b 2 \ = 791.1 


Virtually identical results were obtained with least squares and maximum likelihood 
regressions. 

When setting up multiple regressions in one run, care must be taken to set the 
unused parameters equal to zero because Aspen Plus will use older or stored val¬ 
ues in the regression. The ideas illustrated above are appropriate for any choice of 
activity coefficient equation. Details of the regression runs above are available at 
Examples/chipavanLaar and Examples/chipwilson. 

6.5.2 Regression of LLE Data 

When fitting liquid-liquid equilibrium data the situation is somewhat different. Most 
extraction processes are isothermal; therefore, if binary mutual solubility data are avail¬ 
able, consistent with ternary (or higher) LLE data at the same temperature, the d l: j can 
be obtained by using one data point for regression. Since each data point can generate 
two equations, 


(Yi^i) 1 = (Yi*i) n 


(6.30) 


an exact solution is possible. 
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With suitable data it is possible to estimate the infinite dilution activity coefficients 
by assuming that when the composition of a component is near a mole fraction of 

I. 0, its activity coefficient is 1.0. For example, in the system A-B the composition of 
component A in phase I, x \, is 0.002 and the composition of component B in phase 

II, xf, is 0.01. Applying equation (6.30) for each component yields, approximately at 
infinity, y* = 495 and y° = 998. 

If the bi.j parameters are required because the application involves a variable tem¬ 
perature, several points must be regressed, although an initial value can be obtained 
from one data point using one of the methods above. 

When regressing ternary systems, two possibilities exist: (1) fit the two solvents 
with mutual solubility data and fit the four remaining parameters with ternary data 
while holding the two previously determined parameters constant; or (2) fit all six 
parameters with ternary data. Using the ternary data of Sugi and Katayama (1978) 
shown in Table 6.4, we look next at examples of both methods. 


TABLE 6.4 LLE Data for the System n -Hexane (2)-Ethanol (2)-Acetonitrile (3) 


Temperature °K 

Mole Fraction Phase I 

Mole Fraction Phase II 

A-Hexane 

Ethanol 

A-Hexane 

Ethanol 

313.15 

0.8831 

0.0166 

0.0968 

0.0879 

313.15 

0.8674 

0.0251 

0.1003 

0.1299 

313.15 

0.8546 

0.0332 

0.1054 

0.1622 

313.15 

0.8372 

0.0433 

0.125 

0.1942 

313.15 

0.8101 

0.0567 

0.1367 

0.2169 

313.15 

0.7821 

0.0772 

0.1522 

0.2331 

313.15 

0.7272 

0.1086 

0.2053 

0.2695 

313.15 

0.6912 

0.1269 

0.2249 

0.2748 


First Method An estimate of the mutual solubility data is made by extrapolating a 
plot of the data from Figure 6.6, resulting in x\ = 0.9 and xf = 0.09. The estimates 
calculated for the parameters are 


b l3 = -553.43076 
b 3l = -18.956075 

Holding b\ 3 and b 3 \ constant and fitting the remaining four parameters to the experi¬ 
mental data yields the following results: 


b n = 10.801713 
b 2 i = -169.74093 
b 23 = 215.31305 
b 32 = -178.56 3 57 

Second Method Parameters obtained by calculating the values of six parameters 
using one data point are 


b n = -494.10342 
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b n = 95.2601793 
bn = -530.14039 
b n = -35.602627 
Z? 2 3 = -1.3683628 
b 32 = -253.9741 


Parameters obtained using the full data set initialized with the values above are 

b n = 1892.22192 
b n = 136.837046 
b l3 = -461.75675 
b 3x = -68.070246 
b 23 = 12.3152652 
b 32 = 2285.7669 

All the methods described give passable results, and careful comparison of results 
with the composition range of the process would facilitate a choice of parameters. 
Aspen Plus’s Plot Wizard offers several possibilities for viewing experimental and 
calculated values of the variables. An example is shown in Figure 6.6. Details of the 
Aspen Plus regressions above are available at Examples/lle method one one point, 
Examples/lle method one two parameters fixed, Examples/lle method two one point, 
and Examples/lle method two. 

6.6 WORKSHOPS 

In preparation for the workshops it would be advisable to read the Aspen Plus documen¬ 
tation dealing with the activity coefficient models van Laar, Wilson, NRTL, Uniquac, 


U.&HJM.MIIllIflHIMlBiim’IT_ Jgjxj 


D-l Triangular plot of NIIFXANF ETHANOL ACN 



Figure 6.6 Triangular plot of LLE data of an n-hexane-ethanol-acetonitrile system. 
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Unifac, and Unifac-LL. Note the similarity between Uniquac and Unifac. Note the 
difference in application between Unifac and Unifac-LL. 

Workshop 6.1 Use Property Analysis to estimate the mutual solubility of a 
toluene-water system at temperatures varying from 32 to 86°F at 14.696 psi. Select 
the property option Unifac-LL. The mutual solubility of the toluene-water system 
at 86°F and 14.696 psi is 0.999894 and 0.000237 mole fraction water in each phase, 
respectively. Use engineering (ENG) units. Compare to the results estimated. 

Workshop 6.2a Fit the following single data point of mutual solubility for a 
toluene-water system to the two-parameter form of the van Laar equation. Use the 
parameters and At 25°C and 14.696 psi the solubility of water expressed as 
mole fractions in each phase is 0.999894 and 0.000237, respectively. 

Workshop 6.2b Fit the toluene-water solubility data from the LLE Dechema data 
collection in Table 6.5 to the van Laar equation using the ci l:J and cijj parameters. 
All solubility data are at 14.696 psi and expressed in mole fractions of toluene. Use 
engineering units. 

What are your observations? 


TABLE 6.5 Toluene-Water Mutual Solubility Data 


Temperature 

(°C) 

Mole Fraction 
x l 

Mole Fraction 
x 11 

0 

0.000142 

0.99891 

10 

0.000128 

0.99844 

20 

0.000113 

0.99784 

25 

0.000106 

0.99763 


Workshop 6.2c Repeat Workshop 6.2b using bij and bjj parameters with the van 
Laar equation. What are your observations? 

Workshop 6.3a Repeat Workshop 6.2b using bij and bjj parameters with the Uni¬ 
quac equation. What are your observations? 

Workshop 6.3b Using the data of Washburn et al. (1939) for the water-ethanol- 
toluene system in Table 6.6, fit the data to the Uniquac equation with the bij parameters 
determined from the single mutual solubility data point of Workshop 6.3a held constant. 
Vary the bij parameters for the remaining binary pairs. 

Workshop 6.3c Repeat Workshop 6.3b varying all six bij parameters? Compare 
the results to those of Workshop 6.3b. Which result should be used? 

Workshop 6.4a Fit the data of Suska et al. (1970) given in Table 6.7 to the van 
Laar, Wilson, NRTL, and Uniquac equations. Use the aj parameters for all equations. 
Additionally, use the cj parameter when fitting the NRTL equation. Select and fix the 
appropriate value of the cij parameter. Compare the four results. 
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TABLE 6.6 Liquid-Liquid Equilibrium Data for a Water-Ethanol-Toluene System 


Phase-W Mole Phase-T Mole 

Fraction Fraction 


Temperature (°C) Pressure (psi) Water Ethanol Water Ethanol 


25 

25 

25 

25 

25 

25 

25 

25 

25 

14.696 

14.696 

14.696 

14.696 

14.696 

14.696 

14.696 

14.696 

14.696 

0.95320 0.04621 

0.91360 0.08538 

0.87856 0.11990 

0.82021 0.17712 

0.77323 0.22240 

0.72620 0.26547 

0.69166 0.29572 

0.64420 0.33500 

0.59615 0.37190 

0.00356 0.00390 

0.00506 0.01385 

0.00703 0.02552 

0.01150 0.05210 

0.01960 0.08849 

0.02007 0.09905 

0.02425 0.11970 

0.03111 0.14779 

0.03719 0.16672 

TABLE 6.7 Vapor- 

-Liquid Equilibrium for a Tertiary Butanol (!)■ 

-Water (2) System 

Temperature (°C) 

Pressure (mmHg) 

Liquid Mole Fraction: 
Tertiary Butanol 

Vapor Mole Fraction: 
Tertiary Butanol 

86.50 

760 

0.0190 

0.4040 

83.80 

760 

0.0270 

0.4630 

82.00 

760 

0.0430 

0.4840 

81.70 

760 

0.0630 

0.4890 

81.40 

760 

0.1190 

0.4990 

81.10 

760 

0.1670 

0.5070 

81.00 

760 

0.2340 

0.5190 

80.90 

760 

0.3160 

0.5340 

80.65 

760 

0.3690 

0.5440 

80.50 

760 

0.3880 

0.5460 

80.30 

760 

0.4670 

0.5700 

80.00 

760 

0.5290 

0.5890 

80.00 

760 

0.6090 

0.6220 

80.00 

760 

0.6170 

0.6280 

80.00 

760 

0.6380 

0.6410 

80.10 

760 

0.6710 

0.6960 

80.20 

760 

0.7550 

0.7160 

80.30 

760 

0.8280 

0.7850 

80.90 

760 

0.8740 

0.8330 

81.25 

760 

0.9120 

0.8770 

81.80 

760 

0.9410 

0.9050 


Workshop Notes 

Workshop 6.1 In implementing this workshop an approach similar to one that would 
be used in a laboratory is advisable. For example, if equal quantities of toluene and 
water were placed in a beaker, mixed thoroughly, and permitted to stand for a sufficient 
time to equilibrate, the phases would form and could be sampled. Chemical analysis 
of the samples would yield the result desired. 
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Under the subject Analysis and Tab System, select “generate points along a flash 
curve” and valid phases “vapor-liquid-liquid” even though there will be no vapor. Set 
up a mixture for equilibration as described using the following procedure. Under the tab 
Variable, select pressure as a fixed state variable and temperature as the adjusted vari¬ 
able. Select and highlight temperature and the Range/List button will become available 
to enter temperature data. 

Workshop 6.2a The fit is exact because a single mutual solubility point is described 
by two equilibrium equations like equation (6.30) which will contain as unknowns 
only the two parameters of the activity coefficient equation. 

Workshop 6.2b and 6.2c As expected, use of the aij parameter produces values 
of solubility that are independent of temperature; however, even when using the bij 
parameters, the data are not well represented. 

Workshop 6.3a An exact solution is obtained as in Workshop 6.2a. 

Workshop 6.3b and 6.3c When viewing results using the Plot Wizard’s est vs. exp 
plot, there are small differences between the two approaches, but this is not a universal 
result. 

Use of the mutual solubility data point as in Workshop 6.3b establishes the location 
of the phase envelope termini (on a ternary diagram) on the toluene-water edge of 
the diagram and will force the composition at very dilute phases to approach the 
experimental mutual solubility. If the mutual solubility data are not consistent with the 
other experimental tie lines, the phase envelope in that area of the phase diagram will 
not appear smooth. If the data are fit as in Workshop 6.3c, the phase envelope will 
be smooth, but the mutual solubility predicted from the ternary data may not conform 
to the experimental mutual solubility. Because of this potential problem, an Aspen 
Plus user must be cognizant of the expected compositions of the simulation for which 
the regression is being prepared and make the appropriate choice when the data are 
regressed. 

Workshop 6.4a The fit of the data to the van Laar equation suggests that there is a 
heterogeneous region that lies between 0.1 to 0.25 mole fraction TbOH. A careful look 
at the experimental data shows that this may be so, but it is a judgment call. More 
experimental data may be needed. Extreme caution is advised. 

The Wilson equation fit by its nature cannot have a heterogeneous region and is 
otherwise satisfactory. The Uniquac fit exhibits behavior similar to that of van Laar. 
The NRTL fit with the c parameter = 0.2 exhibits similar behavior to that of van Laar. 

Workshop 6.4b The NRTL fit with the c parameter = 0.3 exhibits behavior similar 
to that of van Laar. 
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CHAPTER SEVEN 


FLASHES AND DECANTER 


Aspen Plus’s model library contains two different rigorous flash blocks that solve the 
appropriate material, energy balance, and equilibrium equations. The Flash2 block is 
designed to produce a single vapor phase and a single liquid phase that are in equi¬ 
librium when the flash conditions are specified such that a multicomponent mixture’s 
state is in the two-phase region. Similarly, the Flash3 block is designed to produce one 
vapor phase and two liquid phases in equilibrium for suitably specified process condi¬ 
tions. The Flash3 block is also capable of solving a liquid-liquid equilibrium problem 
under conditions where no vapor is produced. A similar block is Decanter, which is 
designed to produce two liquid phases in equilibrium in the absence of a vapor phase. 
The Flash2, Flash3, and Decanter blocks can be found in the model library under the 
Separators tab. 

7.1 Flash2 BLOCK 

Figure 7.1 is a graphical depiction of a Flash2 block. The nomenclature used is as 
follows: F represents the total feed flow, in moles/time; fi the flow of component 
i in the feed, in moles/time; V the total vapor flow, in moles/time; v L the flow of 
component i in the vapor, in moles/time; L the total liquid flow, in moles/time; /; the 
flow of component i in the liquid, in moles/time; and n the number of components. 
Including the flash temperature Tf and pressure Pf results in 2n + 2 independent 
variables given the feed state. These are equilibrium temperature; flash pressure; two 
total flows, L and V ; and 2 (n — 1) component flows; or alternatively, 2 n component 
flows, excluding the total flows. Mole fractions are calculated from the independent 


Teach Yourself the Basics of Aspen Plus™ By Ralph Schefflan 
Copyright © 2011 John Wiley & Sons, Inc. 


93 





94 FLASHES AND DECANTER 



variables by an equation such as 



where y, is the mole fraction of component i in the vapor. 

The applicable material balances are n componential equations such as 

fi ~ Vi -h= 0 (7.2) 

or alternatively, n — 1 equations such as equation (7.2) and one overall material balance 
given by 


F -V - L = 0 (7.3) 

Additionally, n equilibrium equations which describe the equality of the fugacities 
of components in each phase are required. When the liquid fugacity is represented by 
an equation of state, where 4>f and $Y are the fugacity coefficients of component i in 
the liquid and vapor phases, respectively, the result is 

-Xi$i = 0 (7.4a) 

When the liquid fugacity is represented by an activity coefficient equation, where 
y i is the activity coefficient of component i and the vapor phase is represented by an 
equation of state 


yitfP-YiXipf=0 (7.4b) 

where p? is the vapor pressure of component /, results. For the sake of simplicity, the 
Poynting correction (see Prausnitz et al., 1999), which has a contribution only for very 
light components, has been omitted. 
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An overall energy balance, where h with a suitable subscript represents the enthalpy 
of the feed, liquid, and vapor, respectively, and Q represents the heat added, given by 

h F F + Q-h L L- h v V = 0 (7.5) 


completes the set of 2n + 1 equations. 

The Flash2 block permits the specification of two of the four possibilities: flash tem¬ 
perature, flash pressure, heat required, and the fraction of the feed vaporized, except the 
combination Q and V/F. Various permutations of the specifications require other mod¬ 
ifications of the equations and the list of unknowns. For example, if V/F is specified, 
equation (7.5) is replaced by 


I) 


specified 


v j 

EU fj 


= o 


(7.6) 


If 2, the heat added or removed, and the flash pressure are specified, the 2 n com¬ 
ponent flows and the flash temperature can be calculated. If the flash temperature and 
pressure are specified, the sequence of the solution of the equations changes. Equations 
(7.1) through (7.4) are solved simultaneously, after which equation (7.5) is solved 
for Q. 

An example of the use of a Flash2 block is given at Chapter Seven Examples/ 
Flash2Example. An important aspect of the setup, which is applicable to all simula¬ 
tions, is the selection of the data bank from which the activity coefficient or equation 
of state parameters were derived. Figure 7.2 shows that the binary interaction Wilson 
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Figure 7.2 Flash2 example Wilson parameters. 
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QSpecifications Flash Options Entrainment Utility 


Flash specifications 


T emperature 

^ | 70 

r 

^3 

Vapor fraction 
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“3] 
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Figure 7.3 Flash2 example specifications. 
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Figure 7.4 Flash2 example results. 


parameters were obtained from the data source apv70vle-lit. The upper and lower tem¬ 
peratures of applicability are also shown. If this set is deemed inappropriate, selection 
of the tab Databanks permits the user to select an alternative source of data. The input 
specifications for this example are shown in Figure 7.3, and the flash results are given 
in Figure 7.4. 

A useful feature of the Flash2 and some other equilibrium stage blocks is the ability 
to specify entrainment; that is, a specified fraction of the liquid phase is carried by 
the flashing vapor into the overhead. A comparison of the performance of a Flash2 
block with and without entrainment is given at Examples/Flash2EntrainmentExample. 
The flowsheet is shown in Figure 7.5. A comparison of the overheads produced by 
the flash is given in Figure 7.6. The composition of stream 22 is lower in the more 
volatile component, methanol, and its fraction vapor is less than 1.0. 

7.2 Flash3 BLOCK 

Figure 7.7 is a graphical depiction of a Flash3 block. The nomenclature used is as 
follows: F represents the total feed flow, in moles/time; ft the flow of component 
i in the feed, in moles/time; V the total vapor flow, in moles/time; v t the flow of 
component i in the vapor, in moles/time; L 1 the total flow of the first liquid phase, in 
moles/time; //, the flow of component i in the liquid, in moles/time; L 2 the total flow 
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Figure 7.5 Flash2 entrainment example. 
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Figure 7.6 Flash2 entrainment results. 
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of the second liquid phase, in moles/time; and if the flow of component i in the liquid, 
in moles/time. Including the flash temperature Tf and pressure Pf results in 3n + 2 
independent variables, where n is the number of components given the feed state. The 
variables consist of the equilibrium temperature, the flash pressure, and the three total 
flows (i.e., and L 1 , L 2 , and V ), and 3 (n — 1) componential flows or, alternatively, 3 n 
componential flows, excluding the total flows. Mole fractions are calculated from the 
independent variables by an equation such as equation (7.1). 

The applicable material balances are n componential equations, such as 

fi - »/ - l\ -if = o (7.7) 

or alternatively, n — 1 equations such as equation (7.7) and one overall material bal¬ 
ance, given by 


F -V -L l -L 2 = 0 (7.8) 

Additionally, 3 n equilibrium equations which describe the equality of the fugacities 
of components in each phase can be written, but only 2 n are independent. The vapor- 
phase fugacity is represented by an equation of state where the are the fugacity 
coefficients of component i in the vapor phase, and the liquid phases are represented 
by an activity coefficient equation where y- and y 2 are the activity coefficients of 
component i and p? is the vapor pressure of component i . Then any two of 


Yi x i — yf x f = 0 

(7.9a) 

yi$P — y! ’ x l pf = o 

(7.9b) 

yi tfP-Y 2 ,xfpf = 0 

(7.9c) 


apply. For the sake of simplicity the Poynting correction (see Poling et al., 2000), which 
has a contribution only for very light components, has been omitted from equations 
(7.9b) and (7.9c). The overall energy balance given by equation (7.5) completes a set 
of 3n T 1 equations: 
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h F F + Q-huL 1 -h L2 L 2 -h v V = 0 (7.10) 

The Flash3 block permits the specification of two of the four possibilities: flash tem¬ 
perature, flash pressure, heat required, and the fraction of the feed vaporized. Various 
permutations of the specifications require other modifications to the equations and the 
list of unknowns. For example, if Q, the heat added/removed and the flash pressure 
are specified, the 3 n component flows and the flash temperature can be calculated. 

When specifying the block it is important to identify the key component in the 
bottom takeoff (liquid 2) such that the block conforms with its physical behavior (i.e., 
the heaviest of the solvents should be selected). An example of the use of a Flash3 
block is given at Examples/Flash3Example. 

7.3 DECANTER BLOCK 

Figure 7.8 is a graphical depiction of a Decanter block. The nomenclature used is as 
follows: F represents the total feed flow, in moles/time; ft the flow of component i 
in the feed, in moles/time; L 1 the total flow of the first liquid phase, in moles/time; ij 
the flow of component i in the liquid, in moles/time; L 2 the total flow of the second 
liquid phase, in moles/time; and if the flow of component i in the liquid, in moles/time. 
Including the decanter temperature T c i and pressure Pd results in 2n + 2 independent 
variables, where n is the number of components given the feed state. The variables 
consist of the equilibrium temperature and the decanter pressure, the two total flows 
L 1 and L 2 , and 2{n— 1) componential flows or, alternatively, 2 n componential flows, 
excluding the total flows. Mole fractions are calculated from the independent variables 
by an equation such as (7.1). 

The applicable material balances are n componential equations, such as 

fi - 1} - if =0 (7.11) 



Figure 7.8 Decant model. 
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or alternatively, n — 1 equations such as equation (7.11) and one overall material bal¬ 
ance, given by 


F — L l — L 2 =0 (7.12) 

Additionally, n equilibrium equations can be written which describe the equality of the 
fugacities of components in each phase, where y \ and y\ are the activity coefficients 
of component i and xj and xf are the mole fractions of component i in phases 1 and 
2, respectively: 


y}x}-y}xf = 0 


(7.13) 


An overall energy balance, 


hfF + Q — hi\L l — hi 2 L 2 = 0 


(7.14) 


completes a set of 2n + 1 equations. 

Most decantation processes are carried out isothermally, and the heat associated 
with solute moving between phases is usually negligible, thus eliminating the need 
for equation (7.14) and reducing the number of unknowns to 2 n when the pressure 
is specified. This is reflected in the block’s input forms, where the only allowable 
specifications are temperature and pressure. As in the Flash3 block, care must be 
taken to specify the key component in the second liquid product, which should be 
the heavier phase. The feed from Flash3Example was used to create an example of 
the use of a Decanter block, which may be found in Examples/DecanterExample. It 
would be useful to compare results of Flash3 and Decanter since both are carried out 
at the same conditions: 1 psia and 62°F. The Flash3 specification is V/F = 0.1, which 
produces a temperature of 61.9°F, which is close enough for a comparison. Figure 7.9 
shows the data sources for the Uniquac parameters, which are LLE data compared 
to VLE data for Flash3. This is an example of a conundrum. Which data produce 
the correct answer? Both seem reasonable. Aspen Plus has solved the appropriate 
equations, but it is up to the user to select the correct result. Perhaps some laboratory 
work will be required. 

The block also permits an alternative formulation of the basic equations as a Gibbs 
free energy minimization problem. This formulation and some methods of solution are 
described by Walas (1985). Additional details are provided in Chapter Ten. 

An interesting aspect of the Decanter block is the possibility of its use to simulate a 
batch extraction process. The describing equations are identical to those of a continuous 
decanter if the flows per unit time are replaced with a batch charge as a feed and the 
batch removal of the two products. A common practice in industry is to place a solvent 
A containing a quantity of C, the product desired, into a vessel. Following this, the 
solvent B is charged to the vessel. After agitation the products are allowed to settle. This 
is the equivalent of feeding both streams as a mixture to the decanter, producing the 
products L 1 and L 2 . Common industrial practice is to remove the equilibrium B stream, 
now containing extracted C, and recharging the vessel with fresh B. This process is 
easily simulated with several Decanter blocks in series, as shown in Figure 7.10. 
The combined feeds are repeatedly settled, the extract and raffinate removed, and the 
resulting product sent to the next stage. Each stage represents the same vessel one 
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Figure 7.9 Alternative data sources. 



Figure 7.10 Batch extraction process. 


time step later in the process, when it contains the remaining raffinate and is charged 
with fresh solvent C. The composition of the mixed extracts is the product of the 
mixer. 

7.4 WORKSHOPS 

Workshop 7.1 Use the Flash2 block to flash the feeds in Table 7.1 adiabatically at 
constant pressure for the purpose of reducing the acetone composition in the bottoms 
product. The feed temperature is 70°C and its pressure is 760 mmHg. 

Be sure to set up appropriate units. Use Aspen Plus’s stored Uniquac, NRTL, and 
Wilson parameters as the basis for the calculations. What do you observe? 
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TABLE 7.1 

Feeds for Workshop 7.1 

Component 

Feed 1 (lbmol/hr) 

Feed 2 (lbmol/hr) 

Ethanol 

40 


Water 


50 

Acetone 

10 



Workshop 7.2 Based on the solution to Workshop 7.1, using 760 mmHg as a base 
case, create a sensitivity study for pressure varying between 10 and 760 mmHg. Show 
the effect of varying flash pressure on fraction vaporized, mole fraction acetone in the 
vapor and liquid, and flash temperature using Aspen Plus’s stored Uniquac parameters. 

Workshop 7.3 The feeds in Table 7.2 are to be fed to a Flash3 block at a pressure 
of 760 mmHg and a temperature of 80°C. Using Aspen’s stored Uniquac parameters, 
use a sensitivity study to find the temperature at which the vapor phase disappears 
when operating at 760 mmHg. The temperature and pressure of both feeds are 25°C 
and 760 mmHg. 


TABLE 7.2 

Feeds for Workshop 7.3 

Component 

Feed 1 (lbmol/hr) 

Feed 2 (lbmol/hr) 

Ethanol 

10 


Water 


40 

Toluene 

50 



Workshop 7.4a Using the feeds from Workshop 7.3, employ the Decant block 
to calculate the composition and quantity of the resulting phases at 25 °C using the 
parameters provided in Table 7.3. Be sure to set up appropriate units. 


TABLE 7.3 Uniquac Parameters for a Toluene-Water-Ethanol System 


Parameter/ 

Element No. 

Component Pair 

Parameter Estimate 
(SI units) 

UNIQ/2 

Toluene-H 2 0 

-926.6623 

UNIQ/2 

H 2 0-Toluene 

-258.8106 

UNIQ/2 

EtOH-H 2 0 

301.8589 

UNIQ/2 

H 2 0-EtOH 

-119.5998 

UNIQ/2 

EtOH-Toluene 

-49.37053 

UNIQ/2 

Toluene-EtOH 

-100.2406 


Workshop 7.4b Develop a sensitivity study to determine the effect of temperature 
on this operation. Vary the temperature from 30 to 100°C, displaying the amount of 
ethanol extracted into the water-rich phase and the ethanol composition of the toluene- 
rich phase. Be sure that the base case is executed last. Note all observations concerning 
this sensitivity study. Are the results reasonable? 
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Workshop Notes 

Workshop 7.1 A duplicator block will permit all three simulations to be done at once, 
and the results will facilitate easy comparison. Two of the results are comparable. The 
third displays a similar trend, but which are correct? 

• Do the activity coefficient models work for this system? 

• Were the parameters in the database fit properly? 

• Should new data be located and custom-fit? 

• Would an “in the field” design be acceptable for controlling the product compo¬ 
sition (e.g., a secondary sensor loop feeding the set point of the primary control 
loop, to make up for model deficiencies)? 

Workshop 7.3 A sensitivity study with a plot should be a good guide for a design 
specification. See the solutions to Workshops 7.4a and 7.4b. 

Workshop 7.4 All Uniquac parameters except Uniq/2 equal must be set to zero. 

• Where did the list of Uniquac parameters come from? How can their quality be 
assessed? 

• Results show that extraction improves as the temperature gets colder. 
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CHAPTER EIGHT 


PRESSURE CHANGERS 


Aspen Plus offers six different types of pressure changers. The Pump, Compressor, 
and Multistage Compressor models actively change pressure as a result of the addition 
of some form of work to the block which is converted to action on a stream. The 
three other models—valve, pipe (single segment with fittings), and pipeline (multiple 
segments with fittings)—are stationary objects that incorporate resistance to flow which 
must be overcome by use of a pump or compressor. 

8.1 PUMP BLOCK 

The Pump block typically is used to calculate the power required to pump a fluid up to 
a specified pressure. The Pump block can be used to model a hydraulic turbine as well. 
In this variation the block calculates the power produced when the discharge pressure 
is specified. Alternatively, the power input can be specified and the pump will calculate 
the discharge pressure. Figure 8.1 shows the Pump and Turbine input options. Details 
of the calculations are well documented in Aspen Plus’s Help under Pump Reference. 
An example of the use of a pump and a hydraulic turbine may be found at Chapter 
Eight Examples/pumps. For cases in which meaningful pumping applications are not 
needed, a Heater block can be used to specify the state of a stream. 

8.2 COMPR BLOCK 

Details of single-stage compressor modeling can be found in Aspen Plus’s Help under 
the Compr Reference, with a variety of links to all variants. An excellent overview 
may be found by following the link Mollier method, Compressor Background. These 
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^/Specifications | Calculation Options | Flash Options | Utility | 
Model 

(* Pump r Turbine 


Pump outlet specification 


(* Discharge pressure: 

|100 

| psia 

^3 

r Pressure increase: 

1 

| psi 

“3 

r Pressure ratio: 

1 



r Power required: 

1 

|hp 

3 


r Use performance curve to determine discharge conditions 


Efficiencies 

Pump: [l Driver: [T 


Figure 8.1 Pump and turbine input. 


sections should be read prior to using the block. The basic idea is given by 


p Ah 



( 8 . 1 ) 


where p is either a polytropic or isentropic efficiency factor which is assumed constant 
for a process, where Ah is the change in enthalpy between the two states, V the 
volume, and p the pressure. The various methods available in Aspen Plus are based 
on assumptions dealing with how the integral is evaluated. 

The polytropic efficiency is given by 


^ = (8 .2) 

n rip 

where n is the polytropic coefficient, with each state described by pV n , k is the heat 
capacity ratio c p /c v , and r\ p is the poly tropic efficiency. 

For an isentropic process h out is obtained from the outlet pressure when the entropy 
at the outlet equals the entropy at the inlet. The isentropic efficiency is used to modify 
an isentropic process to approach a real condition. There are two forms of isentropic 
efficiency: 


for compression and 




h s — h • 

^out "in 
^out ^in 


^out ^in 

~h s — h■ 

"out "m 


(8.3) 


(8.4) 


for expansion. 
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Figure 8.2 Compr options. 


It is also necessary to specify the mechanical efficiency, defined by 

IHP = FAh (8.5a) 

IHP 

BHP =- (8.5b) 

hm 

where F is the molar flow rate, Ah the enthalpy change, IHP the indicated horsepower, 
BHP the brake horsepower, and \} m the mechanical efficiency. 

The Aspen Plus Compr block can be used to model a compressor or turbine. Three 
models of compressors are available: a polytropic centrifugal, a polytropic positive 
displacement, and an isentropic compressor, each with variants as shown in Figure 8.2. 
If a turbine is chosen, isentropic is the only possibility. Some examples of Compr 
applications are given at Examples/compressor. 

8.3 MCompr BLOCK 

Details of multistage compressor modeling can be found in Aspen Plus’s Help under 
the MCompr Reference, and the link specifying MCompr. These sections should be 
read prior to using the block. The MCompr block can be used to model six types of 
compressors, as shown in Figure 8.3. The tabs Specs and Cooler are used to specify the 
various efficiencies and the stage outlet conditions, respectively. The Convergence tab 
is used to select the allowed phases in the compressor and cooler stages. An example 
of a MCompr application is given at Examples/mcompr. 

8.4 PIPELINES AND FITTINGS 

The Pipe block calculates the pressure drop in a single run of pipe. Pipe permits the 
connection of fittings. The flow is assumed to be uniform along the axis and can 




























108 PRESSURE CHANGERS 


^/Configuration | y Material | Heat-Work | ySpecsj y Cooler] y Convergence | 


Number of stages: 



Compressor model 

(* Isentropic C Polytropic using ASME method 

C kentropic using ASME method C Polytropic using GPS A method 
C Isentropic using GPSA method C Positive displacement 





■j|_iciii >_■ □uuri i. y |_ic 

(* Fin discharge pressure from last stage: [200 | psia ^ | 

C Fix discharge conditions from each stage 

C Use performance curves to determine discharge conditions 



ndurm upuun 

I - Use performance curves to determine shaft speed 




(• Calculate at suction conditions C Calculate at average conditions 


u | | - 


ncm uaptduuy utdiuuitduuri 

Calculation method: Rigorous ^ | 




Figure 8.3 MCompr options. 


^Operation | ^Valve Parameters | Calculation Options | y Pipe Fittings | 


Calculation type 

C Adiabatic flash for specified outlet pressure (pressure changer) 

C Calculate valve flow coefficient for specified outlet pressure (design) 
(* [Calculate outlet pressure for specified valve (rating! 


-Pressure specification 
(* Outlet pressure: 

C Pressure drop: 


20 

psia 


psi 


Valve operating specification 
(* % Opening: 

C Flow coef: 


50 


Flash options 




| Vapor-Liquid 

t | Maximum iterations: |30 —H 


Error tolerance: 


y Operation | y Valve Parameters | Calculation Options ypipe Fittings | 


W [Include effect of pipe fittings head loss on valve flow capacity 
Valve and pipe diameters 
Valve inlet diameter: 


Inlet pipe diameter: 
Outlet pipe diameter: 


3 

in ^ | 


ft jd 

1- 

F—3 


Figure 8.4 Pipeline and fittings. 
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perform one, two, or three phase calculations. If there are multiple runs of pipe of 
different diameters or elevations, the Pipeline block should be used. Figure 8.4 shows 
the Pipe input options. An example of usage is given at Examples/pipeline. 

8.5 WORKSHOPS 

Workshop 8.1 The stream whose composition is given in Table 8.1 is to be pumped 
from 25 to 100 psi. The inlet temperature is — 10°C. Use the Pump block to calculate 
the horsepower required if the pump’s efficiency is 80%. 


TABLE 8.1 Feed for All Chapter 8 Workshops 


Component 

Abbreviation 

Lbmol/Hr 

Methane 

Cl 

0.1 

Ethane 

C2 

1 

Propane 

C3 

10 

n -Butane 

NC4 

18.9 

1-Butane 

1C4 

20 

1,3-Butadiene 

DC4 

20 


Workshop 8.2 The feed of Workshop 8.1 is a vapor at 100 psi and 100°C. It is to 
be compressed to 500 psi using a polytropic compressor with an efficiency of 100%. 
The efficiency of the motor driving the compressor is also 100%. Use the Compr block 
to calculate the performance of the compressor and the outlet state of the product. 

Workshop 8.3 Repeat Workshop 8.2 using an isentropic compressor that operates 
reversibly (i.e., its isentropic efficiency is 100%). The efficiency of the motor driving 
the compressor is also 100%. 

Workshop 8.3a Repeat Workshop 8.3 with an isentropic efficiency of 0.8. 

Workshop 8.4 Repeat Workshop 8.2 using a two-stage isentropic compressor where 
both stages and motors have 100% efficiencies. No interstage coolers are to be used. 
All pressure drops are zero. Use the block MCompr to calculate the compressor’s per¬ 
formance and the state of the outlet stream. Compare results to those of Workshop 8.1. 

Workshop 8.5 Using Workshop 8.4 as a base case, run case studies to determine the 
effect of an intercooler between stages 1 and 2 of the compressor. For the case studies, 
specify the heat removed as 0, 50,000, 100,000, 175,000, and 200,000 Btu/hr. Show 
the effect of heat removal on total brake horsepower, total heat removal, temperature 
of the gas at the outlet of stage 1 of the compressor, and temperature of the gas at the 
outlet of the cooler. 

Workshop 8.5a Repeat the 175,000-Btu/hr heat removal case from Workshop 8.5 
but include a liquid removal stream. 
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Workshop Notes 

Workshop 8.2 and 8.3 These two workshops permit comparison of poly tropic and 
isentropic compressors. 

Workshop 8.3 and 8.3a These two workshops permit comparison of ideal and real 
isentropic compressors. 

Workshop 8.4 This is a straightforward two-stage MCompr example. 

Workshop 8.5 and 8.5a Heat removal requires a minus sign. Note the warning 
accompanying the — 17,500-Btu/hr removal case. The effect of an intercooler is shown. 
The warning disappears in Workshop 8.5a when liquid removal is provided for. 


REFERENCE 

Aspen Plus version 7.0, Help: Pump Reference, Compr Reference, and MCompr Reference. 


CHAPTER NINE 


HEAT EXCHANGERS 


The implementation of heat exchanger models is different from that of most Aspen 
Plus models in that some are capable of detail design using very high quality industrial 
programs which are integrated into Aspen Plus. These are the Hetran and Tasc shell- 
and-tube heat exchanger and the Aerotran air-cooled exchanger programs. They are 
documented in Aspen Plus Help: EDR (Exchanger Design and Rating). Most can be 
used for design as well as being used to rate existing exchangers; that is, they can be 
used as sequential modular models. 

To simulate a heat exchanger one must solve the primary equations 


_ hot /rp hot T^hot\ n 

q ~ m hot c p (T in - T oui ) = 0 

_ ™ ^cold/^cold T^coldx 

q m co\d C p (Tout Un ) — ^ 

(j — UAF ATlm — 0 

ATi - AT 2 
LM lnCAri/A^) 


(9.1) 

(9.2) 

(9.3) 

(9.4) 


Here q is the exchanger duty; m a flow rate, c p the heat capacity; T the temperature; 
AT the temperature difference at an end of the exchanger; U an overall heat transfer 
coefficient which depends on temperature, transport properties, and exchanger geom¬ 
etry; and F a correction factor for multiple tube-side and/or shell-side passes. The 
factor F, derived through the work of Nagle (1933) and Underwood (1934), can be 
calculated as 


Teach Yourself the Basics of Aspen Plus™ By Ralph Schefflan 
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Vfl 2 + 1 ln[(l - S)/(l - RS)] 


(9.5) 


(R - 1) In 2 5( * +1 


2-S(^+l+VK 2 +l) 


where 


r r hot 'r'hot 

in 1 out 


(9.6) 


r r cold tp cold 

1 out 1 in 


r r cold nr cold 

1 out 1 in 


(9.7) 


'rhot 'T'cold 

1 in 1 in 


When used in simulation mode, the state of the exchanger feeds must be specified. 
Depending on the complexity of the model chosen, the heat transfer area, A, is either 
specified or calculated from the heat exchanger physical layout. Heat transfer coeffi¬ 
cients are calculated from correlations such as the Hewitt (1992) correlation prepared 
by Gnielinski, which involves the Nusselt Anu = hD/k, Reynolds Are = DG/ |x, and 
Prandtl Apr = c p \i/k numbers, and the Darcy friction factor, given by 



(9.8) 


Id = (l-821og 10 TVre — 1.64) -2 


(9.9) 


Here h t is the inside pipe heat transfer coefficient, D t the inside pipe diameter, k 
the thermal conductivity, G the mass flow rate, |Ji the viscosity, c p the heat capacity, 
and f D the Darcy friction factor. Complete documentation of all correlations used in 
Aspen Plus can be found in Help: Heatx Reference and Model Reference. Depending 
on which model is chosen, U is either specified or calculated iteratively during the 
convergence process. The four heat-exchange-related models can be found in the model 
library under the tab Heat Exchangers. 

9.1 HEATER BLOCK 

An example of the primary input form of the Heater block, which shows the possible 
specifications, is shown in Figure 9.1. The Heater block offers a variety of ways to 
specify the output stream state, all of which result in calculation of the energy required 
to heat (or cool) a stream. Alternatively, one may specify the energy added to or 
removed from a heater, which is used by the block to establish the state of the output 
stream. 

An important capability is the use of two heaters to model a heat exchanger 
bypassing the use of equations (9.3) and (9.4), as shown in Figure 9.2. Note the use 
of a heat stream to connect the two heaters. The heat stream should be aligned in the 
correct direction, which depends on which heater will receive the heat, either positive 
or negative. Care must be taken with the sign of the heat transferred (heat added is 
positive). In this example the outlet temperature of heater H2 is specified and the heat, 
stream 5, flows to heater HI. This example may be found at Chapter Nine 
Examples/Heaters. 
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© Specifications | Flash Options | Utility | 


Flash specifications 

Heat duty ^ 1 

| Btu/hr 

Temperature 

Vapor fraction 

Heat duty 

Temperature change 

Degrees of superheating 

Degrees of subcoolinq 

fpsia 



Pressure-drop correlation parameter 



Figure 9.1 Heater block primary input. 




Figure 9.2 Two heaters modeling a heat exchanger. 


An example of a process in which heat exchangers are replaced with heaters is shown 
in Figure 9.3. The heat exchangers here are used to conserve energy by extracting heat 
from the bottoms of the columns and transferring part of the available energy to the 
column feeds. An important aspect of using heaters rather than heat exchangers is that 
the calculations can be performed once through, and no tear streams are required. This 
is illustrated by Aspen Plus’s analysis of the process, given in Figure 9.4. In most 
cases, whenever a heat exchanger is placed on a flowsheet, a tear stream is introduced. 
This is illustrated in Figure 9.5, where the heaters from the Figure 9.3 process are 
replaced by heat exchangers. Aspen Plus’s analysis, in Figure 9.6, shows that one tear 
stream has been generated for each heat exchanger inserted. Details of the solution to 
the process with four heaters are given at Examples/fourheatersflowsheet. Similarly, 
the two-heat-exchanger process is provided at Examples/twoheatxflowsheets. 
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Figure 9.3 Process with two heaters as heat exchangers. 


COMPUTATION ORDER FOR THE FLOWSHEET: 
B3 B1 BG B2 B6 B4 

-Calculations begin 


Black: 

B3 

Model: HEATER 

Black: 

El 

Model: SEP2 

Block: 

B5 

Model: HEATER 

Black : 

B2 

Model: SEP2 

Black: 

B5 

Model: HEATER 

Black: 

B4 

Model: HEATER 

- Circulation 

calculations completed 


No Errcra or Warnings Generated * 


Figure 9.4 Computation sequence for a process with heaters. 
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Flowsheet Analysis : 

Block $OLUER01 { Method:: WEGSTEIN) has been defined to converge 
streams: 3 

Blcck $0LUERQ2 4Methed: WEGSTEIN) has been defined tc converge 
streajr.s : 4 

COMPUTATION ORDER FOR THE FLOWSHEET: 

$OLUER01 B3 B1 
{RETURN $OLVER01) 

$OL¥ERQ2 B4 B2 
^ RETURN $OLUER02 3 


Figure 9.6 Tear streams for a process with heat exchangers. 


9.2 HEATX BLOCK 

The Heatx block has very wide applicability, as shown in its primary input form shown 
in Figure 9.7. To use a Heatx exchanger in a simulation, it must have been designed so 
that it can be used as part of the sequential modular approach. Note that Heatx permits 
a design, a rating, and a simulation mode. The simulation mode requires that the inputs 
to the exchanger be either specified or calculated from a previous block, because the 
block is expected to behave as a sequential modular model. The rating mode permits 
the specification of a block output as shown in Figure 9.8 and behaves as if it were a 
sequential modular model with a design specification imposed on an output. 

The simplest approach for using the block for design is to select a shortcut design, 
in which case an output of an exchanger is specified. The available specifications for 
design are shown in Figure 9.7. The solution to the design problem, using the shortcut 
method, yields a value for the UA product which can be employed as part of the 
specifications when setting up a simulation problem. An example of this approach 
is given in the solution to the flowsheet shown in Figure 9.5, which incorporates 
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© Specifications | LMTD | Pressure Drop ] Methods | Film Coefficients | Utilities 


Calculation 

(* Shortcut C Detailed 
r ShellUube 
C AirCooled 
r Plate 
r Hetran 
r Aerotran 

r tasc toshortcut 


Flow arrangement 

Hot fluid: |~ ” | 

Flow direction: itercurrent ▼ | 

No. shells in series: pi 
|” Calculate number of shells 


Type: —ZI 

Exchanger Type: | Heat exchanger ▼ | 


Exchanger specification 

Specification: 

Z! 

Value: 

Exchanger area: 

Constant UA: 

Minimum tempera 

Hot stream outlet temperature 

Hot stream outlet temperature decrease 

Hot outlet-cold inlet temperature difference 

Hot stream outlet degrees subcooling 

Hot stream outlet vapor fraction 

Hot inlet-cold outlet temperature difference 

Cold stream outlet temperature 



Cold stream outlet temperature increase 

Cold stream outlet degrees superheat 

Cold stream outlet vapor fraction 

Exchanger duty 

Hot / cold outlet temperature approach 



Figure 9.7 Heatx block primary input form. 


exchanger design. The specification used for design, the outlet temperature of the cold 
stream, which was used for both exchangers, was transferred from the four-heater 
flowsheet to the two-exchanger flowsheet. The UA values calculated for the shortcut 
design of two exchangers were transferred to the simulation flowsheet to produce 
the completed process simulation. Details may be found at Examples by comparing 
twoheatxflowsheetd (design) and twoheatxflowsheets (simulation). When using rating 
or simulation calculations, both the shortcut and detailed options permit the user to 
enter estimates of the overall heat transfer coefficients from sources such as Perry and 
Green (1999). In addition to an overall heat transfer coefficient, the heat transfer area 
must be provided. 

The detailed option permits only rating and simulation modes, both of which can be 
used as part of a sequential modular simulation. In both cases the selection of geometry 
prompts a user for the detailed layout of a shell-and-tube heat exchanger. A detailed 
design must be completed prior to the selection of this option. 

As can be seen from the calculation frame in Figure 9.7, Aspen Plus provides a 
rich selection of detailed design calculation methods. These are referred to as Aspen 
Exchanger Design and Rating. An example of the use of “shell and tube” is given by 
redesigning the first exchanger at Examples/twoheatxflowsheetd. The input is virtually 
identical to that at Examples/twoheatxflowsheetd except for the identification of the 
hot fluid and the specification of an input file, as shown in Figure 9.9. This example 
is given at Examples/B3shellandtube. The file B3shellandtube.txt gives a report of the 





























HEATX BLOCK 117 



LMTD | Pressure Drop | Methods | Film Coefficients 


Calculation 
(* Shortcut 
r ShellScT ube 
C AirCooled 
r Plate 
C Hetran 
C Aerotran 
r TASC 


C Detailed 


T ransfer UA 
to shortcut 


Flow arrangement 
Hot fluid: 

Flow direction: 

No. shells in series: 

I” Calculate number of shells 


| Countercurrent ^ 



Type: 


| Rating 


C3 


Exchanger Type: |Heat exchanger “3 


Exchanger specification 

Specification: 

| Cold stream outlet temperature ▼ | 


Value: 

Hot stream outlet temperature 


Exchanger area: 

Constant UA: 

Minimum tempera 

Hot stream outlet temperature decrease 

Hot outlet-cold inlet temperature difference 

Hot stream outlet degrees subcooling 

Hot stream outlet vapor fraction 

Hot inlet-cold outlet temperature difference 

Cold stream outlet temperature 



Cold stream outlet temperature increase 

Cold stream outlet degrees superheat 

Cold stream outlet vapor fraction 

Exchanger duty 

Hot / cold outlet temperature approach 



Figure 9.8 Heatx block’s rating specifications. 


^Input File | Property Curves | Analysis Parameters 


Shell-Tube Geometry 


EDR Input File 

Hetran/Aerotran input file with full path and file extension ( K .bjt): 

TASC input file with full path and file extension ( K .tai): 
Shell&Tube/AirCooled/Plate input file with full path and file extension ( K .edr): 

|B3SHELLANDTUBE.EDR Browse... I 


F Save hot stream and curve properties 
F Save cold stream and curve properties 
Save analysis parameters 
W Save final geometry from design mode 
1^ Save all to Shell&Tube /AirCooled /Plate input 


T ransfer geometry to 
Hetran 


T ransfer geometry to 
ShellScT ube 


Utilities 


Figure 9.9 Mheatx primary input form. 







































118 HEAT EXCHANGERS 


1=) Exchanger Geometry ▼ | S 

£3 Shell&Tube 
$ Cl Input 
Results 

S Cl Input Summary 
ffl-Pl Result Summary 
El Cl Thermal / Hydraulic Summary 
El ^ Mechanical Summary 

1^1 Exchanger Geometry 
(=) Setting Plan & Tubesheet Layout 
|=| Cost l Weights 
S Cl Calculation Details 


US 3 M | <=> | « » | N-> | # □ ► ▼ | 

Basic Geometry | Tubes | Baffles | Supports-Misc. Baffles| Bundle| Enhancements! 


Unit Configuration 


Exchanger Type 


BEM 

T ube number 



48 

Position 


Hor 

T ube length actual 


ft 

12 

Arrangement 


1 par 1 ser 

T ube passes 



2 

Baffle type 


Single segmental 

T ube type 



Plain 

Baffle number 


54 

Tube O.D. 


in 

0.75 

Spacing (center-center) 

in 

2.5 

T ube pitch 


in 

0.9375 

Spacing at inlet 

in 

4.25 

T ube pattern 



30 



Shell 

Kettle 

Front head 


Rear head 

Outside diameter 

in 

8.G25 


8.896 


8.696 

Inside Diameter 

in 

8.071 


8.071 


8.071 




Shell Side 

Tube Side 

Nozzle type 


Inlet 

Outlet 


Inlet 

Outlet 

Number of nozzles 


1 

1 


1 

1 

Actual outside diameter 

in 

1.315 

1.315 


2.375 

1.9 

Inside diameter 

in 

1.049 

1.049 


2.067 

1.61 

Height under nozzle 

H in| 

0.552 

0.5226 




Dome inside diameter 

in 






Vapor belt inside diameter 

in 






Vapor belt inside width 

H in| 






Vapor belt slot area 

in 2 






Impingement protection 


No 

impinqement 

No 

impinqement 

No 

impinqement 



Distance to tubesheet 

in 

140 

3.75 





Figure 9.10 Mheatx zone output. 


performance of the exchanger; however, mechanical and other details are available 
only after executing the B3shellandtube.bkp file and using the EDR browser in the 
report mode of the results. An example is given in Figure 9.10. Note that there are 
many detailed types of outputs which are selected from the left-hand-side column under 
“shell and tube” by clicking on the plus sign, which exposes the reports available. The 
choice of a different detailed design program involves the same techniques as those 
described above. 

9.3 MHEATX BLOCK 

The Mheatx block Is unique in that it is able to handle multiple hot and cold streams. 
The block is reorganized internally into a collection of Heater blocks connected by heat 
streams. An example is given in Examples/mheatxexample. Here there is one Mheatx 
block with two hot and two cold feeds. The primary input form is shown in Figure 9.11. 
Note that a specification is required on all but one output stream. Selection of the tab 
Zone Analysis permits the assignment of zones within the heat exchanger. Results are 
shown in Figure 9.12. The analysis summary gives a breakdown of the important heat 
exchange factors, including identification of the locations where streams should enter 
and leave the exchanger. 

9.4 WORKSHOPS 

Workshop 9.1 Use the shortcut feature of the Heatx block to calculate the state of 
the cold stream and the heat transfer area required when the streams in Table 9.1 are 
passed through a countercurrent heat exchanger with an overall heat transfer coefficient 
of 100 Btu/(hr-ft 2 -°R). The outlet temperature desired for the hot stream is 70°F. The 
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^Specifications] Zone Analysis j Flash Table | Options | 


Specifications 



Inlet stream: 

1 

3 

5 

7 


Exchanger side: 

HOT 

COLD 

HOT 

COLD 


Outlet stream: 

2 

4 

6 

8 


Decant stream: 






Valid phases: 

Vapor-Liquid 

Vapor-Liquid 

Vapor-Liquid 

Vapor-Liquid 

► 

Specification: 

Temperature 

T emperature 

Temperature ^ | 



Value: 

100 

110 

T emperature 

Vapor fraction 

T emperature change 

Deg Super 



Units: 

F 

F 



Pressure: 





Units: 

psia 

psia 

Heat duty 

Molar heat duty 

Mass heat duty 

Fraction of heat duty 

psia 


Duty estimate: 





Units: 

Btu/hr 

Eltu/hr 

Eltu/hr 


Max. iterations: 

30 

30 

30 

30 


T olerance: 

0.0001 

0.0001 

0.0001 

0.0001 


Figure 9.11 Shell-and-tube EDR input file specification. 


Streams | Balance | Exchanger Zone Profiles | 


Profiles 

View: | Summary 



Point 

Cumulative 

Duty 

T-Hot 

T -Cold 

Delta-T 

LMTD 

Zone duty 

Zone UA 

Cumulative 

UA 

Pinch point 

Stream/ 

Phases 

| Btu/hr ▼ | 

k jlI 

l F jd 

If jd 

|f jlJ 

| Btu/hr 

|Btu/hr-F ▼ | 

iBtu/hr-f ▼ | 

1 d 

1 l 


1 

0 

100 

70 

30.0000002 








2 

33126.0809 

104.730108 

80 

24.7301085 

27.2802724 

33126.0609 

1214.28629 

1214.28629 


IN 7 


3 

92295.4629 

113.113353 

90.1095895 

23.0037635 

23.8565265 

59169.402 

2480.21865 

3694.50494 




4 

184590.926 

126.048069 

105.681058 

20.3670106 

21.6586435 

92295.4629 

4261.36858 

7955.87352 




5 

210448.487 

129.637541 

110 

19.6375407 

20.0000585 

25857.5616 

1292.87429 

9248.74781 


OUT 3 


8 

213064.851 

130 

111.016887 

18.9831131 

19.3084785 

2616.36399 

135.503374 

9384.25119 


IN 5 


7 

276886.389 

140 

135.489822 

4.51017819 

10.0701382 

63821.5374 

6337.70224 

15721.9534 

GBL 



Figure 9.12 EDR browser in report mode. 


TABLE 9.1 Feed for Workshop 9.1 


Component 

Hot Feed (lb/hr) 

Cold Feed (lb/hr) 

Methanol 

200 


Water 

1800 


Ethylene glycol 


400 
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temperatures of the hot and cold feeds are 140 and 40°F, respectively. Both streams 
are at 14.696 psi. Use the Wilson property method. 

Workshop 9.2 Use two Heater blocks with a heat stream transferring the heat 
between the blocks for the specifications given in Workshop 9.1. 

Workshop 9.3 Create an input file duplicating the input of Workshop 9.1. Change 
the problem type to simulation. Remove the exchanger specification, that is, the outlet 
temperature of the hot stream. Insert the UA value calculated in Workshop 9.1 and 
insert the overall heat transfer coefficient of 100 Btu/(hr-ft 2 -°R) in the U Methods 
form. Run and compare the design to the simulation. 

Workshop 9.4 Repeat Workshop 9.3 using a cocurrent heat exchanger. Compare the 
results to Workshop 9.3. How do the results compare? What is the primary source of 
the differences? 

Workshop 9.5 An existing heat exchanger is to be used to cool 145,000 lb/hr of 
benzene at 390°F and 400 psia. The coolant is 490,000 lb/hr n-dodecane at 100°F 
and 200 psia. Aspen Plus’s Heatx block is to be used with the detailed option. Search 
Aspen Plus’s Help for exchanger configuration for general information on exchanger 
geometry. 

Details of the exchanger design are as follows: 

• Tema type E: one shell pass, two tube passes 

• Countercurrent flow 

• Horizontal alignment 

• Hot fluid in the shell 

• 2.75-ft inside shell diameter 

• 0.5-inch shell-to-bundle clearance 

• No sealing strips 

• U calculated from film coefficients 

• No fouling 

• 200 tubes, length 32 ft, pitch 1.25 inches, square layout, carbon steel, inside 
diameter 0.875 inch; outside diameter 1 inch 

• 24 segmental baffles; baffle cut 0.2 

• Tubesheet-to-baffle spacing and baffle-to-baffle spacing 1.2 inches 

• Tubes in baffle window 

• Shell-side nozzles 6 inches 

• Tube-side nozzles 8 inches 

Use Aspen Plus’s defaults when necessary. 

Workshop 9.6 Repeat Workshop 9.5 using the Rating Option. 

Workshop 9.7 Repeat Workshop 9.5 using the Simulation Option. 
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Workshop Notes 

Workshop 9.1 This is a design problem using shortcut methods. After providing an 
exchanger specification (i.e., hot stream outlet temperature and an overall heat transfer 
coefficient, both defined in the workshop description), the design is straightforward. 

Workshop 9.2 The results of Workshop 9.1 are duplicated. 

Workshop 9.3 As this is a simulation workshop it is necessary to change the spec¬ 
ification of Workshop 9.1 to simulation and to insert the value of UA calculated by 
Workshop 9.1. The hot stream specification must also be removed. The simulation 
duplicates the results of Workshop 9.1. 

Workshop 9.4 Use of the exchanger design from Workshop 9.1 but with cocurrent 
flows results in poorer performance. The heat transferred in Workshop 9.3 is about 
131,000 Btu/hr, compared to 99,000 Btu/hr in Workshop 9.4. The log mean temperature 
difference for the two cases is 35.1 and 26.6°F, respectively, which explains the poorer 
performance of Workshop 9.4. 

Workshop 9.5 The solution is straightforward. 

Workshop 9.6 The hot stream specification is taken from the solution of Workshop 
9.5. The results are virtually identical to those of Workshop 9.5. 


REFERENCES 

Aspen Plus version 7.0, Heaters documentation. 

Aspen Plus version 7.0, Help: EDR, Exchanger Configuration, Heatx Reference, and Model Reference. 
Hewitt, G. F., Ed., Handbook of Heat Exchanger Design, Begell House, New York, 1992. 

Nagle, W. M., Mean temperature differences in multipass heat exchangers, Ind. Eng. Chem., 25, 604-609 
(1933). 

Perry, R. H., and Green, D. W., Perry’s Chemical Engineers’ Handbook on CD- ROM, 7th ed., McGraw-Hill, 
New York, 1999, Table 11.2. 

Underwood, A. J. V., The calculation of the mean temperature differences in multipass heat exchangers, 
J. Inst. Pet. Technol. , 20, 145-158 (1934). 


CHAPTER TEN 


REACTORS 


Aspen Plus provides seven different models of reactor blocks that deal with chemical 
reactors in different ways: 

1. RStoic: a stoichiometry-based reactor with specified extents of reaction. 

2. RYield: a reactor based on specified yields. 

3. REquil: a rigorous equilibrium reactor based on reaction stoichiometry. 

4. RGibbs: a rigorous reactor which includes phase equilibrium using Gibbs free 
energy minimization. 

5. RCSTR: a rigorous stirred-tank reactor with rate-controlled reactions based on 
specified kinetics. 

6. RPlug: a rigorous plug flow reactor with rate-controlled reactions based on spec¬ 
ified kinetics. 

7. RBatch: a rigorous batch and semibatch reactor with rate-controlled reactions 
based on specified kinetics. 

The reactor blocks can be found in the model library under the tab Reactors. 

10.1 RStoic BLOCK 

The RStoic block is the simplest of Aspen Plus’s reactor blocks. It permits the use 
of several reactions with the molar extent of conversion or fractional conversion of a 
component, specified for each reaction. For the reaction 

aA T —> cC T dY) (10.1) 
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the molar extent of reaction X is defined by 


A n a A nt> An c Arid 

a bed 


( 10 . 2 ) 


where A n t are the moles of component i either created or consumed. The applicable 
material balance equations with the feeds and products of component /, defined as F L 
and Pi , respectively, are for a reactant. For example, if i is component A, 


Pi = F t - aX 


(10.3) 


and for a reaction product, for example, if i is component D, 

Pi = Fi+dX (10.4) 


This is illustrated at Chapter Ten Examples/RStoicExample and is shown in Figure 10.1. 
Here a duplicator block feeds each RStoic reactor block. The reaction involved is the 
formation of ammonia from hydrogen and nitrogen. The reactor ST1 is configured for 
an extent of reaction of 90 lbmol/hr and reactor ST2 with a fraction converted of 0.50 
nitrogen. The details of the extent of reaction configuration is shown in Figure 10.2. 



Figure 10.1 RStoic example. 



Figure 10.2 Configuration of extent of reaction. 
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Material | Vol.% Curves | Wt. % Curves | Retro. Curves | Poly. Curves | 


Display: | Streams Format: | FULL -3. Stream Table | 









[1 

F 

-1 

~w 


Substream: MIXED 







Mole Flow lbmol/hr 







H2 

400.0000 

265.0000 

250.0000 


d 


N2 

100.0000 

55.00000 

50.00000 



NH3 

0.0 

30.00000 

100.0000 



T otal Flow lbmol/hr 

500.0000 

410.0000 

400.0000 



T otal Flow Ib/hr 

3607.700 

3607.700 

3607.700 



T otal Flow cuft/hr 

3524.412 

2302.163 

2631.437 



Temperature F 

332.0000 

332.0000 

332.0000 



Pressure psia 

2204.332 

2204.332 

2204.332 



Figure 10.3 Comparison of extent and yield results. 


Note that the shaded fraction conversion is 0.45 of N 2 , which is equivalent to 90 
lbmol/hr of conversion. Aspen Plus does not provide this value. It was placed here for 
the convenience of the reader. The results of the simulation are given in Figure 10.3, 
where it may be observed that the fraction converted, 0.5 in ST2, results in 5 additional 
lbmol/hr of N 2 converted relative to ST1. A value of 0.45 would produce results 
identical to those of ST1. 

10.2 RYield BLOCK 

The RYield block is used when the distribution of products is known. No reaction 
stoichiometry is involved. Two primary options are available: specification of compo¬ 
nent yields and specification of component mapping. The other two options, a user 
subroutine and petro specifications, are not considered here. 

For the specification of component yields, the yield of each reaction product is spec¬ 
ified as either moles or mass per unit mass of feed. A component may be specified as an 
inert. The yields will then be based on a unit mass of noninert feed. Calculated yields 
are normalized to maintain an overall material balance. Because of this, yield specifica¬ 
tions are actually expressed as ratios of components in the product of the reactor. Yield 
specifications result in a yield distribution rather than absolute yields. An example of 
the use of the component yield option may be found at Examples/RYieldExamplel. 
This example involves the reaction of benzene and chlorine to form mono- and diben¬ 
zyl chloride with hydrochloric acid as a by-product. In the absence of kinetic data 
the ratio of products is specified as yields. The yields for this example are shown in 
Figure 10.4. It should be noted that the multiplication of the yield data by a factor of 10 
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^Specifications ^Yield | Flash Options j PSD | Comp. Attr. 
Yield specification 

Yield options: | Component yields ^ | 

Component yields 



Component 

Basis Yield 


CHLORINE 

Mass 

0 


HCL 

Mass 

0.001 


C6H6 

Mass 

0 


C6H5CL 

Mass 

0.04 


C6H4CL2 

Mass 

0.04 

* 





Figure 10.4 Component relative yield specifications. 

produces the same results as Example 1, which demonstrates that the yield distribution 
is independent of the numerical values of the yields, but merely on the component 
ratios. This is easily verified by temporarily changing the yield data in RYieldExam- 
plel, executing, and comparing results to those in the file RyieldExamplel.txt, which 
is based on the original specifications. 

The same set of reactions is used to demonstrate the component mapping option 
of the RYield block. No reaction stoichiometry is employed. A reaction can be char¬ 
acterized as lumped, in which case a list of reactants are lumped together to form 
a specified product. Each of the reactants is assigned a coefficient which defines the 
weight fraction of the reactant converted to the product. 

When a reaction is characterized as delumped, the reactant selected produces a 
specified list of product components, and similar to lumping, each product is assigned 
a coefficient which represents the weight fraction of the reactant that is converted to 
each product. An example of the use of the component yield option may be found 
at Examples/RYieldExample2. Figure 10.5 shows the implementation of component 
mapping for the benzene chlorination example. Selecting a reactant and choosing Edit 
permits selection of the lumping or delumping coefficients. 

In both block options selection of inconsistent yield factors or component mapping 
factors may produce material balance discrepancies. In this case, Aspen Plus normalizes 
the results to maintain the material balance while maintaining the component ratios 
calculated. 

10.3 REquil BLOCK 

The REquil block is used when the reaction stoichiometry is known and some of 
the reactions reach chemical equilibrium. Additionally, the block imposes phase 
















REquil BLOCK 127 


y Specifications | ^Yield j Flash Options | PSD | Comp. Attr. ^/Comp. Mapping | Utility | 


Lumping / delumping reaction 



Reaction 

TyP e 


► 

CGH5CL 

LUMP 

0.5 CHLORINE + 0.5CGHG --> CGH5CL 


CGH4CL2 

LUMP 

0.5 CGH5CL + 0.25 CHLORINE --> CGH4CL2 


HCL 

LUMP 

0.25 CHLORINE -> HCL 


Figure 10.5 Lumping and component mapping. 


equilibrium on the products of reaction. Two primary options are available for each 
reaction: specification of the extent of reaction as described by equations (10.1) 
through (10.4), or calculation of the equilibrium constant for the reaction K &q from 
the standard state change Gibbs free energy, A G°, of the components of the reaction, 
using the equation 


K ea=e -(A g °/RT) 


(10.5) 


where 


a o o o 

^ ^products ^ ^reactants 


( 10 . 6 ) 


The effect of temperature on the reaction equilibrium constant can be calculated as 


In 


K t 


t 2 


eq 

k t 1 


4 / 


A h° 
~T^ 


dT 


(10.7) 


Subsequently, the compositions of the products and reactants at equilibrium are 
calculated using equation (10.7): 


K q q = 


fl /product 
n /reactants 


( 10 . 8 ) 


where / is the fugacity of a component. 

An example of REquil, the chlorination of benzene, using the approach to the equi¬ 
librium option may be found at Examples/REquilExamplel. The reaction setup and 
specifications are shown in Figure 10.6. Note that the selection “temperature approach” 
is 0°F, which indicates equilibrium. An example of REquil using the “extent of reac¬ 
tion” option, with the same reactions, may be found at Examples/REquilExample2. 
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</Specifications "/Reactions Convergence | Entrainment | Utility | 


Reactions 




Rxn No. 

Specification type 

Stoichiometry 



1 

T emp. approach 

CGH6+ CL2 <-•> C6H5CL + HCL 



2 

T emp. approach 

CSH5CL + CL2 <-> P-CGH4CL + HCL 



^dEdit Stoichiometry 


Reaction No.: 
Reactants 





Component 

Coefficient 

Solid 

*1 


CGH6 

-1 

No 


CL2 

-1 

No 

zi 

* 






Component 

Coefficient 

Solid 



CGH5CL 


No 

H 


HCL 


No 


* 




zi 


Products generation 


C Molar extent: 

r=r 

3] 

(• T emperature approach: 11 

F“ 

jd 

Extent estimate: |3 

| Ibmol/hr 

d 


Reaction number. 


Figure 10.6 Equilibrium reaction stoichiometry. 


10.4 RGibbs BLOCK 

The RGibbs block can be used to establish the chemical equilibrium composition 
between reactants and products. Consider the reaction of isomer A to isomer B at a 
given temperature and pressure. For an n -component mixture where X[ is the mole 
fraction of component i , the Gibbs free energy g m | x is given by 


girtix = £ Xi + RT £ Xi In Xi (10.9) 

Differentiating equation (10.9) with respect to x a and recognizing that for a two- 
component mixture x b = 1 — x a , one obtains 

dgmix o o . -w~% m -I /1 r\ i r\ \ 

- 7 — =g a -g b + RT In —— (10.10) 

dx a t x a 

The minimum of g m j x can be found by setting equation (10.10) to zero and solving 
for x a . The same approach can be taken for complex chemical and phase equilibria. 
This method is especially useful for liquid-liquid equilibrium, as it does not produce 
spurious solutions as may be the case for models such as Decant and Flash3, which 
solve the describing equations with the appropriate activity coefficient equations. 

An example of the use of the RGibbs block for the chlorination of benzene is given at 
Examples/RGibbsExamplel. Figure 10.7 shows the basic setup and calculation options 
available. Figure 10.8 shows detailed specifications and product options. It is important 
to note that no reactions need to be defined. One may optionally choose to eliminate 
from consideration components from the mixture. Examples/RgibbsExample2 is used 
to illustrate the possibility of excluding a possible product. In this case the component 
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^Specifications! Products | Assign Streams | Inerts | Restricted Equilibrium | Utility | 


Operating conditions 


Pressure: |l 4.7 

| psia 


Temperature: |80 

F“ 

T 

Heat duty: | 

| Btu/hr 

"ZJ 


Calculation options 


Calculate phase equilibrium and chemical equilibrium ▼ 


| Calculate phase equilibrium only 


9 Calculate phase equilibrium and chemical equilibrium 

Restrict chemical equilibrium - specify temperature approach or reactions 
Restrict chemical equilibrium - specify duty and temp, calc temp approach 


W Include vapor phase Solid Phases 


Figure 10.7 RGibbs basic specifications. 



Figure 10.8 RGibbs product options. 


chlorobenzene has been removed, as shown in Figure 10.9, and indeed, a look at the 
solution shows that the chlorination benzene to chlorobenzene does not take place. For 
this system of reactions, the block offers the possibility to consider additional products: 
for example, the inclusion of other dichloro isomers in the list of components. 

10.5 REACTIONS FOR THE RIGOROUS MODELS 

Prior to using the rigorous models, it is necessary to completely define the reac¬ 
tions to be considered. The Reactions setup may be found at the Data Browser 
display. When the Reactions folder is opened, two possibilities are presented, Chem¬ 
istry for electrolytes, and Reactions for all others. When reaction(s) are configured 
they are assigned an identifier, such as R-l, for use in a reactor. As an example, 
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^ Specifications ^Products] Assign Streams | Inerts j Restricted Equilibrium | Utility | 

C' RGibbs considers all components as products 
(* Identify possible products 
C Define phases in which products appear 


Products 



Component 

Valid phases 


HCL 

Mixed 


C6H6 

Mixed 


PC6H4CL2 

Mixed 


C6H3CL3 

Mixed 

* 




Figure 10.9 RGibbs product identification. 


Figure 10.10 shows the configuration of one of the reactions that occurs during the chlo¬ 
rination of benzene. Note that reactant coefficients require a minus sign and products 
a plus sign. Additionally, a check box is available to define the reaction as reversible. 
Four reaction class options are displayed: Equilibrium, Powerlaw, LLHW [i.e., the 
Langmuir-Hinshelwood-Hougen-Watson kinetics model (Perry and Green, 1999)], 
and GLLHW (a generalized LLHW model). Various tabs become available, such as 
equilibrium based on the reaction model selected. It is possible for different reaction 
classes to exist simultaneously within a reactor. 

In all cases, when any rigorous reactor block is used, a Reactions tab appears. This 
will permit the selection of a reaction group for the current reactor. 

10.5.1 Equilibrium Class 

When the Equilibrium class is chosen, the Equilibrium tab becomes available and 
Figure 10.11 is displayed. The choice of equilibrium constant by Gibbs energies or by 
a built-in empirical expression is available, as is the choice of the phase in which the 
reaction takes place. 

10.5.2 Powerlaw Class 

When the Powerlaw class is chosen, the Kinetic tab becomes available and Figure 10.12 
is displayed. It is important to note that there are six different choices for the concen¬ 
tration of a component. For example, when molarity is chosen from the drop-down list, 
the appropriate unit, kmol/cubic meter, is displayed but is hidden temporarily under¬ 
neath the drop-down list. The Powerlaw kinetic expression is composed of two parts, 
the kinetic factor and the driving force. For all cases the kinetic factor is given by 

r =k(—Y e (~E/R)a/T-l/T 0 ) 

\ToJ 


( 10 . 11 ) 
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^Configuration] ^Kinetic] Equilibrium | Activity | GLHHWAdsorption | 



New... 


Figure 10.10 Powerlaw kinetic options. 


where k is the reaction rate constant, T and Tq are the reaction and reference tem¬ 
peratures, E , is the activation energy, and n, is an arbitrary coefficient. The input 
for the kinetic factor can be configured without a reference temperature by setting n 
and To equal to zero. Setting the activation energy to zero leaves only the reaction 
rate constant. The driving force for the reaction depends on the concentration basis 
selected. For example, when the Driving Force button is pushed and molarity is chosen 
as the concentration basis, the driving force is the continuous product of the compo¬ 
nent concentrations raised to exponents defined in Figure 10.13. The rate expression 
becomes 


r =kU Cf (10.12) 

Searching Aspen Plus’s Help for kinetic factors will yield an entry Rate Con¬ 
trolled Reactions, which provides much detailed documentation—most important, the 
power-law expressions used for each possible concentration basis. 
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yConfiguration j Kineti ^Equilibrium | Activity | 



Figure 10.11 Equilibrium constant options. 


yConfiguration ^Kineticj Equilibrium | Activity ] 


|l) CL2+ BENZENE-> CC6H5CL + 

HCL 



Reaction class: | 

Reacting phase: 

Liquid 

“31 

Powerlaw kinetic expression 

[Ci] basis: 

Molarity 

“31 


r = [Kinetic factor][Driving force] 


[Ci] units: 
Rate basis: 


Kinetic factor 


0.001 


w 

F 

| kcal/mol ▼ | 

— 

p— 


Rate units: 


Molarity 
Molality 
Mole fraction 
Mass fraction 
Partial pressure 
Mass concentration 



Driving Force 


Adsorption 


Figure 10.12 Reacting phase options. 
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11) CL2 + BENZENE ~> CC6H5CL + HCL Reaction is reversible 


Reacting phase: 


Liquid 


[Ci] basis: 



- Reverse rate calculation option- 

<* Compute reverse rate using microscopic reversibility 
C Specify rate const, and cone, exponents for forward and reverse reactions 

Concentration exponents 

Specify parameters for | Forward ■*"[ reaction rate 



Reactants Comp. 

Exponent 


CL2 

i 


BENZENE 

1 



Products comp. 

Exponent 


CC6H5CL 



HCL 



- Coefficients for driving force constant 


i 

B: l 

i 

C: l 

i 

D:| 

r 


Close 


Figure 10.13 Reaction concentration coefficients. 


An example of the use of the Reactions entry on the Data Browser with the 
RCSTR block may be found at Examples/RcstrExamplel. This example involves a 
simplification of the chlorination of benzene to chlorobenzene. The reaction setup for 
this example is shown in Figures 10.12 and 10.13. 

10.5.3 Langmuir-Hinshelwood-Hougen-Watson Class 

The LHHW model (Perry and Green, 1999) involves reactions that include adsorption 
terms in the reaction rate expression 

(kinetic factor) (driving force) 

rate =- (10.13) 

adsorption term 

The kinetic and driving force terms are identical to those in the Powerlaw class. The 
adsorption expression is given by 

adsorption expression = [E Ki(Yl C”)]J 


(10.14) 
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</Configuration Kinetic | Equilibrium | Activity J GLHHW Adsorption 


) CL2 + BENZENE --> CC6H5CL + HCL 


Reaction class: | LHH\ 

LHHW kinetic expression 

[Kinetic factorJ[D riving force] 


Kinetic factor 
k: 


|0.001 

I 5 ” 


l°"~ 

1 kcal/mol ▼] 

1 

|F Jd 


[Ci] basis: 
[Ci] units: 
Rate basis: 
Rate units: 


| Liquid 


(Molarity 

L 

|kmol/cum 

|Reac (vol) 

3] 

| kmol/cum-s 


Solids 

J 

Driving Force | 

Adsorption | 


^Adsorption Expression 


Ads.=[Terml +Term2+...pAds. Expression exponent; Termi = Ki”[C1 pnl *[C2] A n2. 

Adsorption parameters 

Reacting phase: | Liquid 

[Ci] basis: | Molarity 

Adsorption expression exponent: | 

Concentration exponents 



Component 

T erm no. 1 

* 




Adsorption constants 
Term No.: f 

A|— 


23 


b r 


□ r 


Figure 10.14 LHHW kinetic and adsorption expressions. 


where i and j are component indices, v is a concentration exponent, and m is an 
adsorption expression. The K { are adsorption equilibrium constants, given by 

In Ki = A; + y- + Ci In 7 + D t T (10.15) 

where T is temperature in kelvin and A, B, C, and D are user-provided parameters. 

An example of the reaction setup is given in Figure 10.14. When the adsorption 
button is pushed, the right-hand display appears and facilitates the entry of adsorption 
data. 

10.5.4 Generalized-Langmuir-Hinshelwood-Hougen-Watson Class 

The GLHHW class is designed to enable the user to prepare a customized adsorp¬ 
tion term, which is best illustrated by Figure 10.15, which is displayed when the 
GLHHW tab is selected. Note the equation near the bottom of the figure which 
begins with Ads.=. This is the complete denominator: that is, the adsorption term 
in the rate equation. The equation is composed of terms each containing a product 
of Ki [equation (10.13)] and exponentiated component concentration terms. The exact 
details of the resulting combination can be seen from the artificial example shown in 
Figure 10.15. Note that the complete series of products is raised to the power of the 
adsorption exponent. 

10.6 RCSTR BLOCK 

The basic ideas describing the analysis of steady-state continuously stirred reactors 
is well described by Levenspiel (1962). A key element of the analysis is that the 
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</Configuration |y Kinetic] Equilibrium | Activity ^GLHHW Adsorption| 
GLHHW adsorption parameters 

GLHHW No: |yi t | Adsorption exponent: [05 


Concentration exponents 



Component 

T erm no. 1 

T erm no. 2 

Term no. 3 

Term no. 4 


CL2 

1 

2 




HCL 

1 

0 

3 


* 







Ads.={K1 K [CL2] K [HCL] + K2 K [CL2]^2 + K3 K [HCLp3r0.5 



Figure 10.15 GLLHW adsorption parameters. 

composition of the reactor is at all times equal to the composition of the reactor 
product. In essence, for each component entering and leaving a CSTR, the product is 
calculated by the material balance by accounting for disappearance of the reactant, and 
accumulation of a product, by the rate of reaction multiplied by reactor volume adjusted 
using appropriate stoichiometric coefficients. Similarly, an overall energy balance is 
employed which takes into consideration the energy change due to the reaction. 

The primary RCSTR input shown in Figure 10.16 gives all the reactor options 
available. The minimum additional data required represent the association of the reactor 
with the reactions involved by using the Reactions entry on the Data Browser. The 
complete example setup and solution are available at Examples/RCSTRExamplel. 

10.7 RPlug BLOCK 

Unlike the RCSTR block, the composition of an RPlug reactor changes along its length. 
The basic describing equations involve the integration of appropriate composition and 
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^Specifications! Streams | <f Reactions | PSD | Component Attr. | Utility | Catalyst | 


Operating conditions 


Pressure 


1 

|atm 

| T emperature 

jd 

|55 

l c -d 


Holdup 
Valid phases: 
Specification type: [*R& 


| Liquid-Only 


U 2nd Liquid 

- 


Reactor 

Volume: |0.5 

Resi. time: 


Reactor volume 


Residence time 

Reactor volume & Phase volume 
Reactor volume & Phase volume fraction 
Reactor volume & Phase residence time 
Residence time & Phase volume fraction 
Phase residence time & volume fraction 

n 


time: 




|cuft_; 




|hr lJ 


Figure 10.16 Primary RCSTR Input. 


rate terms along the reactor length. An overall energy balance is also employed. Details 
may be found in Levenspiel (1962). 

The primary specifications for an RPlug block are shown in Figure 10.17. Each of 
the reactor types available has an auxiliary set of inputs. For example, in Figure 10.18, 
a reactor with a constant coolant temperature, the overall heat transfer coefficient 


^Specifications | ^ Configuration | Strea | Reactions | *f Pressure | ioldup ] Catalyst | 


Reactor type: Reactor with specified temperature 


Reactor with specified temperature 


Operating coi Adiabatic reactor 


r 

r 


Constant 
Constant 
T emperat 


* 


Location 


Reactor with constant coolant temperature 
Reactor with co-current coolant 
Reactor with counter-current coolant 
Reactor with specified coolant temperature profile 
Reactor with specified external heat flux profile 


Temperature 

[T^ 


Figure 10.17 RPlug options. 
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& Specifications | ^ Configuration | Streams | <f Reactions | *f Pressure | Holdup j Catal 

Reactor type: | Reactor with constant coolant temperature —3 

Operating condition 
Heat transfer specification 


(* Specify heat transfer parameters 
U (coolant-process stream): 

C Calculate in user subroutine 


|Watt/sqm-K ^ | 



Figure 10.18 Rplug reactor type and heat transfer specifications. 


^/Specifications ^/Configuration | Streams | ^Reactions | ^Pressure | Holdup j Catalyst j Diameter | 



I” Multitube reactor 

F Diameter varies along the length of the reactor 



Figure 10.19 RPlug size and layout. 


must be provided or, alternatively, calculated in a subroutine. Additionally, the cooler 
temperature must be provided. When the tab Configuration is selected, Figure 10.19 
is presented. Here data dealing with the physical layout of the reactor, as well as 
the reaction phase are input. Selection of the Pressure tab permits specification of the 
process stream inlet pressure and various pressure-drop options. The complete example 
setup and solution appear at Examples/RPlugExamplel. 
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10.8 RBatch BLOCK 

The RBatch block is designed to evaluate the performance of a batch or semibatch 
reactor. The basic describing equations involve the equality of the rate of accumu¬ 
lation of a product or reactant within the reactor to the rate of gain or loss due to 
chemical reaction. Solution is by integration as a function of time. An overall energy 
balance as a function of time is also employed. Details may be found in Levenspiel 
(1962). 

The initial setup for Rbatch is shown in Figure 10.20, in which the reactor oper¬ 
ating and pressure specifications may be selected. Each of the specifications brings 
up unique entry requirements, as shown in Figure 10.21. It is also important to note 
that a specification for valid phases is required. If vapor-liquid-liquid is chosen, the 
second liquid button provides a means of identifying the second phase, as shown in 
Figure 10.22. 

Selection of the tab Stop Criteria permits a wide variety of possibilities for the 
selection of a means of stopping the reaction as shown in Figure 10.23. More than 
one stop criterion may be selected. Various drop-down lists, such as variable type, 
are available as needed. Figure 10.24 shows the options dealing with batch cycles 
and control of the number and interval of the results of the simulation. For this 
example, Figure 10.25 shows mole fractions of each component as a function of 
time. Additional results, as a function of time, may be found under the tabs Over¬ 
all and Feed. The complete example setup and solution are available at Examples/ 
RBatchExample 1. 


^Specifications | ^Reactions | V Stop Criteria | ^Operation Times | Continuous Feeds | Controllers | 


Reactor operating specification 


Constant temperature 


Constant temperature 
T emperature profile 
Constant heat duty 
Heat duty profile 
Constant coolant temperature 


Heat transfer user subroutine 


Pressure specification 
|Specify reactor pressure 



Valid phases 


Reactor: 

"73 

Vent accumulator: Vapor-Only 



2nd Liquid 


Figure 10.20 RBatch initial setup. 
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^Specifications | ^Reactions | Stop Criteria | ^Operation Times | Continuous Feeds | Controllers | 


Reactor operating specification 
| Constant coolant temperature ^ | 


Temperature: | 


u: r 

Watt/sqm-l ▼ | 

Area: |~ 

|sqm ▼ | 


Valid phases 
Reactor: 


Pressure specification 


| Calculate reactor pressure 
Reactor volume: 


"3 


|cum 

"3 

Vent opening pressure: 

|N/sqm 

3 ] 

Vent closing pressure: 


|N/sqm 

3 

Catalyst loading: 

[kg 

31 


| Liquid-Only 

T 


| Vapor-Only 

3 

2nd Liquid | 


Figure 10.21 RBatch operating and pressure specifications. 


Q Specifications | y Reactions | y Stop Criteria | y Operation Times | Cont 


Reactor operating specification 
Constant coolant temperature 


Temperature: 

U: |- 

Area: | 


|Watt/sqm-l ▼ | 
[sqm ^ 



Pressure specification 
| Calculate reactor pressure 


|cum ▼ | 

Vent opening pressure: 


t closing pressure: 

—peer 


Catalyst loading: 


[kg 


"3 

“3 


Valid phases 


Reactor: 

| Vapor-Liquid-Liquid 

r| 

Vent accumulator: 

I Vapor-Only 

"3 


Select key component(s) to identify 2nd liquid phase 


2nd Liquid | 


Available components 


CL2 

HCL 

BENZENE 

CC6H5CL 


Selected components 


Cutoff mole fraction for 2nd liquid phase 


Cutoff value for sum of mole fractions of key components: 


Use to select key component(s) for identifying the 2nd liquid phase. 


Figure 10.22 RBatch phase specifications. 


10.9 WORKSHOPS 

Workshop 10.1 Two reactors are connected in series. In the first reactor two reac¬ 
tions take place in series: 
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*f Specifications | y Reactions Slop Criteria | Operation Times | Continuous Feeds | Controllers | 


Stop criteria 



Criterion no. 

li "31 

1 “d 

1 3 


Location 

Reactor 

Reactor 



Variable type 

Mole fraction 

Time ^ | 



Stop value 

0.1 

Time 

Mole fraction 

Mass fraction 
Conversion 

T otal moles 

T otal mass 

T otal volume 

T emperature 



Component 

BENZENE 



Substream 

MIXED 



Property set ID 




Approach from 

Above 



Pressure 
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Figure 10.23 RBatch stop criteria. 
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Figure 10.24 RBatch operation time. 


The first reactor operates at 300°F and 300 psi. The first reaction conversion, based on 
nC 4 H 10 , is 50%. The second reaction conversion, based on C 6 H 6 , is 90%. In the second 
reactor, which operates at 400°F and 300 psi, the remaining ethylene is completely 
converted to 1 -butene by the reaction 


2C 2 H 4 -> C 4 H 8 
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Overall Composition | Feed | Properties | 


Composition profiles 

View:! Reactor liquid molar composition 



Time 
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HCL 

BENZENE 

C6H5CL 
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0 

0 

0 

1 

0 


30 
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0.05277989 
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240 

0.31109388 

0.22736766 
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0.24092574 
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Figure 10.25 Product distribution at various times. 


TABLE 10.1 Reactor Train Feeds 


Ethane, 

Ethylene, 

n -Butane, 

Benzene, 

1-Butene, 

Ethylbenzene, 

c 2 h 6 

C 2 H 4 

C 4 H 10 

c 6 h 6 

c 4 h 8 

CsHjo 

Feed 1 

10 

90 




Feed 2 



50 




The feeds for the reactor train, in lbmol/hr, and the component list are given in 
Table 10.1. The pressure of both feeds is 300 psia, and both are saturated liquids. 
Calculate the complete material and energy balance. 

Workshop 10.2 The following describes a process for making toluene from n- 
heptane. A fresh feed of 20 lbmol/hr of pure n -heptane at 77°F and 1 atmosphere is 
combined with a solvent recycle from an extractor and heated to 425°F at 1 atmosphere. 
The hot stream is fed into a reactor in which the following reaction occurs: 
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The conversion based on n -heptane is 15%. The products of reaction are cooled to 
180°F, after which the hydrogen is completely separated from the reactor products in 
the first separator. A feed of 100 lbmol/hr of benzene at 180°F and 1 atmosphere is 
combined with the remaining products of reaction to extract the toluene. All of the 
toluene and benzene leave as the product of the process. The unreacted n -heptane is 
recycled to the mixer. A sketch of the process is given in Figure 10.26. Using the 
RStoic model for the reactor, calculate the complete material and energy balance. 

Workshop 10.3a Production of hydrogen from methanol involves two steps, 
methanol re-forming and the shift reaction: 

CH 4 0 -> 2H 2 + CO 
CO + h 2 o -> h 2 + co 2 

Using an REquil block, find the temperature that maximizes hydrogen formation by 
re-forming at 1 atmosphere. 

Workshop 10.3b Repeat Workshop 10.3a using a RGibbs block. Compare the 
results. 

Workshop 10.3c Using 450°F as the operating temperature, determine the effect of 
pressure on the reaction. Select the optimal operating conditions. Is there a problem 
with this analysis? 

Workshop 10.3d Choi and Stenger (2002) have suggested that there is a possibility 
that dimethyl ether and methyl formate may be side products of the re-forming and 
shift reactions in Workshop 10.3a. Assess this possibility using the RGibbs block at 
450°F and a variety of pressures. 

Workshop 10.4 The formation of chlorobenzene as described at Examples/RPlugl 
Example in Section 10.8 is not limited to chlorobenzene formation but has the possi¬ 
bility of forming both p -dichlorobenzene and trichlorobenzene. Belfiore (2003) reports 
that the relative rates of reaction of the dichloro- and trichlorobenzen relative to the 
monobenzen are 1/8 and 1/243, respectively. Determine the quantities of the minor 
products formed for the conditions of RBatchExamplel. 
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Workshop Notes 

Workshop 10.1 The choice of serial or parallel reaction is made via a check box in 
the lower left-hand corner of the reaction setup page. 

Workshop 10.2 Initialization is with Aspen Plus defaults (i.e., zeros); however, after 
one iteration the error messages disappear and the process converges. 

Workshop 10.3a Use a sensitivity study, varying the temperature. 


Sensitivity S-1 Results Summary 



VARY 1 MEGIBBS PARAM TEMP F 
Figure 10.27 Temperature for maximum hydrogen production. 



VARY 1 MEGIBBS PARAM PRES PSIA 


Figure 10.28 Pressure effect on hydrogen production. 
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Workshop 10.3b The results are virtually identical to those of Workshop 10.3a. One 
may observe from Figure 10.27, a plot of the results from a sensitivity study, that 
450°F is the optimal operating temperature. 

Workshop 10.3c The effect of pressure shown in Figure 10.28, a plot of the results 
from a sensitivity study, is pronounced. Based on these results the optimal operating 
conditions would be 450°F and 14.7 psia. This analysis did not consider pressure 
effects at other temperatures, and one might find another optimal solution. This point 
is revisited in Chapter Thirteen. 

Workshop 10.3d Results of the case study show that low pressure favors hydrogen 
formation and high pressure favors the formation of small quantities of dimethyl ether. 
Methyl formate formation is negligible. 
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CHAPTER ELEVEN 


MULTISTAGE EQUILIBRIUM SEPARATORS 


Distillation and extraction are the most commonly used separation methods in chemical 
processing, and models of these are at the heart of Aspen Plus. All the Aspen Plus 
distillation and extraction models may be found in the model library under the tab 
Columns. 

For a distillation column, simulation implies the existence of a fully configured 
column (i.e., a number of theoretical stages) and the feed and sidestream locations. 
When the state of the feeds, and the overhead, bottoms, and sidestream specifications 
are provided, the model calculates the state of the products, and the reboiler and 
condenser heat loads, if they were not explicitly specified. For an extraction column 
the same ideas apply except that no overhead and bottoms specifications are required. 
There is limited capability in Aspen Plus to design a column directly to yield products 
of specified composition. Thus, the technique is to search the solution space for a trial 
column with the characteristics desired, and therefore many rating solutions may be 
required until one that matches the desired design specifications is obtained. 

11.1 BASIC EQUATIONS 

All models of distillation and extraction processes involve the solution of the material, 
equilibrium, and energy balance equations or an applicable subset. All models are 
based on the analysis of a single stage, such as that shown in Figure 11.1. There are 
several ways to write the describing equations, based on the choice of independent 
variables. For the set of equations below, the choices v[ and //, the component i vapor 
and liquid leaving stage j ; f j , the componential feed to stage j ; sj , a componential 
sidestream leaving stage j , Ti and P ; , the temperature and pressure of stage j, and <2 ; , 
the total energy entering or leaving stage j, are used. The phase equilibrium constant 
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F J fJ; 


hi lj-1. 


Vf yjj 


Q j . 


V J+ ' V;j +1 


U, |ij 


Stage j at V and P ] 

Figure 11.1 Theoretical stage. 


Sh sh 


K- depends on the stage temperature, pressure, and liquid and vapor mole fractions, 
xj and yj . Mole fractions can be expressed in terms of componential molar flows; for 
example, a liquid mole fraction can be calculated by 


rJ - 


J2k=i,j 


A 


(n.i) 


where m is the number of components. V 7 and L 7 are auxiliary variables which refer 
to the total molar vapor and liquid leaving stage j and are calculated by summing the 
componential flows in the vapor and liquid, leaving a stage as shown in 


Vj = E v i 

k=l,m 

l j = ^2 l l 


H j and h j are the molar enthalpies of the vapor and liquid leaving stage j, which 
depend on composition, temperature, and pressure. 

The material balance for component /, on stage j, is given by 

.// + '■/ 1 I //" -■/-//• v/ = 0 (11.3) 

where sj refers to the molar flow of component i in a sidestream, if it exists. There are 
nm material balances, where n is the number of theoretical stages and m is the number 
of components. The introduction of a sidestream adds one degree of freedom since 
the state of a sidestream is equal to the state of its source (i.e., the phase of the stage 
from which it is removed), but its flow rate is unknown; however, the composition, 
temperature, and pressure of the sidestream are known. A specification for a liquid 
sidestream, such as the ratio of the sidestream to the liquid, removes the extra degree of 
freedom and permits the added calculation of the componential flows of the sidestream. 
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The basic equilibrium equations for component i , on stage j , are 

yj = Kjxj (11.4) 

When the independent variables and equations (11.2) are substituted into (11.4), we 
obtain 




lj J2k=l,n 

;=1,j 


A 


(11.5) 


There are nm equilibrium equations. 

The enthalpy per mole of liquid and vapor flows H j and h j refer to the vapor and 
liquid leaving stage j , respectively. For each stage, one energy balance, given by 

h j f E + E <~ x +h j+l E l i +1 ~ Hj E v i~ hi E l L 

k=l,m k=\ ,m k=\ ,m k=\,m k= 1 ,ra 

- h J E V +Q j = 0 (11.6) 

k=l,m 

may be written, where S 7 represents the total molar flow of a sidestream. The 
enthalpy/mole of a sidestream depends on its source and may be written as either Hi 
or h{ depending on the state of the sidestream. If there are no feeds, sidestreams, or 
heat exchangers associated with stage j, the F, S , and Q terms disappear. There are 
n energy balance equations. 

Models usually do not include stage pressure calculations, and pressures are nor¬ 
mally assigned a priori. The summation of equations for a distillation column is as 
follows: 


Material balances nm 

Equilibrium equations nm 

Energy balances m 

Total m(2n + 1) 

In the absence of sidestreams or interstage heat exchangers, a count of unknowns 
for a column with a condenser and a reboiler is as follows: 


Componential vapor flows 
Componential liquid flows 
Stage temperatures 
Condenser and reboiler duties 
Total 


nm 

nm 

m 


2 

m(2n + 1) + 2 


Thus, there are typically two degrees of freedom, which correspond to the column’s 
end specifications. These may be reboiler and condenser duty specifications. Other 
specifications that can be used, such as product rate and reflux ratio, require minor 
modifications to the basic equations. 

Each sidestream adds one degree of freedom; therefore, it is necessary to provide 
two end specifications for each sidestream, plus a specification such as a sidestream 
flow rate. 
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For extraction columns the equations are identical to those of distillation columns if 
one recognizes that the vapor terminology now refers to the lighter of the two phases. 
Additionally, since there are no reboilers or condensers, the degrees of freedom are 
zero except that one additional degree of freedom is introduced for each sidestream, 
as in distillation columns. 

The set of equations described above lend themselves to solution simultaneously 
by means of the Newton-Raphson method, such as the algorithm of Naphtali and 
Sandholm (1971). This method has the advantage of permitting the use of all possible 
end and sidestream specifications; however, it has the disadvantage of requiring starting 
estimates for all independent variables. 

There are many arrangements of the equations with an accompanying convenient 
selection of variables that lead to a variety of algorithms to generate a solution. For 
example, one could employ total flows and mole fractions. There are a group of algo¬ 
rithms that alternately solve the locally linear material balances and then the nonlinear 
energy balances, adjusting the K values between iterations. The method used in Aspen 
Plus’s Distil block is the inside-out algorithm of Boston and Britt (1974). Details of 
many of the methods may be found in Seader and Henley (1998). 

11.2 THE DESIGN PROBLEM 

To simulate a distillation column within a process it is necessary initially to solve a 
design problem: that is for a given service, to establish the number of stages required 
and the location of the feeds and sidestreams, and to define suitable end specifications. 
Since the rigorous distillation models in Aspen Plus do not calculate directly such 
design parameters as the number of stages and the feed location, it is necessary to 
explore many possible rating solutions. Prior to the use of a rigorous model, it is useful 
to estimate the design by less rigorous methods, including the manual McCabe-Thiele 
method (1925). Prior to any design work, it is imperative that the vapor-liquid equilib¬ 
rium diagrams and/or liquid-liquid equilibrium diagrams be reviewed. In many cases 
the data stored in the Aspen Plus database will be suitable, but this may not be the 
case and a literature search or experimental work may be required. 

A brief review of the McCabe-Thiele method is given below. Unlike the models 
in Aspen Plus, the McCabe-Thiele method can be used to solve the design problem: 
that is, given the feed composition, flow rate, thermal state, and product composition 
targets, determine feed location and the number of stages above and below the feed 
stage. The following features are at the heart of the method. 

1. All stages are at steady state and the stage products are at equilibrium. 

2. The column operates under equal molal overflow conditions; that is, within a 
section of a distillation column where there is no feed, no product takeoff, and 
no interstage heat exchangers, the flow rates of all liquids leaving a stage are 
identical, as are the flow rates of all vapors. 

3. After a solution has been obtained, the reboiler and condenser duties are calcu¬ 
lated by a manual energy balance around each device. 

The McCabe-Thiele method can also be used to rate columns by a trial-and-error 
procedure. 
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An overview of the method is given with reference to Figure 11.2, the 
McCabe-Thiele diagram. The coordinates are the vapor and liquid mole fractions. 
An equilibrium curve based on experimental data is plotted in the figure. There are 
operating lines above and below the feed which are derived from the material balance 
equations in terms of the vapor and liquid mole fractions. The equations of the two 
operating lines are given in the figure. The two operating lines intersect with the 
y = x line at the column’s respective end compositions and intersect with a feed line 
that emanates from the feed composition, Zf, located on the y = x line. The slope of 
the feed line, q , is calculated by 


Hf-Hp 
H f - h f 


(11.7) 


where Hp represents the molar enthalpy of the feed, Hf is the molar enthalpy of 
saturated vapor, and hf is the molar enthalpy of saturated liquid. For example, if Hp 
is a saturated liquid, the value of q is zero, which represents a horizontal line on the 
McCabe-Thiele diagram. 

The number of theoretical stages is obtained by drawing a vertical line from an 
operating line to the equilibrium curve. This is equivalent to finding the equilibrium 
compositions for a stage. A horizontal line from the equilibrium curve back to the 
operating line gives the compositions of the passing streams (i.e., the material balance). 
This is repeated from one end of the column to the other. 

An important feature is the ability to determine the minimum reflux ratio, LIV , which 
occurs at an infinite number of stages. This is known as the pinch point and usually 
occurs when a line drawn from the overhead composition specification is connected 
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to the feed line/g-line intersection, as shown in Figure 11.3. If the equilibrium curve 
has an unusual shape, as in a methanol-water system, a different pinch point may be 
located, as shown in Figure 11.4. A rule of thumb commonly used for the design of 
distillation columns is that the operating reflux should be between 1.25 and 2.0 times 
the minimum. An estimate of the minimum number of stages required, which occurs 
at total reflux, is shown in Figure 11.5. 

If a sidestream is required, a third operating line is introduced, as shown in 
Figure 11.6. This operating line starts at the feed/g-line intersection and has a slope 
L'/V'. After determining the operating reflux ratio, a McCabe-Thiele diagram can 
be constructed and theoretical stages stepped off as described above. A detailed 
McCabe-Thiele example follows. 

11.3 A THREE-PRODUCT DISTILLATION EXAMPLE 

A 20 mol % ethanol-water solution is to be distilled. Two products, 80 mol % and 
40 mol %, are required. The feed is available as a saturated liquid. The composition 
of the bottom product should not exceed 0.25 mol % ethanol. The external reflux ratio 
is to be four times the minimum. The shape of the equilibrium curve necessitates a 
relatively high reflux ratio, The sidestream draw-off is 10 lbmol/hr. Using a basis of 
100 lbmol/hr, determine the number of theoretical stages required, the feed location, 
the sidestream location, and the quantities of the products. The experimental data of 
Kirschbaum and Gerstner (1939) are to be used. 

The data are plotted in Figure 11.7, as is the minimum L/V of 0.5, which corresponds 
to a minimum external reflux ratio LID of 1.0. The enriching section operating line’s 





y, mole fraction vapor y - mole of fraction 
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Figure 11.4 Minimum reflux with nonideal equilibrium. 



Figure 11.5 Minimum number of stages. 














mole fraction 
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Figure 11.6 Operating lines with a sidestream. 



Figure 11.7 Minimum and operating L/V. 
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L/V of 0.8 is derived from the external reflux ratio specification, which is also shown 
and corresponds to an operating external reflux ratio of 4. 

An overall material balance yields the following results: 

Distillate, D = 24.8 lbmol/hr, xp =0.8 
Sidestream, S = 10 lbmol/hr, = 0.4 
Bottoms, B = 65.2 lbmol/hr, xb = 0.25 

The various liquid and vapor flow rates throughout the column, based on equal molar 
overflow, are given in Figure 11.8, from which the operating line slopes are derived: 

Enriching section: L/V = 99.2/124 = 0.8 
Sidestream section: L'/V' = 89.2/124 = 0.72 
Stripping section: L"/V" = 189.2/124 = 1.53 

The preliminary design is given in Figure 11.9. There are a total of eight theoretical 
equilibrium contacts, one of which is the reboiler, as well as a total condenser. The 
feed is on stage 7 and the sidestream is on stage 6, counting from the top down. This 



Figure 11.8 Three-product distillation example with equal molal overflow. 
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is one of several possible solutions, since there is a trade-off between the number of 
stages, the reflux ratio, the equipment cost, and energy considerations. 


11.4 PRELIMINARY DESIGN AND RATING MODELS 


11.4.1 DSTWU 

The DSTWU model is applicable to both binary and multicomponent systems and 
was developed for systems that exhibit approximately constant relative volatility. The 
DSTWU model is based on the Gilliland correlation shown in Figure 11.10. This 
correlation was developed from the operating data of many columns and relates the 
minimum and operating reflux ratios, R m i n and R , to the minimum and actual number 
of theoretical stages, N m [ n and N. Gilliland’s data were correlated by Molokanov 
et al. (1972), which resulted in an analytical version of the Gilliland correlation: 


N - N min / 1 + 54AX X-l \ 

N + 1 eXP V11 + 117.2X X 0 - 5 ) 


( 11 - 8 ) 


where 


g _ R ~ ^min 

” R +1 

The minimum reflux ratio is obtained by applying Underwood’s (1946) method, which 
is applicable to binary or multicomponent systems. In the latter case it is necessary to 
select light and heavy key components and to solve for the roots of 
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E UjZj 
a i — 0 


-q 


(11.9) 


Here a * is the relative volatility of component i, Zi is the mole fraction in the feed of 
component i, q is the thermal state of the feed, using the McCabe-Thiele definition, 
and 0 is a root of the equation. To calculate the minimum reflux ratio, the value of 0, 
which lies between the light and heavy keys, is substituted into the equation 


Rn 


+ > = £ 


QLjXt 


( 11 . 10 ) 


where xP refers to the mole fraction of component i in the distillate. 

The minimum number of stages is calculated using Winn’s (1958) equation, an 
improvement on Fenske’s (1932) equation, where the usual definition of relative volatil¬ 
ity is replaced by 


K t = UKut 


( 11 . 11 ) 


and the parameters ft and 0* are determined from the equilibrium data for each com¬ 
ponent. Winn’s equation is 


In 


N min = 


(xd/xb)lk ( x b/xd)hk j 


In $lk 


( 11 . 12 ) 


The procedure for estimating the distribution of nonkey components is as follows. 
Assuming that the equation 


dj_ 

bi 



(11.13) 
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applies where d t , b t , d r , and b r are the distillate and bottoms componential flows of 
components i and r, and a^ r is their relative volatility, and when equation (11.13) is 
combined with the overall componential material balance equation, 


fi — di + bi 

where / is the feed of component i , the equations 

, fi 

bi = - 

1 + (dy jby} Oif y 

J fi (dy/by) Oii y 

di = - 

1 + (dy j by') di^y 


(11.14) 


(11.15a) 

(11.15b) 


are obtained. 

An example of the use of DSTWU is provided at Chapter Eleven Examples/DSTWU 
Example. The model’s primary input is shown in Figure 11.11. Selection of the Cal¬ 
culation Options tab enables a user to specify a table of reflux versus the number of 
theoretical stages. 

11.4.2 Distl 

The Distl model is applicable to both binary and multicomponent systems and makes 
use of the equal molal overflow and constant relative volatility approximations. The 
model employs Edmister’s method and rates single-feed, two-product columns. Details 
may be found in Seader and Henley (1998). 

The primary input form is shown in Figure 11.12. The DSTWU example was used 
as input to the Distl model. Results are given at Examples/DistlExample. 


^Specifications | ^/Calculation Options Convergence 


Column specifications 

. Number 
of stages: 
f* Reflux ratio: 


W 


Pressure 


Reboiler: 


|l 4.7 

psia ^ 

|l 4.7 

psia 


Key component recoveries 
Light key: 


Comp: 

METHANOL 

T3 

Recov: 

|0.95 


Heavy key: 



Comp: 

ETHANOL 

J 

Recov: 

|0.05 



Condenser specifications 

(* T otal condenser 

Partial condenser with all vapor 
distillate 

Partial condenser with vapor 
and liquid distillate 

Distillate vapor fraction: 10 


Figure 11.11 DSTWU primary input. 
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^/Specifications 


Convergence 


Column specifications 


Number of stages: 

14 


Feed stage: 

7 


Reflux ratio: 

6 


Distillate to feed mole ratio: 

0.295 


Condenser type: 

j Total 

3 


Pressure specifications 
Condenser pressure: 
Reboiler pressure: 


14.7 


14.7 


psia 
| psia 


73 

3J 


Figure 11.12 Distl primary input. 


11.5 RIGOROUS MODELS 

After developing a preliminary design, it should be possible to develop a rating model 
of the system of interest. If the system is nearly ideal, or the phases can be represented 
by an equation of state, the Aspen Plus database is probably suitable. If, on the other 
hand, the system is sufficiently nonideal that an activity coefficient equation is required, 
it is advisable to assess the equilibrium data within Aspen Plus by comparing them with 
experimental data. If the data do not compare favorably, fit suitable data to an activity 
coefficient model. If the Aspen Plus database is lacking, acquire suitable equilibrium 
data either by searching the literature, developing data experimentally, or estimating 
data using Unifac. Once suitable data have been found, they should be fit to an activity 
coefficient model. 

As an example, the vapor-liquid equilibrium data of Wilson and Simons’ (1952) 
for the isopropanol-water system at 95 mmHg is to be considered for a distillation 
application. The experimental data were compared against results calculated using a 
property analysis run “points along a flash curve” using Aspen Plus’s stored data fit 
with the Wilson equation. Figure 11.13 shows a plot of both data sets. 

The usefulness of Aspen Plus’s stored data should be assessed by an engineer in 
light of the precision required for the application. For example, If the application 
required that the product composition not be in the region approaching the azeotrope, 
the Aspen Plus data would suffice since the data calculated differ very little from the 
experimental data. If however, the proposed column required that the product be near 
the composition of the azeotrope, Figure 11.13 shows that the azeotrope composition 
calculated is about 0.05 mole fraction higher than the experimental value. This would 
produce serious problems, for example, if this column were part of an azeotropic 
distillation scheme involving two columns operating at different pressures. The quality 
of the 95 mmHg data calculated from the stored Wilson parameters is suitable since 
the data that were fit ranged between 77 and 212°F, as can be seen from Figure 11.14, 
and the experimental data range from a low of about 35°C (95°F) to a high of 49°C 
(120°C), but unfortunately, do not accurately model the composition of the azeotrope. 
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IPA - H 2 0 - 95 mmHg 



If no experimental data are available for a system, or data cannot be developed in 
a timely manner, the Unifac group contribution method may be used to estimate the 
data. Although these data may be used for preliminary design, it would be wise to 
set in motion means for experimental verification. For example, Figure 11.15 shows 
a comparison of experimental and Unifac estimated vapor-liquid equilibrium for an 
isopropyl alcohol-water system at 95 mmHg. The prediction of the azeotrope compo¬ 
sition is remarkable, but the lower composition range is acceptable at best, but again, 
use of these data is application driven. 

11.5.1 RadFrac 

RadFrac is the rigorous distillation modeling workhorse of Aspen Plus. It is based on 
the inside-out formulation of Boston and Sullivan (1974), which employs a reorgani¬ 
zation of the basic equations and the sequence of calculation. A detailed description 
of the method may be found in Seader and Henley (1998). For problems in which 
convergence difficulties are experienced, five variations of the basic algorithm are 
available. 

The first example of RadFrac use here is the same as that used for DSTWU and 
Distl. Since RadFrac is a rating or simulation program, the necessary specifications 
are obtained from the results of the DSTWU example (i.e., feed stage location, reflux 
ratio, and distillate rate). The primary input form for RadFrac is shown in Figure 11.16. 
Most of the remaining process specifications are straightforward and accessed from the 
tabs available. The tabs are Streams, in which the feeds and products are specified; 
Pressure, in which the stage pressures are specified; Condenser, in which the condenser 
details are specified; and Reboiler, in which the reboiler details are specified. If in 
Figure 11.16 the reboiler specified is thermosyphon, the Thermosyphon Config and 
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^Input ^Databanks 

Parameter: |WILS0N Dataset: fl 


T emperature-dependent binary parameters 



Component i 

Component j 



C3H80-01 

T 

WATER 

T 


T emperature units 

F 



Source 

APV70 VLE-IG 



AIJ 

■6.226100000 



AJI 

-.0133000000 



BIJ 

2418.577001 



BJI 

-239.2651701 



CIJ 

0.0 



CJI 

0.0 



DIJ 

0.0 



DJI 

0.0 



T LOWER 

77.00000338 



TUPPER 

212.0000023 



EIJ 

0.0 



EJI 

0.0 



Property units: 




Figure 11.14 Stored isopropyl alcohol-water Wilson parameters. 


Rebolier tabs become available. When the Thermosyphon Conhg tab is selected, the 
thermosyphon reboiler baffle configuration may be selected, as shown in Figure 11.17. 
When the tab Reboiler specification is chosen, Figure 11.18 appears. This permits 
various reboiler specifications to be made for either kettle or thermosyphon reboilers. 
If the Reboiler Wizard button is pushed, the Reboiler Wizard display, on the right 
of Figure 11.18, appears and enables details of an associated heat exchanger to be 
entered, including a link to the Heatx block. Results for this simulation may be found 
at Example/RadFracExamplel. 

The RadFrac block has an associated design specification function. In the example 
above, the DSTWU solution shows a 0.95 recovery fraction of methanol; however, 
the rigorous RadFrac solution shows that the recovery is 0.937. One possibility for 
improving the recovery is to vary a column specification, in this case either the distillate 
rate or the reflux ratio. Use of the design specification, varying the distillate rate, 




































160 MULTISTAGE EQUILIBRIUM SEPARATORS 


I PA- WATER- 95 mmHg 
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Convergence: 
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Figure 11.16 Primary RadFrac input. 
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V Configuration | </Streams | Pressure j </Condenser Thermosiphon Config. j© Reboiler | 3-Phase 



Figure 11.17 Thermosyphon reboiler options. 


y Configuration | </Streams| ■/Pressure | ■/Condenser | Thermosiphon Config. Reboiler 


Thermosiphon reboiler options 

(• Specify reboiler flow rate 
C Specify reboiler outlet condition 
C Specify both flow and outlet condition 


Reboiler Wizard | 


Flow rate 

|Mole || |Ibmol/hr ▼ | 

(Temperature J 


1-F-3 


Optional 
Reboiler pressure: 

Reboiler return feed convention: 
Utility: 


| psia ▼ | 


Above-Stage 3 



f Reboiler Wizard 


Reboiler Type: 

Pseudo Stream ID: 14 

rHeatx specifications 

- 1 


Block ID: | 


Type: | giHIHMUgg 

Mode: | Design 

Circulation rate: 

Shell&Tube input file with full path and file extension: 

Browse... | 


Hash Block ID: | ▼ 

1st Calculator Block ID: |c-l 


I - Move to hierarchy block ["" 

I* Maintain link to these blocks 

OK | Cancel | 


Help | 


Figure 11.18 Thermosyphon reboiler specifications and wizard. 


is illustrated at RadFracExample2. Figure 11.16 shows a list of options under the 
configured block named Distil. The entry Design Specs, followed by Object Manager, 
New, and Mole recovery target of 0.95, will create the specification, which for this 
example is shown in Figure 11.19. The selections of Components and Feed/Product 
Streams allow completion of the specification. Figure 11.20, created by using the Vary 
entry, shows the bounds within which to search for a result. Results are given at 
Examples/RadFracExample2. 

















































































162 MULTISTAGE EQUILIBRIUM SEPARATORS 


^Specifications ^/Components Feed/Product Streams 


Options 


Design specification 
Type: 


Mole recovery 


"3 


Specification 
T arget: |0.95 


Figure 11.19 RadFrac design specification. 


7 S pecif ications | 


Adjusted variable 


Type: 

|Distillate rate ^ | 


l-'l- 1 - 1 al lu IUVYGI IU-j- 

Lower bound: 

Upper bound: 


100 | Ibmol/hr 

200 | Ibmol/hr 


Optional 

Maximum step size: 


Figure 11.20 RadFrac vary form. 


RadFracExample3 illustrates the use of design specifications in conjunction with the 
normal end specifications of RadFRac. In this case the distillate rate and reflux ratio are 
permitted to vary, with the aim of achieving composition targets of the distillate and 
bottoms composition of 0.95 mole fraction methanol and 0.01 mole fraction methanol, 
respectively. The solution shows that for the limits of the perturbed variables the 
bottoms composition is obtained and the overhead specification is close and probably 
acceptable. Note that the Aspen Plus results show many errors, but in a practical sense 
they are meaningless. These results are obtained by using two implementations of the 
Design Spec and Vary functions. The final values of the reflux ratio and the distillate 
rate are given in the report. 

The remaining design questions are: Are the number of stages nearly optimal? Where 
is the best feed location? RadFrac addresses these questions by selecting NQ Curves 
under the list of options for the entry Distil. The setup of this capability is shown in 
Figure 11.21 and is illustrated at Examples//Example3a. In this example, many errors, 
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^/Specifications Feed / Product j PA/Decanter | 
T otal stage optimisation 

F 


Lower limit: 
Step size: 


F 


Upper limit: |15 — 7 -j 

Start with the maximum number of stages 


Feed tray optimization 
Feed stream: 

Minimum reflux ratio: 
Stage limit from top: 
Stage limit from bottom: 


>i 


F~ 


Objective Function 

Minimize: jQreb-Qcond U Reboiler/Condenser cost ratio: 


Qcond 

Qreb 

Mole-Rr 

Mass-Rr 

Stdvol-Rr 


|T otal heat load of condenser and reboiler (adjusted by reboiler/condenser cost ratio). 

Figure 11.21 NQ curves specifications. 


some of which are shown in Figure 11.22, are generated. The various reasons for the 
errors include limitations at the bounds of perturbed variables. In this case the optimal 
feed location and number of stages required were determined correctly from DSTWU 
and are given in the base case solution. This will usually not occur unless the species 
involved form relatively ideal mixtures. 

Aspen Plus’s RadFrac block provides many detail design functions, which are acces¬ 
sible through the Tray and Packing Sizing and Rating folders under the RadFrac block. 
These functions are based on vendor-provided information or well-known methods 
described in the literature. Details of the calculations, correlations, and so on, are 
given in Aspen Plus’s Help documentation, Appendix A. An example of the tray design 
input is given in Figure 11.23. Note that five different tray types are available. In this 
example the sizing function was used for three sections of the column: enriching, feed, 
and stripping. Results are shown at Examples/Example3b. When using RadFrac for 
rating or simulation, the design specifications should be removed and operating spec¬ 
ifications such as reflux ratio and distillate rate should be used. Additionally, the tray 
or packing design specifications should be removed and replaced with the equivalent 
tray and packing rating functions. An example of such an implementation is shown at 
Examples/Example3c. Figure 11.24 shows an example of the packing input. 
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Iteration = 2 Feed Stage = 6 Objective Function = 0_37617E-K}8 


Convergence 
OL ML 

1 4 

2 2 

3 2 


iterations: 

IL Err/Tol 

12 73.239 

4 3.9229 

3 0.96062 


DESIGN SPEC IS NOT SATISFIED BECAUSE ONE OR MOBE MANIPULATED 
VARIABLE IS AT ITS BOUND. 


Iteration = 3 Feed Stage = 4 Objective Function = 0 _ 37618E-I-G8 

Feed tray cptiir.isaticn converged for NSIAGI = 8 

Optimum feed stage = 5 Objective Function = 0.37611E+U3 

WARNING 

NQ CURVE ™1™ STOPPED AT STAGE 11 8" SINCE THEBE IS NO SIGNIFICANT 
CHANGE IN OBJECTIVE FUNCTION WHEN INCREASING NUMBER OF STAGES . 


** Converging the base case for NSTAGE = 15 ** 

Convergenee iteratiens: 

OL ML IL Err/Tol 

1 5 14 0.3S203E-03 

*** END GENERATION OF NQ CURVES *** 

->Simulaticn calculations completed . . 


Figure 11.22 Errors generated by NQ run. 
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pH 
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Minimum, l etter Float Valve 


Cap slot area to active area ratio: 
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: 1 i 

- 

Traytype: Sieve 

▼ | Number of passes: 

j 

, 1 1 


2 
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1 

ft 

0.12 


jo. 12 


Figure 11.23 Tray design input. 


11.5.2 Extract 

Extract is a rigorous extraction column simulation block. The material, energy, and 
equilibrium equations given in equations (11.1) through (11.6) apply with the under¬ 
standing that the vapor-phase nomenclature is fully appropriate for interpretation as 
the second liquid. Because there are no phase changes, there are rarely large energy 
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^Specs Design / Pdrop Layout Downcomers 


T rayed section 


Starting stage: 2 


Ending stage: 

F 

Tray type: Sieve 


t | Number of passes: 

F 



T ray geometry 


Diameter: 

T ray spacing: 

7 

P Zf 

Deck thickness: |l[l GAUGE ^ 

2 
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Panel B 
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|! 3 

ft 

d 

d 
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1 
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Figure 11.24 Packing input. 


'/input [^Databanks] 


Parameter: |UNIQ Dataset: [i 

T emperature-dependent binary parameters 



Component i 

Component j 








|METHANOL ▼ | 

| METHANOL ^ | 

Imethanoi -1 

1 ETHANOL »J 

| ETHANOL *J 

|WATER 

1 lJ 

|ETHANOL 

jwATER 

|toluene 

jwATER *J 

|TOLUENE 

|T0LUENE ^J 

1 JlJ 


T emperature units 

C 

C 

c 

C 

C 

C 



Source 

APV70 VLE-IG 

APV70 VLE-IG 

APV70 LLE-LIT 

APV70 LLE-LIT 

APV70 LLE-LIT 

APV70 LLE-LIT 



AIJ 

•2.650900000 

-1.066200000 

0.0 

0.0 

0.0 

0.0 



AJI 

1.289100000 

.6437000000 

0.0 

0.0 

0.0 

0.0 



BIJ 

651.4882000 

432.8785000 

-178.5900000 

185.3000000 

189.1900000 

-350.2100000 



BJI 

-273.6917000 

-322.1312000 

-397.5800000 

167.3800000 

-357.0300000 

-950.6000000 



CIJ 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 



CJI 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 



DIJ 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 



DJI 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 



T LOWER 

20.00000000 

24.99000000 

20.00000000 

20.00000000 

20.00000000 

20.00000000 



TUPPER 

78.40000000 

100.0000000 

40.00000000 

40.00000000 

40.00000000 

40.00000000 



EIJ 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 



EJI 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 



Property units: 









I - Estimate missing parameters by UNIFAC 


Figure 11.25 Uniquac parameters and data sources. 


swings in extraction columns, and thus they operate virtually isothermally. There are 
no condensers or reboilers, and therefore the number of degrees of freedom is equal 
to the number of available equations, and no operating specifications are required. The 
Extract block is a rating or simulation model, and its use requires the specification of 
feed rates, states, and locations as well as the number of equilibrium stages. 

When preparing an extraction column simulation, a review of the available data 
in Aspen Plus’s databases is imperative. Figure 11.25 shows the Uniquac parameters 
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for the components of the extraction example given at Examples/ExtractionExample. 
These parameters are located in the subheadings under Parameters/Binary Interaction/ 
Uniq-1. Figure 11.25 shows that the source of most of these parameters is APV70LLE- 
LIT (liquid-liquid equilibrium data) but that the data for the methanol-water and 
ethanol-water systems are derived from vapor-liquid equilibrium data. The choice of 
the source of the data is shown in Figure 11.26. The previous caution about data source 
applies here. Are quaternary or ternary data available in the literature for this system? 
Do multicomponent results generated by Aspen Plus’s stored data agree with experi¬ 
mental data? Note the temperature limitations on the data in Figure 11.25. Simulation 
work carried on outside the applicable temperature ranges of the data have potential 
for substantial error. 


/Input /Databanks 


Select databanks and their search order 


Available databanks 


APV70 VLE-FiK 
APV70 VLE-HOC 


<< 


Selected databanks 


>> 


APV70 VLE-IG 
AFV70 LLE-ASPEN 
APV70 VLE-LIT 


JL 


*1 

♦ I 


[Databank from which you want to retrieve parameters. 


Figure 11.26 Data source choices. 
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11.6 BatchSep 

Batch distillation is commonly used in processing specialty chemicals and pharmaceuti¬ 
cals and has a substantial role in the chemical industries. Aspen Plus’s batch distillation 
block is BatchSep. The equations describing batch distillation are very similar to those 
describing steady-state distillation; however, the material and energy balances are time 
varying, and therefore differential equations are required. An example of a differential 
material balance analogous to equation (11.3) is 



Here the alternative form of the material balance is used, with the independent variables 
being mole fractions and total flows. If M 7 is the volumetric stage holdup and p is 
the molar density, equation (11.16) results. Similarly, a differential energy balance 
analogous to equation (11.6) is 


d(pM J c p T) 
dt 


= H J F F j + HV j+l +hL j ~ i - HV j ~hV - hS J + Q (11.17) 


For an ft-component column, there are n differential material balances and one 
differential energy balance per stage. Additionally, n nonlinear algebraic equations 
describing phase equilibrium and the auxiliary summation of vapor and liquid mole 
fractions equations are also required per stage. This large set of combined differential 
and nonlinear equations can be solved numerically. Initial conditions are required for 
all the variables described in differential equation terms. Targets during the distillation, 
such as overhead temperature, are used to identify the beginning and end of a cut, and 
receivers must be defined for the collection of each cut produced. 



{U o 





Figure 11.27 Batch distillation column. 
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Setup 
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Initial Conditions 
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Operating Steps 
Final Conditions 
Solver Options 
Charge Stream Result 
Pot Results 
Profile Results 
Condenser Results 
Second Condenser Re 
Reflux Results 
Distillate Results 
Holdup Summary Resc 
Controller Results 
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Multiple Batches Resul 
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Configuration Pot Geometry PotHeatTransfer Condenser Condenser 2 Reflux Reflux Level 
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13 
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C3H8-01 
C^H 10-01 
ETHANOL 
METHANOL 



Figure 11.28 BatchSep basic setup. 


An example of a batch distillation application is given at Examples/BatchSep 
Example. Figure 11.27 shows a sketch of a batch distillation column. Stream 1 is 
the batch charge, streams 2, 3, and 4 are the products (cuts), and stream 5 is the 
pot contents at the conclusion of the distillation. Each of the product streams has an 
associated receiver. BatchSep applications require Aspen Plus release 7.1 or higher. 

The Setup form for the BatchSep block is shown in Figure 11.28. Note the possi¬ 
bility of two liquid phases in equilibrium with a vapor phase. Each of the available 
tabs enables definition of the physical dimensions of the still as well as the column’s 
overhead specification, such as the reflux ratio. The Jacket Heating tab enables the 
choice and value of the bottoms specification. The tab Receivers permits definition 
of the number of receivers and the association of each product stream with a par¬ 
ticular receiver. The Side Draws tab permits the entry of all necessary side draw 
information. 

The Controllers folder facilitates the use of a simulated PID-type controller during 
batch distillation. Figure 11.29 shows the two tabs that permit a user to choose a 
variable to be controlled and the controller settings that operate the controller. This 
capability enables one to simulate the real-time behavior of the batch still. 

The input forms for the Initial Conditions tab are illustrated in Figure 11.30. Both 
the Main and Initial Charge tabs are shown. 

Control of the products of the distillation is done with the Operating Steps entry. 
In this example there are two operating steps to the process. Figure 11.31 shows the 
setup of Stepone. When all the steps have been executed, the Final Conditions entry 
permits the choice of disposition of the contents of the reflux drum, the packing, and 
the pot. 
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Q Connections Parameters 


Process Variable(Measured Variable) 
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■j 
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Stage: 
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3 
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Z3 
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Figure 11.29 Batch distillation column controller setup. 



Figure 11.30 Initial conditions for BatchSep. 
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Changed Parameters 


End Condition 


Step end condition: 
Duration: 


Trigger value 


Trigger 

Variablelocation: 

Distillate receiver: 

Variable: 

Component 
Approach from: 
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Distillate receiver 


±1 
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V 
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V 
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V 
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V 
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Figure 11.31 Operating step input, Stepone. 


11.7 WORKSHOPS 

Workshop 11.1a A mixture of 50 mol/hr of benzene, 50 mol/hr of toluene, and 50 
mol/hr of paraxylene is to be distilled at 1 atmosphere. The feed is a saturated liquid 
at 1 atmosphere. The column is to operate at a reflux ratio of 2.8, and a 1-psi pressure 
drop across the column is anticipated. Ninety-five percent of the benzene fed and no 
more than 5% of the toluene fed is to be taken overhead. Use the shortcut distillation 
block DISTWU to estimate the column performance. Generate a table of reflux ratio 
versus theoretical stages required for this separation. Create a mixed unit set derived 
from ENG, replacing °F with degrees °C and psi with mmHg for all Workshop 11.1 
variations. 

Workshop 11.1b Repeat Workshop 11.1a using the physical properties option 
NRTL. Compare to Workshop 11.1a. 

Workshop 11.1c Solve Workshop 11.1a rigorously using Radfrac. Employ 12 the¬ 
oretical stages with feed on stage 6. Specify that 50 mol/hr of distillate is required. 
Solve using the ideal property option. Compare to the results of Workshop 11.1a. 

Workshop 11. Id Repeat Workshop 11.1c, but use the property option NRTL. Com¬ 
pare to the Workshop 11.1c results. 

Workshop 11.2a A mixture of 50 mol/hr of methanol, 50 mol/hr of ethanol, and 
50 mol/hr of water is to be distilled at 1 atmosphere. The feed is a saturated liquid at 
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1 atmosphere. The column is to operate at a reflux ratio of 1.2, and a 1-psi pressure 
drop across the column is anticipated. Ninety-five percent of the ethanol fed, and no 
more than 5% of the water fed, is to be taken overhead. Use the shortcut distillation 
block DISTWU with the physical properties option Wilson to estimate the column 
performance. Generate a table of reflux ratio versus theoretical stages required for this 
separation. Create a mixed unit set derived from ENG, replacing °F with °C and psi 
with mmHg for all Workshop 11.2 variations. 

Workshop 11.2b Solve Workshop 11.2a rigorously using the block RadFrac. Specify 
that 50 mol/hr of bottoms is required. Use 15 theoretical stages with feed on stage 10. 
Plot the liquid compositions for the column. What do you observe? Compare to the 
results of Workshop 11.2a. 

Workshop 11.2c Repeat Workshop 11.2b rigorously. Specify that 100 mol/hr of 
distillate is required. Use 11 theoretical stages with feed on stage 7. Compare to the 
results of Workshop 11.2b. 

Workshop 11.3a The data for this problem, from Udovenko and Mazanko (1964), 
are available in the file DataEleven-3bkp in the Workshops folder. Create a set of units, 
NEWSET, which employs pressure in mmHg and AP in atm, copied from the data 
set MET. Fit the data to the NRTL equation, assigning 0.0 to the A; j parameters and 
assigning 0.2 to each C parameter. Produce a triangular plot of the data. 

A mixture of 500 lb/hr of benzene and 500 lb/hr of 2-propanol is to be extracted 
with 1000 lb/hr of water in a continuous countercurrent extraction column operated 
at 1 atmosphere and 30°C. The feeds are liquids at 16 psi and 30°C. A 1-psi pressure 
drop across the column is anticipated. The column has four theoretical stages. Use 
the extraction block to calculate the column’s performance. Create a sensitivity study 
varying the quantity of water used for extraction from 700 to 1800 lb/hr in increments 
of 100 lb/hr. Create a set of units copied from ENG using °C. 

Workshop 11.3b Use a design specification to determine the quantity of water 
needed to recover 96% of the alcohol. 

Workshop 11.3c Use the decanter block to develop a four-stage batch extraction 
solution to Workshop 11.3b. Divide the amount of water calculated in Workshop 11.3b 
into four equal quantities so that each can be added to the batch four times. Compare 
the quantity of alcohol extracted to the solution for the continuous case. What do you 
conclude? 

Workshop Notes 

Workshop 11.1a through 11. Id All results are similar because this system is nearly 
ideal. Workshop 11.1c has an error message relating to the pressure of the feed stream. 
Note that the feed pressure, 14.696, is lower than the feed stage pressure. Workshop 
1 l.ld’s feed is set at 16 psi to fix the problem. 

Workshop 11.2b Results at reflux ratio of 1.2 with a bottoms takeoff of 50 lbmol/hr 
are similar to Workshop 11.2a. Figure 11.32 shows a plot of liquid composition as a 
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DATAO Triangular plot of BENZENE 2-PROH WATER 



Mole fraction BENZENE 
Figure 11.32 Triangular plot. 


function of the stage numbers. One may observe that the methanol composition barely 
changes between stages 6 and 11, which suggests that there are too many stages in the 
current design. 

Workshop 11.2c The column is reconfigured to 11 stages with the feed on stage 7. 
The results are very close to those of Workshop 11.2b. The column design as proposed 
by DSTWU has too many stages. 

Workshop 11.3a The triangular plot is shown in Figure 11.32. 

Workshop 11.3c To obtain the total quantity of 2-propanol extracted, it is convenient 
to set up a Mix block and feed it the aqueous phases from each decanter. When this 
is first run, one observes that with one-fourth of the required water added to the the 
organic feed, no phase splitting occurs. It is necessary to feed both products of the first 
decanter to the second decanter, and then everything proceeds as expected. If an error 
message referring to insufficient iterations arises, double the default value of 30. 

When comparing the results to Workshop 11.3b, one observes that less than half of 
the alcohol is extracted by batch extraction. 
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CHAPTER TWELVE 


PROCESS FLOWSHEET DEVELOPMENT 


Prior to this chapter, most of the Aspen Plus blocks have been introduced and the 
reader has gained experience in their use. This chapter deals with the integration of 
blocks into a process flow diagram (PFD), which is the source document for the process 
engineering function of producing process and instrument diagrams (P&IDs). A typical 
PFD contains process details of each processing unit and all reaction units, connected 
by the main process lines. Each unit operation and reaction unit will contain basic 
design data, such as operating temperature and pressure, and limited equipment design 
data, such as the number of theoretical stages required. Additionally, the presence of 
trace components will be identified, although it is not unusual that some surprises 
arise. Very few detail design data are included on a PFD, nor are instrumentation 
and process control requirements, operating instructions, equipment redundancy, safety 
considerations, shutdown and startup requirements, and the myriad other details that 
are the responsibility of the P&ID developers. 

The PFD development process begins with a process description typically based on 
small-scale chemical experiments and larger-scale minipilot and pilot-plant scale work. 
Additionally, a collection of physical and thermodynamic property data is developed 
by experimentation, collected from the literature, and estimated if necessary. These are 
used to produce a block diagram of the process, which is usually nothing more than a 
material balance that includes many performance estimates based on laboratory work 
and, in some cases, guesstimates. 

12.1 HEURISTICS 

When developing a process flow diagram, the following experience-based ideas are 
worthwhile: 
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1. A large flowsheet should be divided into a subset of smaller flowsheets. 

2. Each subset should have its own appropriate properties. 

3. Initially, in solving a subset, basic blocks without energy balances should be 
employed and converged. 

4. All initial attempts to converge should use Aspen Plus defaults (i.e., Wegstein’s 
method). The Wegstein approach can solve a wide variety of problems. 

5. Initial convergence should be with loose specifications. 

6. As the development proceeds, the basic blocks should be replaced, one or two 
at a time, with rigorous blocks, and energy balances turned on. 

7. The rigorous blocks should be tested on a stand-alone basis with feeds from the 
streams calculated using the basic blocks prior to integration. 

8. When a subset is executed using rigorous blocks, the independent variables from 
all tear streams that are feeds to the blocks should be transferred from the con¬ 
verged initial flowsheet where the basic blocks were used. 

9. If the Aspen Plus-selected tear streams within the subsets are not suitable, new 
tear streams should be identified. Convergence blocks should be defined for the 
new tear streams. New calculation sequences may also be required. 

10. When the subsets are completed, assembly of a fully rigorous flowsheet should 
proceed one subset at a time with intermediate testing. Tear streams may need to 
be changed. A gradual tightening of specifications is required. This process con¬ 
tinues until the combined subsets have been integrated and the entire flowsheet 
has converged. 

12.2 EXAMPLE: THE PRODUCTION OF STYRENE 

A flowsheet of a simplified version of the commercial process for the production of 
styrene by catalytic dehydrogenation of ethylbenzene is shown in Figure 12.1. The 
process description is as follows. 



Figure 12.1 Modified styrene production flowsheet. 
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1. Ethylbenzene is converted to styrene by the catalytic dehydrogenation reaction 

C 8 H 10 (g)^C 8 H 8 (g)+H 2 (g) 

The presence of steam in the reactor is known to suppress side reactions. 

2. Fresh ethylbenzene, stream 1, and a recycle, stream 15, containing primarily 
liquid ethylbenzene, combine at A, producing stream 3, which is heated to 500° C 
at B, producing stream 4. 

3. A recycle, stream 11, which is primarily water, is combined with stream 14, a 
freshwater makeup, at E, producing stream 13 at 50°C. 

4. Stream 13 is heated to 700°C, producing stream 5, which is combined with 
stream 4 at C, producing stream 6 at 560°C. 

5. Stream 6 feeds the reactor, F. Stream 7 exits the reactor at a temperature and 
pressure of 560°C and 14.7 psi, respectively. The reaction conversion is 35%. 

6. Stream 7 is cooled to 50°C at G, producing stream 9, which is primarily H 2 , 
and goes to another part of the plant. Stream 8 is cooled further to 25°C at H, 
to separate the water, stream 11, and organics, stream 10. 

7. Ethylbenzene and styrene in stream 10 are separated at J, producing the prod¬ 
uct stream 12, containing primarily styrene, and a stream containing primarily 
ethylbenzene, which is condensed at K, producing stream 15. 

The process model is to be based on the following: 

Stream 1: 100 lbmol/hr of pure ethylbenzene at 25°C and 14.7 psi 
Stream 14: 40 lbmol/hr of pure water at 25°C and 14.7 psi 

The entire process operates at approximately 1 atmosphere. 

12.3 A MODEL WITH BASIC BLOCKS 

The first model of this process was developed using split blocks for the separations. 
The process description above was used to estimate the split values and a simple 
conversion based block was used to describe the reaction. The model results are given 
at Chapter Twelve Examples/Example Twelve-1. Analysis of the results seem to agree 
with the process description. 

12.4 PROPERTIES 

There are three separations that require rigorous modeling in the development of 
a process flow diagram. All the necessary pure component data are available in 
the Aspen Plus data bank. Vapor-liquid, equilibrium and liquid-liquid equilibrium 
data are required to model the flash, the decantation, and the distillation elements 
of the process. The Uniquac binary parameters stored in Aspen Plus will be used 
as a source for the equilibrium data. Property analysis runs designed to present 
vapor-liquid equilibria for the three possible systems exclusive of hydrogen are 
shown at Examples/ExampleTwelvel-a. Figure 12.2 shows the binary parameter data 
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*/lnput VDatabanks 


Select databanks and their search order 


Available databanks 


APV70 VLE-RK 
APV70 VLE-HOC 


Selected databanks 


AFV70 LLE-LIT 
APV70 VLE-LIT 
APV70 VLE-IG 
AFV70 LLE-ASPEN 




J 

J 


Figure 12.2 Data bank sources. 


source choices for the systems, Figure 12.3 shows the Uniquac binary parameters 
for the three systems, and Figure 12.4 shows a plot of the ethylbenzene-styrene 
vapor-liquid equilibrium and hand calculations indicate relative volatilities between 
1.2 and 1.4. The other two systems exhibit mutual solubility in a very narrow range of 
composition, and because of this, only liquid-liquid equilibrium is available. These 
values should be suitable for use in the distillation calculations. 


12.5 RIGOROUS FLASH AND DECANTER 

A more rigorous version of the process may be found at Examples/ExampleTwelve-2. 
A two-phase flash and a decanter block were placed onto the flowsheet replacing the 
splitters at G and H, using the properties from above. Several exploratory runs were 
made, with the following results. With the two-phase flash at G at 50°C, the water in 
the process was pushed out at product stream 9, along with the hydrogen. Subsequently, 
the flash temperature was changed to 15°C. 
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0.0 
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Figure 12.3 Uniquac binary parameters. 



mole fraction liquid 


Figure 12.4 Ethylbenzene (l)-styrene (2) vapor-liquid equilibrium. 
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Summary Define Variable 
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2 

OK 
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3 

OK 

7 
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4 

OK 

9 
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0.69421955 
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0.07797347 


5 

OK 

11 

0.99919029 

0.G941559G 

299.294932 

0.05730372 

0.09494642 


6 

OK 

13 

0.99293443 

0.69410549 

303.592594 

0.07391347 

0.11121941 


7 

Errors 

15 

0.97323275 

0.69419504 

327.277253 

0.14095932 

0.120924 


Figure 12.5 Sensitivity study varying the water feed. 


A sensitivity study was made varying the amount of water entering the process for 
the purpose of determining the effect on the products and the composition and flow 
rate of the reactor feed. The water feed rate was varied from 3 to 15 lbmol/hr. The 
results are shown in Figure 12.5. The variables are defined below. 

• H2PROD is the mole fraction hydrogen in stream 9. 

• ORGPRD is the flow rate of water in lbmol/hr in stream 12. 

• FLOS6 is the flow rate of stream 6 in lbmol/hr. 

• XH20S6 is the mole fraction of water in stream 6. 

• XH20S9 is the mole fraction of water in stream 9. 

After about 13 lbmol/hr of water fed per 100 lbmol/hr of ethylbenzene, the material 
balance is unstable. 

1. The mole fraction of water in stream 6, the reactor feed (variable XH20S6), 
varies from about 0.011 to 0.074. 

2. The mole fraction of water in stream 9, the hydrogen product (variable XH20S9), 
varies from 0.023 to 0.11, which corresponds to a change in the composition of 
hydrogen (variable H2PROD), varying from 0.97 to 0.88. 

3. The total reactor feed flow rate (variable FLOS6) varies from 284 to 327 lbmol/hr. 

4. The flow rate of the organic product varies very little; however, rigorous distil¬ 
lation has not yet been implemented. 

The question of how much water to use in the process must be reviewed in light 
of the foregoing results. Assuming that it is desired to produce as pure a hydrogen 
product as possible, the remainder of the analysis will use a 3-lbmol/hr feed of water, 
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which corresponds to about 1% water in the reactor feed. The process description did 
not indicate what quantity of water is required in the reactor. Thus, 40 lbmol of water 
per 100 lbmol of feed is not suitable with the process description given. 

12.6 ANALYZING THE RIGOROUS DISTILLATION 

The distillation feed from Example 12.2 using a 3-lbmol/hr water feed was employed 
for the preliminary calculations in preparation for designing and then simulating the 
column. The feed with appropriate column specifications was employed with a DSTWU 
block to obtain an estimate of the number of stages, feed stage location, and reflux 
ratio required. These results may be found at Examples/ExampleTwelve-3a. The results 
were used to explore the design parameters of the column with several RadFrac 
runs. The results of these runs culminated in a Sensitivity study with Radfrac at 
Examples/ExampleTwelve-3b. The desired styrene product composition was about 0.98 
mole fraction. Figure 12.6 shows the effect of reflux ratio on the two product compo¬ 
sitions. Here S12XE and S2XS are the mole fractions of ethylbenzene in stream 12 
and styrene in stream 2, respectively. The column operates with 60 total stages, feed 
on stage 30, and a distillate rate of 183 lbmol/hr. The external reflux ratio required is 
about 9. Note that there is a slight amount of hydrogen in the feed; therefore, the use 
of a partial condenser is imperative. 

12.7 INTEGRATING THE RIGOROUS DISTILLATION 
INTO THE FLOWSHEET 

Introducing the converged rigorous distillation results of ExampleTwelve-3b into the 
converged flowsheet of ExampleTwelve-2 with the sensitivity study removed produces 
Example Twelve-4a. This example does not converge and produces several errors, 


S ummary D efine Variable 
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OK 

S 
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2 
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3 

OK 

8 

0.01544934 

0.01511328 


4 

OK 

9 

0.01126580 

0.0128G552 


5 

OK 

10 

0.00858002 
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Figure 12.6 Effect of reflux ratio on product composition. 
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Figure 12.7 Aspen Plus run control panel. 


some of which are displayed in the run control panel, a part of which is shown in 
Figure 12.7. The execution history shows that the RadFrac block and the Heater block 
that follows it have been bypassed. It would be wise to refer to the process flowsheet 
shown in Figure 12.8, where the failed blocks are clearly identified. What could be 
misleading is the message at the end of the execution history, which states that the 
loop $OLVER01 has converged. A look at the feed to the RadFrac block shows that 
it is significantly different from that of ExampleTwelve-3b. It is important that the 
remainder of the flowsheet produce the proper feed. The difficulty is that the sequence 
of calculations and the initial estimate of zeros for the variables of stream 6, the tear 
stream chosen by Aspen Plus, do not assure that the feed to the distillation column 
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from ExampleTwelve-3b enters the column. Because of this it will be necessary to 
alter the sequence of calculations and provide good estimates for the new tear stream. 


12.7.1 Selection of a Tear Stream 

If the process flow diagram (Figure 12.8) is examined, one can see that guessing the 
initial values of either stream 6 or stream 8 permits the sequential modular calculations 
to pass through each block, finally arriving at the stream values guessed; however, using 
stream 8 assures that only the organics pass through to the distillation column. One may 
then copy the converged values from Example Twelve-3b stream 8 onto the appropriate 
input of Example Twelve-4b. The definition of the tear stream, however, requires the 
definition of a convergence block for stream 8. This is accomplished by opening the 
Convergence entry under the list of specifications, opening the Convergence subentry, 
and selecting New with the object manager. Figure 12.9 shows the selection of tear 
stream with its convergence block. 


12.7.2 Sequence of Calculations 

The list of specifications also shows the entry Sequence, which when selected opens the 
object manager to enable the creation of a sequence of calculations by selecting New. 
The new execution sequence is shown in Figure 12.10. When this is executed, a satis¬ 
factory solution is obtained for each operation. Results are given at Examples/Example 
Twelve-4b. 


12.8 REACTOR FEED 

The process description, item 4, states: “Stream 13 is heated to 700°C, producing 
stream 5, which is combined with stream 4 at C, producing stream 6 at 560°C”. In 
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Figure 12.9 Tear stream selection. 
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Figure 12.10 User-defined execution sequence. 
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actuality it is necessary merely to provide 560°C as the set point to heaters B and D. 
This is implemented at Example Twelve-4b. 

12.9 MISCELLANEOUS CONSIDERATIONS 

Pumps are not used at this level of development. It is understood that the use of a mostly 
constant 14.7 psi pressure in this flowsheet is technically incorrect but introduces only a 
small error. The feed to the distillation column was changed to 18 psi to accommodate 
the column’s internal pressure. The actual size and location of the equipment has not 
been defined since the physical layout there has not yet been developed; therefore, 
pump selection cannot proceed since the actual service is unknown. For example, the 
distillation column has 60 theoretical stages, and if the efficiency is 100% and the tray 
spacing is 2 ft, the column is at least 120 ft tall (a 12-story building). If the site is in 
a windy area, this column may be built as two 60-ft columns with pumps in between. 

Similarly, the type of equipment used for heating some streams to high temperatures 
and cooling off others is unknown; thus, detailed design will be carried out by others. 

The amount of water used to control the formation of by-products is unknown. The 
reaction kinetics is unknown, and therefore at this point in the flowsheet development, 
a reaction conversion of 35% is used. Scale-up would be from laboratory data. 

12.10 WORKSHOPS 

Workshop 12.1 Figure 12.11 shows a section of a process for converting toluene 
to benzene in a hydroalkylation reactor. The main reaction is 


c 7 h 8 + h 2 ^ c 6 h 6 + ch 4 



Figure 12.11 Process for production of benzene. 
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TABLE 12.1 Process Feeds 


Component 

Feed 

(lbmol/hr) 

Recycle 

(lbmol/hr) 

Gas Recycle 
(lbmol/hr) 

Hydrogen 



2045.9 

Methane 



3020.8 

Benzene 


3.4 

42.8 

Toluene 

274.2 

82.5 

5.3 

Biphenyl 


1.0 


Temperature (°F) 

75 

250 

121 

Pressure (psi) 

569 

569 

569 


TABLE 12.2 Process Operating Conditions 


Location 

Pressure (psi) 

Temperature (°F) 

Furnace inlet 

569 

1000 

Furnace outlet 

494 

1200 

Reactor outlet 

494 

1268 

Pump outlet 

494 


Heater outlet 

480 

180 

Flash 

480 



An unavoidable side reaction that produces biphenyl is 

2C 6 H 6 ^Ci 2 H 10 + H 2 

The conditions for the feed and two recycle streams are given in Table 12.1. The 
specifications in Table 12.2 are the pressures and temperatures at various points in the 
process. 

The flow rate of the quench stream is to be such that the temperature of the com¬ 
bination of reactor effluent and the quench stream should be 1150°F. Conversion of 
toluene is 75%. Two percent of benzene present after the first reaction occurs is con¬ 
verted to biphenyl. The overall heat transfer coefficient in the heat exchanger is 60 
Btu/(hr-ft 2 -°F). Develop an approximate material and energy balance for the process. 

Workshop 12.2a Figure 12.12 shows a process to create a mixture suitable for use as 
a natural gas and recover a light hydrocarbon liquid from the feed given in Table 12.3. 
The feed is at 80°F and 237.7 psi. Develop a first estimate material and energy balance 
for this process. The equipment is to be modeled simply. The specifications for the 
processing equipment are as follows: 

Chiller : outlet at —40°F and 230 psi. 

Gas/gas heat exchanger : modeled as two heaters with a heat stream used to pass 
heat between them. 

Cold-side heater : outlet at 160 psi. A 4°F approach at the cold-stream outlet of the 
heat exchanger (i.e., the temperature of the inlet of the hot-side heater minus the 
temperature of the outlet of the cold-side heater). 
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TABLE 12.3 Component List 


ID 

Formula 

Databank Name 

Lb-mol/hr 

N2 

n 2 

Nitrogen 

52.0 

C02 

CO 2 

Carbon dioxide 

89.12 

Cl 

ch 4 

Methane 

6033.50 

C2 

c 2 h 6 

Ethane 

297.08 

C3 

c 3 h 8 

Propane 

141.62 

IC4 

3(CH 3 )CH 

Isobutane 

31.84 

NC4 

2(CH 3 )2(CH 2 ) 

n -Butane 

45.11 

IC5 

(2CH 3) CHCH 2 CH 3 

2-Methylbutane 

17.60 

NC5 

CH 3 • 3(CH 2 )CH 3 

n -Pentane 

14.90 

NC6 

CH 3 4 • (CH 2 )CH 3 

n -Hexane 

19.24 


Hot-side heater : outlet at 225 psi. 

Primary flash: adiabatic at 160 psi. 

Secondary flash: adiabatic at 100 psi. 

Pump : outlet at 260 psi. 

Distillation column : use DSTWU with NC4 and IC5 as the light and heavy key 
components, respectively; 99% of the light key is to be in the distillate and 99% 
of the heavy key is to be in the bottoms of the column. A partial condenser is 
required. The top of the column is at 240 psi and the bottom at 255 psi. 

Recycle splitter : 20% of the column bottoms are to be returned to the process. 

Air cooler : outlet at 100°F and 250 psi. 

Gas/liquid heat exchanger : modeled as two heaters with a heat stream used to pass 
heat between them. 

Hot-side heater: outlet at 80°F and 237.7 psi. 

Cold-side heater: outlet at 100 psi. 

A partial setup for this flowsheet exists at Workshops/Twelve-2. It is complete 
except for implementation of the specifications for the gas/gas heat exchanger. 
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Complete this implementation and run the workshop. The results are given at 
Workshops/Twelve-2a. 

Workshop 12.2b Replace all of the heaters in Workshop 12.2a with countercurrent 
heat exchangers, using the shortcut method to implement the design specifications 
given above. Run Aspen Plus to determine the tear streams. Compare the number of 
tear streams to those in Workshop 12.2a. Where necessary, substitute the appropriate 
values obtained from Workshop 12.2a to aid in converging the model. The results are 
given at Workshops/Twelve-2b. 

Workshop 12.2c Convert the heat exchanger design from above to shortcut simula¬ 
tion models and rerun the simulation. The results are given at Workshops/Twelve-2c. 

Workshop 12.26 Design the distillation column using RadFrac at twice the min¬ 
imum reflux ratio and a distillate rate specification such that the column meets the 
specifications above. Select a suitable design to prepare for using it to replace the 
existing DSTWU model. The results are given at Workshops/Twelve-2d. 

Workshop 12.2e Replace the DSTWU model with the RadFrac model from 
Workshop 12.2d. Make the necessary Aspen Plus run(s) required to develop a 
complete material and energy balance. If necessary, substitute the appropriate values 
obtained from Twelve-2a to aid in converging the model. The results are given at the 
Workshops/Twelve- 2e. 

Workshop Notes 

Workshop 12.2a The placement of a pump is appropriate because of the large 
changes in pressure. 

Workshop 12.2b The shortcut design method does not provide any design details but 
merely estimates the area required. It uses a constant overall heat transfer coefficient 
with a default of about 150 Btu/(hr-ft 2 -°F). The best of poor choices would be to obtain 
an estimate of U from a literature source such as Perry and Green (1999). A proper 
design would involve use of the EDR models referred to in Chapter Nine. 

Workshop 12.2c Conversion of the design to a rating model consists merely of 
providing the area calculated in 12.11. Convergence of a small exchanger will proba¬ 
bly require the specification of a smaller area convergence tolerance than the default 
value. 

Workshop 12.2d Using the DSTWU as an approximate solution, the RadFrac model 
solution space was searched with a few trial runs settling on about 35 trays, feed 
onstage 18, and a reflux ratio of about 1.3. For a duplicator block fed a DSTWU (for 
comparison) and a RadFrac block with the same feed, after settling on a DIF value of 
0.641, a sensitivity study varying the reflux ratio was made. Results show that a reflux 
ratio of 1.2 meets the recovery specifications. 
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Workshop 12.2e Caution on defining the overhead product state of the distillation 
column. Try using a partial condenser and observe that the process converges and the 
results are what one might expect. Change to a total condenser and observe that the 
process converges. But is it correct, and does it meet process requirements? 
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CHAPTER THIRTEEN 


OPTIMIZATION 


Aspen Plus’s focus in the area of process optimization is limited primarily to appli¬ 
cations involving process models constructed using Aspen Plus. As examples, process 
optimization can refer to selection of the optimal operating conditions of process equip¬ 
ment, or optimal design of equipment, each subject to the definition of “optimal.” These 
conditions, or design variables, could be driven by “lowest cost,” or “easiest to con¬ 
trol,” or any other basis of evaluation. The applications are typically characterized by a 
function that is to be maximized or minimized subject to variables that are meaningful 
within limits; for example, a negative mole fraction has no meaning within a process. In 
the classical sense, such problems are formally described by an objective function and 
constraints on the independent variables. A formal example of an objective function 
is given by 


<J> = <D(xi,* 2 , (13.1) 

where x\, X2, ..., x n refer to the independent variables. An example of an equality 
constraint, f\ , is given by 


fi = /i(*i, X 2 ,..., x n ) = 0 (13.2) 

and an inequality constraint, / 2 , by 

h = fi(xuX2, ...,x n )<0 (13.3) 

The variables that define the constraints and the objective function need not necessarily 
be continuous but may take on only integer values. The objective function has one 
value for a given set of independent variables. The problem is to find the set of 
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independent variables, subject to the applicable constraints, that yield the optimum 
value of the objective function. In many circumstances there may be several solutions 
to an optimization problem. 

There are many algorithms for solving linear and nonlinear optimization problems. 
Good summaries may be found in Beveridge and Schechter (1970) and in Edgar and 
Himmelblau (2001). Aspen Plus provides the following proprietary methods: 

• SQP: successive quadratic programming 

• Complex: a black-box pattern search 

Each is described in Aspen Plus’s Help. 

13.1 OPTIMIZATION EXAMPLE 

A simple example is the trade-off between the number of theoretical stages and the 
reflux ratio in a distillation column required to obtain a target product composition. 
For this case the independent variables are the number of theoretical stages, a discrete 
variable that takes on integer values, and the reflux ratio, a continuous variable that is 
constrained between a minimum value which corresponds to a practical value near an 
infinite number of theoretical stages, and a maximum which corresponds to a minimum 
number of stages. The objective function has contributions from the following factors. 

• Number of theoretical stages, N 

• Installed cost of the column, where each stage adds an increment of height to the 
column, 1$ (dollars/theoretical stage) 

• Expected life of the column, L (years) 

• Annual maintenance of the column, M$ (dollars/year) 

• Annual reboiler load, R (Btu) 

• Annual condenser load, C (Btu) 

• Cost of heating, H$ (dollars/Btu) 

• Cost of cooling, C$ (dollars/Btu) 

• Operating labor, 0$ (dollars/year) 

A simplified objective function, <D, which represents the total cost of the column 
over its life is defined by 

O = N*I$ + M$*L + R*H$*L + C*C$*L + 0$*L (13.4) 

The implementation of the objective function is carried out using the entry Model 
Analysis, Optimization. The initial input form on which the objective function is defined 
is given in Figure 13.1. In this example, variables defined in the figure are taken from 
the process during each iteration, and are used in calculation of the objective function. 
The various elements that are required for optimization are given in the tabs across 
the top of the figure. Figure 13.2 provides for the formal definition of the objective 
function and its constraints. Note that it is possible to code the objective function in the 
objective function frame in Fortran. Check boxes to select to maximize or minimize 
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^Define -/Objecti 


:ive & Constraints 


</Vary ] Fortran Declarations 



Flowsheet variable 

Definition 


ENSTGS 

Block-Var Block=DSTW Variable=ACT-STAGES Sentence=FiESULTS 


RDUTY 

Block-Var Block=DSTW Variable=REB-DUTY Sentence=RESULTS Units=Btu/hr 


CDUTY 

Block-Var Block=DSTW Variable=COND-DUTY Sentence=RESULTS Units=Btu/hr 

* 




Figure 13.1 Definition of variables. 


^/Define ^/Objective t Constraints ^/Vary ^/Fortran Declarations | 
Objective function 

r Maximize i—-- 

PHI 

(* Minimize 


Constraints associated with the optimization 


Available constraints 


>> 

_< 

<< 


Selected constraints 


Figure 13.2 Definition of objective function. 


the objective function are also available. Alternatively, if the tab Fortran is selected, 
a box appears, for entering the objective function as a series of statements in Fortran, 
as shown in Figure 13.3. Figure 13.4 provides the mechanism to vary an independent 
variable, in this case the reflux ratio, and to define the limits that apply to the variable. 
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Define] Objective & Constraints </Vary ^Fortran Declarations 
Enter executable Fortran statements 


c 

CHGT 

INSTALLED COST $/STAGE 


CHGI 

= 

10000 

c 

TYRS 

= 

LIFE OF PROJECI IN HOI 


TYRS 

= 

10 _ 0*300.*24. 

c 

CS 

= 

COST OF STEAM $/BTU 


CS 

= 

10.0/1000000. 

c 

CC 

= 

COST OF COOLANT $/BTU 


CC 

= 

2 . *CS 

c 

CO 

= 

OPERATING COST $/HR 


CO 

= 

25 „ 

c 

CM 

= 

COST OF MAINTENANACE i 


CM 

= 

.. 0 5 *ENSTCS * CHGI * 10 

c 

ENSIGS 

= NUMBER OF STAGES 


PHI = ENSTGS *CHGT + -CDUTY*CC*TYRS + RDUTY*CS *TYR3 * TYRS*CO + CM 


Figure 13.3 Fortran coding of objective function. 


y Define] Objective & Constraints ^Vary VFortran Declarations 



Disable variable 

Manipulated variable limits 
Lower: 

Upper: 




Figure 13.4 Setup to vary independent variable. 


Having defined the objective function, it is necessary to develop a model of the 
process, which, given the feed, the desired product composition, the reflux ratio, the 
feed location, and the number of theoretical stages, will calculate the reboiler and 
condenser duties. As an example, the DSTWU block will be used to create the model. 
A characteristic of the block is that the number of stages is a continuous variable. If 
a constraint is required, the folder Constraints is selected and variables are selected in 
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^"Define ^Spec Fortran Declarations 


Constraint expressions 
Specification: 


TBOTM 


| Greater than or equal to ^ | 130 


T olerance: 


F 


Vector constraint information 
This is a vector constraint 

First element: pi ~H 

Last element: 


Figure 13.5 Constraint. 

a display similar to Figure 13.1. Although this example does not require a constraint, 
an example of the formulation of a constraint is shown in Figure 13.5. Results of 
the example above are given at Chapter Thirteen Examples/ExampleOne. In this case 
the nature of DSTWU does not permit simultaneous variation of the reflux ratio and 
the number of stages required; however, for the case of varying the reflux ratio only, 
the results show that a tall column with a relatively low reflux ratio produces the lowest 
value of the objective function. 

Although results are available at the Optimization block, additional results at Con¬ 
vergence/Convergence under the tab Iterations show values of the objective function 
and the independent variable through the solution process. 

13.2 WORKSHOPS 

Workshop 13.1 Propane at 14.7 psi and 200°F is to be compressed to 120 psi 
with two single-stage polytropic compressors, using the ASME method, separated by 
a cooler that reduces the first compressor’s outlet temperature to 200°F. Develop an 
optimization that will select the outlet pressure of the first compressor such that the 
work required for the system is minimized. The polytropic efficiency for both com¬ 
pressors is 0.72 and their mechanic efficiency is 1.0. The minimum allowable work 
for both compressors is 10 hp. 

Workshop 13.2a The process shown in Figure 13.6 is to be optimized with respect 
to the pressure of the flash, the column internal reflux ratio, and its distillate to feed 
ratio. In preparation it is necessary to develop a series of process models in steps. 
A key element of the process is the requirement that the column vapor product be 
compressed to the operating pressure of the flash, which for model development is 150 
psi. The flash is to operate adiabatically. It is necessary that a Fortran statement in a 
calculator block be written to accomplish this. Once the code has been written, the 
sequence tab is selected and the appropriate option to be employed upon execution is 
configured. The column has 10 theoretical stages with the feed on stage 5. The column 
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TABLE 13.1 Process Feed 


Component 

Flow Rate (lbmol/hr) 

Propane, C3 

10 

n -Butane, NC4 

20 

n -Pentane, NC5 

20 

n -Hexane, NC6 

10 

Cyclohexane, CC6 

40 


reflux ratio is 1.0 and the distillate/feed ratio is 0.5. The overhead product is saturated 
vapor at a pressure of 50 psi, and the column bottoms pressure is 51 psi. This model 
should not incorporate optimization. The feed to the process is a saturated liquid at 
250 psi and its composition is given in Table 13.1. 

Workshop 13.2b Develop an optimization, maximizing profitability, of the flow¬ 
sheet above by varying the flash pressure between 150 and 250 psi. The objective 
function needs to be coded in Fortran and will incorporate the following terms: 

Stream 5 n -hexane value: $0.1/lbmole 
Stream 5 cyclohexane value: $0.20/mole 
Coolant cost: $1.0 x 10 -6 /Btu 
Heating media cost: $1.0 x 10 -7 /Btu 
Net work cost: $0.020/hp 

Workshop 13.2c Develop an optimization of the flowsheet in three dependent vari¬ 
ables. In addition to flash pressure, the column distillate/feed ratio varies between 0.4 
and 0.7. The internal reflux ratio varies between 0.4 and 0.7 using the same objective 
function as above. 






















WORKSHOPS 197 


Workshop Notes 

Workshop 13.1 A Calculator block is required to pass the outlet pressure of the 
compressor to the pressure specification of the Heater block. Constraints on both com¬ 
pressor horsepower limitations have to be set up and selected on the Objective Function 
and Constraints tab of the Optimization display. 

Workshop 13.2a The Calculator setup is given in Figures 13.7, 13.8, and 13.9. 


^Define ]-^Calculate | ^Sequence Tears Stream Flash EO Options 



Variable name 

Info, flow 

Definition 


PFLASH 

Export 

Block-Var Block=FLASH Variable=PRES Sentence=PARAM Units=psia 


PS 

Import 

Stream-Var Stream=6 Substream=MIXED Variable=PRES Units=psia 

9k 





Figure 13.7 Defining calculator variables. 


Define Calculate Sequence ] Tears Stream Flash EO Options 


Calculation method 
(* Fortran C Excel 

Enter executable Fortran statements 
E FLASH = P6 


Fortran Declarations 


Figure 13.8 Defining calculator calculations. 


Define| Calculate ^Sequence Tears Stream Flash EO Options 


Calculator block execution sequence 

Execute: Block type: Block name: 





List variables as import or export 


Import variables: 

PS 




Export variables: 

PFLASH 



Figure 13.9 Defining calculator sequence. 
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Workshop 13.2b To facilitate the passing of pressure from the compressor to the 
flash, it is necessary to set the flash pressure to zero on its input form. This problem 
was run several times at starting compressor pressures, varying from 250 to 150. In 
all cases the SQP algorithm converged, sometimes with output pass errors, but in all 
cases a value of the final objective function was displayed under the Convergence tab. 
Results showed that a pressure approaching 150 psi yielded the lowest value of the 
objective function. Results of the SQP algorithm need careful checking. 

Workshop 13.2c Although the optimization with three variables gave a con¬ 
verged solution, runs at different starting points exhibited the same behavior as in 
Workshop 13.2b. 
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CHAPTER FOURTEEN 


COMPLEX EQUILIBRIUM 
STAGE SEPARATIONS 


Complex equilibrium stage separations typically involve highly nonideal systems and 
frequently require a literature search to locate vapor-liquid and/or liquid-liquid equi¬ 
librium data when they are not available in Aspen Plus’s data bank. Often data are not 
available in the literature, so that estimation by Unifac is the last resort. In many cases, 
data at more than one set of conditions are involved: for example, azeotrope composi¬ 
tion at two or more pressures. Sometimes it is necessary to employ both vapor-liquid 
and liquid-liquid equilibrium data, and it may be necessary to define unique prop¬ 
erty sets for each column, decanter, and so on. These types of separations involve the 
following general areas: 

• Energy-integrated systems 

• Homogeneous azeotropic distillation 

• Extractive distillation 

• Heterogeneous systems 

The basic ideas and analytical techniques applicable to the applications above are 
described in the following sections. 

14.1 ENERGY INTEGRATION APPLICATIONS 

A common technique is to make use of the energy available in the bottoms product of 
a distillation column. Since the bottoms temperature is almost always higher than that 
of the feed, a bottoms-feed heat exchanger is frequently employed. Modeling of this 
idea is trivial, and a modest amount of energy can be saved. 
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A method that has much greater potential is applicable when two columns operate 
at temperatures such that the overhead temperature of one column is at a higher tem¬ 
perature than the bottoms of the other, and the net reboiler duty of the first is on the 
same order of magnitude as the condenser duty of the other. The procedure is to use 
one side of a single heat exchanger as the reboiler of column 1 and the other side as 
the condenser for column 2. This is illustrated in Figure 14.1. In this example a small 
heat exchanger is used to make up for the insufficiency of the net heat flow between 
heaters 1 and 2. 

One may note that both columns’ configurations have changed considerably. The 
first column is reconfigured as a refluxed stripper, and only one degree of freedom, 
such as the reflux ratio, is available. The bottoms product of the column is split in a 
ratio that is equivalent to the second degree of freedom if the column were not energy 
integrated: for example, the value of the vapor/distillate ratio, or the flow leaving the 
splitter. 

The second column is configured as a reboiled absorber, which also has only one 
degree of freedom. This could be the boil-up ratio or the bottoms product rate. The 
reboiler-condenser is modeled as a pair of heaters with a heat stream flowing from 
the condenser half to the reboiler half. The overhead splitter returns reflux after the 
overhead flow is condensed by the reboiler-condenser heat exchanger. 






Figure 14.1 Shared reboiler-condenser flowsheet. 
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Figure 14.2 Shared reboiler-condenser flowsheet with reference column. 


When implementing this approach it would be prudent first to have developed a 
model of each column without energy integration to facilitate setting up specifications 
for the degrees of freedom. An example of such a reference column, Dist, is illustrated 
in Figure 14.2. This example employs a benzene-toluene distillation as its basis, and 
the demonstration of the energy savings uses two columns each, virtually identical to 
the reference column. Although the original column, Dist, is not part of the scheme, it is 
convenient to have it available for debugging and as an aid to setting the specifications 
for the degrees of freedom. In this example the reboiled absorber operates at a higher 
pressure than the refluxed stripper, to facilitate heat transfer between the bottoms of 
the refluxed stripper and the reboiled absorber’s distillate. Unifac is used merely as a 
convenience to provide estimated vapor-liquid equilibrium and is not recommended 
for a real application. When the model is executed using Aspen Plus’s defaults, it 
does not converge. Initial estimates were provided for both secondary column feeds to 
facilitate convergence. Details are given at Chapter Fourteen Examples/ExampleOne. 
The difference in performance in columns Dist and Dist2 is due to the higher pressure 
in Dist2. The higher pressure reduces the relative volatility of the benzene-toluene 
system. 

The energies used in the reference column (Dist) are as follows: 

• Condenser duty: —2,005,480 Btu/hr 

• Reboiler duty: 2,019,100 Btu/hr 

• Heat transferred through stream Ql: 1,831,371 Btu/hr 

• Reboiler energy use saved (Distl): 91.3% 

• Condenser energy use saved (Dist2): 100% 

The reboiler duty will still need to be provided for the source column, Dist2, and 
the sink column, Distl, will still require that its condenser duty be provided. This type 
of configuration for several columns is analogous to a multieffect evaporator. 
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Figure 14.3 Homogeneous azeotropic distillation flowsheet. 


14.2 HOMOGENEOUS AZEOTROPIC DISTILLATION 

Figure 14.3 shows an example of the configuration of two distillation columns set up for 
homogeneous azeotropic distillation. The intent of the separation is to produce nearly 
pure ethanol from a mixture of ethanol and toluene whose composition is below the 
azeotrope’s composition. The arrangement of the two columns takes advantage of the 
difference in azeotrope composition at different pressures. To facilitate the following 
example, Aspen Plus’s stored Wilson equation binary parameters, which are shown 
in Table 14.1 are used to estimate the composition of the azeotrope as a function of 
pressure. If these were to be used for an industrial application, the results would need 
to be checked against experimental data. 

The basic idea of the separation is explained with reference to Figure 14.3. A 
simple material balance shows that the larger the pressure differential between the two 
azeotropes, the smaller the recycle flow ETAZL. For illustration purposes, the higher- 
pressure column, Bl, is selected to operate at 760 mmHg and the lower-pressure 
column, B2, at 50 mmHg. The stream ETAZH’s composition is approximately 0.815 
mole fraction ethanol and feeds the column B2. The stream ETAZL’s composition 
is approximately 0.685 mole fraction ethanol and is a second feed to column Bl. 
The bottoms product of column B2 is nearly pure ethanol, and the bottoms product 
of column Bl is nearly pure toluene. Other pressures can surely be selected with an 
analysis of the trade-off between recycle flow and associated operating costs against 
the costs associated with purchase and installation of the columns. An initial estimate 
of the design of the two columns can easily be done using the McCabe-Thiele method 
on both columns. The low-pressure column can be analyzed with the toluene as the 
more volatile component. 

When such a design is undertaken, the following procedure is recommended: 

1. Locate a source of experimental data and fit it to an activity coefficient equation. 

2. Determine the azeotrope composition as a function of pressure. 

TABLE 14.1 Ethanol-Toluene Azeotrope at Various Pressures 

Pressure (mmHg) 50 100 400 760 1420 

Mole fraction ethanol 0.685 0.715 0.77 0.815 0.85 
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3. Perform a material balance for the connected columns to estimate the approximate 
flows and compositions of the intercolumn streams. 

4. Design each column rigorously with loose specifications. 

5. Tighten the specifications gradually. 

6. Connect the columns, loosening the specifications. 

7. If convergence is difficult, provide starting estimates for the intermediate streams. 

8. If necessary, change the default convergence method. 

9. When convergence is achieved, tighten the specifications gradually. 

Examples of some of the steps above are given at Examples/etohtolueneazeos, to 
calculate azeotrope composition as a function of pressure, ExampleTwoMB the material 
balance-simple blocks calculation, and ExampleTwo the rigorous solution. 

14.3 EXTRACTIVE DISTILLATION 

Benzene-cyclohexane is an example of a close-boiling system. Kojima et al.’s (1968) 
vapor-liquid equilibrium data for this system are shown in Figure 14.4. Note that the 
vapor-liquid equilibrium curve is almost the same as the y = x line, which makes it 
virtually impossible to separate using standard distillation methods. Systems such as 
these can be separated by introducing a third component, the solvent, with a lower 
volatility than either component, into a distillation column. The solvent must have the 



x - mole fraction benzene 


Figure 14.4 Vapor-liquid equilibrium benzene-cyclohexane at 760 mmHg. 
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ability to change the activity coefficients of the binary system such that the relative 
volatility of the lighter component is increased. Typically, the solvent is introduced 
near the top of the column and is withdrawn at the bottom along with the heavier 
component. The solvent is separated from the heavier component in another column, 
and the solvent is recycled. Such a process is shown in Figure 14.5. 

When designing such a process, in addition to determining the number of theoretical 
stages required and the primary feed location, one must deal with the following issues: 

• Solvent selection 

• Solvent use per mole of feed 

• Stage for introduction of the solvent into the column 

There are no hard and fast rules for solvent selection. A list of potential solvents is 
given by Van Winkle (1967), and general guidelines for solvent selection based on the 
polarities of the components to be separated are given by Wankat (1988). 

The initial step in designing a system is to evaluate a proposed solvent. Does 
it increase the relative volatility of the lighter component? How much solvent is 
required? As an example, Seader and Henley (2006) suggest that aniline would be 
a good solvent for the benzene-cyclohexane system. This can be evaluated by cre¬ 
ating a series of x-y diagrams on a solvent-free basis using the ratio of solvent 
to nonsolvents as a parameter. The calculations to generate such data are given at 
Examples/ExampleThreeData. These data were transferred to the Excel spreadsheet 
Examples/benzenecyclohexaneaniline.xls, where the calculations to produce x -y data 
on an aniline-free basis are given. The results are displayed in Figure 14.6. One may 
observe that as the ratio of aniline to cylohexane plus benzene increases, there is a 
dramatic effect on the relative volatility of the cyclohexane-benzene system. 

Prior to modeling the complete system it is important that each column’s design be 
undertaken on a stand-alone basis. For the first column, initial calculations to estab¬ 
lish the number of stages, reflux ratio, bottom specification, and feed locations are 
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x - liquid mole fraction cyclohexane 

Figure 14.6 Effect of aniline addition on relative volatility in a benzene-cyclohexane 
system. 

incorporated in a sensitivity study to develop the effect of solvent flow rate on prod¬ 
uct compositions. This study is given at Examples/ExampleThreeS. The results of the 
study are shown in Figure 14.7. An optimum occurs at a solvent flow rate of 300 
lbmol/hr. 

The design of the complete system to produce products that are at least 95% pure 
products is given at Examples/ExampleThree2. It is important to note that the amount 
of makeup solvent is difficult to determine in advance; therefore, a series of runs 
ExampleThreel, ExampleThree2, and ExampleThree3 were executed to determine the 
effect of solvent makeup on the composition of the products. These were transferred to 
a spreadsheet, with the results shown in Figure 14.8. An approximate optimum appears 
at a solvent makeup of 0.4 mol/hr per 100 mol of feed. 

14.4 HETEROGENEOUS OPERATIONS 

These types of separations usually involve a decanter in the process. A system such as 
that shown in Figure 14.9, production of purified ethanol at a concentration higher than 
the ethanol-water azeotrope composition, involves the inclusion of a third component, 
some times called an entrainer , which when combined with the azeotrope forms two 
liquid phases, one rich in ethanol and the other rich in water. The water-rich product 
of the decanter is fed to the azeotropic column, Distl, and the entrainer-rich product 
to the alcohol-producing column, Dist2. The aqueous product issues from the bottom 
of Distl and ethanol is the bottom product of Dist2. 
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Figure 14.7 Effect of solvent flow rate on product composition. 
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Figure 14.8 Effect of solvent makeup flow rate on product composition. 



Figure 14.9 Heterogeneous azeotropic distillation flowsheet. 
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The key to solving such problems is a thorough understanding of the liquid-liquid 
phase equilibrium. The binary parameter values that determine the decanter perfor¬ 
mance are critical and must be checked carefully to be sure that the calculations yield 
the result expected. Excellent sources of liquid-liquid equilibrium data are Arlt et al. 
(1979-1987) and Stephen and Stephen (1964). When ternary liquid-liquid equilibrium 
data for a particular component pair is not available, it may be necessary to determine 
its binary parameters from vapor-liquid equilibrium data. 

It is very difficult to merely specify the appropriate blocks in Aspen Plus and expect 
the process to converge. As an example, an analysis of Figure 14.9 shows that two tear 
streams are required. The streams about which one potentially knows the most about 
are the decanter products, streams C and F. It is likely that Aspen Plus’s choice of tear 
streams will need to be changed. The design and simulation of such a system begin 
with the development of a ternary diagram which shows the various tie lines. The use 
of a selected tie line and performance specifications with simple blocks permits the 
development of a material balance. An example is given at Examples/ExampleFour. 
Starting from this point, rigorous blocks are piecewise inserted into the process, and 
each of the columns is designed separately. The complete simulation requires consider¬ 
able refinement, starting with loose specifications and will probably require an increase 
in the number of theoretical stages required and the reflux ratio as the specifications are 
tightened. An exact solution to tighter specifications becomes more and more difficult 
to achieve and should be approached by starting with loose specifications, high reflux 
ratios, and excessive number of stages. 

14.5 WORKSHOPS 

Workshop 14.1 A saturated liquid feed of 50 mol % methylcyclohexane, the remain¬ 
der being n -heptane, is to be separated by extractive distillation using phenol as the 
extracting agent. It is required that the methylcyclohexane product be at least 95 mol 
% pure and the toluene product be at least 95 mol % pure. The phenol recycle stream 
is to be at least 99.5 mol % pure. The columns are to operate at 760 mmHg. The 
binary parameter values stored in the Aspen Plus data bank are adequate. Design the 
system. 

Workshop 14.1a Determine the optimum amount of phenol to be used as an auxil¬ 
iary feed to the first column, Distl. Determine the number of theoretical stages required 
and the location of the feeds. 

Workshop 14.1b Place the second column, Dist2, on the process flowsheet of Work¬ 
shop 14.1a. Use the bottoms product of Distl as the feed to Dist2. Design Dist2. 

Workshop 14.1c Connect the bottoms of Dist2 to Distl. Be sure to provide a 
makeup feed stream to account for the phenol losses from each column. Execute the 
completed model calculations until convergence is achieved. 

Workshop 14.2 A mixture containing 28.1 mol % methanol, 70.9 mol % water, and 
1 mol % mesity oxide is to be distilled in a column designed to produce a distillate 
with at least 99.5 mol % methanol and a bottoms with at least 99.5 mol % water. 
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Figure 14.10 Methanol-water distillation with mesityl oxide sidestream. 


The feed is available at its bubble point. A sidestream whose purpose is to remove 
mesity oxide is required. The sidestream will be fed to a decanter with an auxiliary 
feed of water. Two phases will be formed, one mesity oxide-rich, the other primarily 
water. The mesity oxide phase will be removed and the aqueous phase returned to the 
distillation column. A sketch of the process is shown in Figure 14.10. The following 
is known about the ternary system: No ternary liquid-liquid equilibrium data exist 
but a heterogeneous azeotrope has been observed; no vapor-liquid equilibrium for the 
binary mesity oxide-methanol exists. Design the process concepts starting with the 
steps given in Workshops 14.2a through 14.2h. 

Workshop 14.2a A mesity oxide-water azeotrope exits at 760 mmHg and 91.8°C. 

MSO in the vapor phase 65.2 wt% 

MSO in the liquid phases 93.9 and 2.0 wt% 

Determine the van Laar and Uniquac equations binary parameters. 

Workshop 14.2b Estimate the vapor-liquid equilibrium of the methanol-mesity 
oxide system at one 760 mmHg with a property analysis run using Unifac. 

Workshop 14.2c Regress the results of Workshop 14.2b to the a L:J parameters of 
the van Laar and Uniquac equations. 

Workshop 14.2d Fit the vapor-liquid equilibrium data for the methanol-water 
system of Dunlop (1948) given in Table 14.2 to the van Laar and Uniquac equations 
using the a^j parameters. 

Workshop 14.2e Develop the liquid-liquid equilibrium tie lines for the ternary 
system and create a triangular ternary equilibrium diagram at 90°C using the van Laar 
equation. 
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TABLE 14.2 Methanol (1)-Water (2) Vapor-Liquid Equilibrium 

Temperature (°C) Pressure (mmHg) Vapor Mole Fraction, y Liquid Mole Fraction, x 


96.4 

760 

0.02 

0.134 

93.5 

760 

0.04 

0.23 

91.2 

760 

0.06 

0.304 

89.3 

760 

0.08 

0.365 

87.7 

760 

0.1 

0.418 

84.4 

760 

0.15 

0.517 

81.7 

760 

0.2 

0.579 

78 

760 

0.3 

0.665 

75.3 

760 

0.4 

0.729 

73.1 

760 

0.5 

0.779 

71.2 

760 

0.6 

0.825 

69.3 

760 

0.7 

0.87 

67.5 

760 

0.8 

0.915 

66 

760 

0.9 

0.958 

65 

760 

0.95 

0.975 


Workshop 14.2f Make preliminary runs of the column to aid in locating the feed 
and sidestream. Experiment with product and sidestream flows to develop appropriate 
distillation column specifications. 

Workshop 14.2g Link the column and the decanter and develop a solution with 
loose specifications. 

Workshop 14.2h Refine the specification to develop a final conceptual design. 

Workshop Notes 

Workshop 14.2a Aspen Plus provides a unique form of binary parameter equation 
for almost all activity coefficient equations, which is given in 

In Aij = Oij + -jr + Cij In T + dijT + (14.1) 

When regressing data, one usually selects ci^j or bij or both. The other parameters 
are rarely used. For a case where only one data point is available b t j may be used, 
but extrapolation to other temperatures may cause difficulties. 

Workshop 14.2b A plot of vapor-liquid equilibrium for mesityl oxide-water based 
on Unifac estimates is shown in Figure 14.11. 

Workshop 14.2e After performing the property analysis, a triangular plot is avail¬ 
able from Aspen Plus’s Plot Wizard (see Figure 14.12). Note that a methanol-water 
distillation, without mesity oxide, would take place along the left side of the triangle; 
however, the presence of a small quantity of mesity oxide places the composition to 





- mole fraction vapor methanol 
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Ternary map for MSO/METHANOL/WATER 



Figure 14.12 Ternary diagram methanol-water-mesityl oxide at 760 mmHg. 
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the right. The strategy employed here is to locate a vapor composition, within the 
distillation column, which is to be withdrawn as a sidestream, to the right of a tangent 
drawn from the left base of the phase envelope along the envelope (see the heavy 
line on Figure 14.10). When water is added to such a location, the composition of 
a condensed vapor is within the envelope and forms two phases, one being mesity 
oxide-rich, which is removed. The aqueous phase is recycled and mixed with the feed 
to the column. 

Workshop 14.2f Many trial runs are required to estimate the number of stages, the 
locations of the feed and sidestreams, the value of the reflux ratio, and flow specifica¬ 
tions for a product and the sidestream. A good look at the composition profile of the 
column vapors should provide an estimate of the sidestream location. Use of sensitivity 
studies will aid in establishing column specifications. 

Workshop 14.2g When linking the column and the decanter, one may experience 
convergence failure. This may result due to the fact that the flow rate of the recycle 
stream will vary during the iterations. This is an important factor to take into account 
when establishing or altering the column specifications to accommodate the recycle 
stream. For example, a total distillate flow specification may need to be changed to a 
distillate/feed specification. Occasionally, a block within a process will fail to converge 
during the iterations but is converged at the conclusion of a run. This is not a problem 
as long as all blocks and recycles are converged satisfactorily at the conclusion of the 
run. 

Workshop 14.2h Reducing the reflux ratio and making small adjustments to other 
column details, such as the number of stages, to the result of Workshop 14.2g will 
produce a completed conceptual design. 
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